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1 - Introduction 



The Cauchy problem for a system of conservation laws in one space dimension takes the form 

u t + f(u) x = 0, (1.1) 

u(0,x)=u(x). (1.2) 

Here u = {u\ , . . . , u n ) is the vector of conserved quantities, while the components of / = (/i, . . . , f n ) 
are the fluxes. We assume that the flux function / : TR n h- M n is smooth and that the system 
is strictly hyperbolic, i. e., at each point u the Jacobian matrix A(u) = Df(u) has n real, distinct 
eigenvalues 

Ai («)<■■■< An(«). (1.3) 
One can then select bases of right and left eigenvectors Vi(u), k(u), normalized so that 

N-l Irl-fj % (1.4) 

Several fundamental laws of physics take the form of a conservation equation. For the relevance 
of hyperbolic conservation laws in continuum physics we refer to the recent book of Dafermos [D]. 

A distinguished feature of nonlinear hyperbolic systems is the possible loss of regularity. Even 
with smooth initial data, it is well known that the solution can develop shocks within finite time. 
Therefore, global solutions can only be constructed within a space of discontinuous functions. The 
equation (1.1) must then be interpreted in distributional sense. A vector valued function u = u(t, x) 
is a weak solution of (1.1) if 

j j [u</h + f(u) <f> x ] dxdt = (1.5) 

for every test function cj) 6 C\, continuously differentiable with compact support. When disconti- 
nuities are present, weak solutions may not be unique. To single out a unique "good" solution of 
the Cauchy problem, additional entropy conditions must be imposed along shocks [Lx], [LI]. These 
are often motivated by physical considerations [D]. 

Toward a rigorous mathematical analysis of solutions, the lack of regularity has always been 
a considerable source of difficulties. For discontinuous solutions, most of the standard tools of 
differential calculus do not apply. Moreover, for general n x n systems, the powerful techniques of 
functional analysis cannot be used. In particular, solutions cannot be obtained as fixed points of a 
nonlinear transformation, or in variational form as critical points of a suitable functional. Dealing 
with vector valued functions, comparison arguments based on upper and lower solutions do not 
apply either. Up to now, the theory of conservation laws has thus progressed largely by developing 
ad hoc methods. In particular, a basic building block is the so-called Riemann problem, where the 
initial data is piecewise constant with a single jump at the origin: 



u(0, x) 



u if x < , 
u + if x > . 



Weak solutions to the Cauchy problem (1.1)-(1.2) were constructed in the celebrated paper of 
Glimm [G]. This global existence result is valid for small BV initial data and under the additional 
assumption 

(H) For each i G {1, . . . ,n}, the z-th characteristic field is either linearly degenerate, so that 

D\i(u) ■ n(u) = for all u, (1.6) 
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or else it is genuinely nonlinear, i. e. 



D\i(u)-ri(u) > 



for all u . 



(1.7) 



In [G], an approximate solution of the general Cauchy problem is obtained by piecing together 
solutions of several Riemann problems, with a restarting procedure based on random sampling. 
The key step in Glimm's proof is an a priori estimate on the total variation of the approximate 
solutions, obtained by introducing a wave interaction potential. In turn, the control of the total 
variation yields the compactness of the family of approximate solutions, and hence the existence 
of a strongly convergent subsequence. Alternative constructions of approximate solutions, based 
on front tracking approximations, were subsequently developed in [DPI], [B2], [Ri], [BaJ]. 

The above existence results are all based on a compactness argument which, by itself, does not 
guarantee the uniqueness of solutions. The well posedness of the Cauchy problem has now been 
established in a series of papers [B3], [BC1], [BCP], [BLY] , [BLF], [BG], [BLe]. The main results 
can be summarized as follows: 

- The solutions obtained as limits of Glimm or front tracking approximations are unique and 
depend Lipschitz continuously on the initial data, in the L 1 norm. 

- Every small BV solution of the Cauchy problem (1.1)-(1.2) which satisfies the Lax entropy 
conditions coincides with the unique limit of front tracking approximations. 

For comprehensive account of the recent uniqueness and stability theory we refer to [B5]. 

A long standing conjecture is that the entropic solutions of the hyperbolic system (1.1) actually 
coincide with the limits of solutions to the parabolic system 



letting the viscosity coefficient e — > 0. In view of the recent uniqueness results, it looks indeed 
very plausible that the vanishing viscosity limit should single out the unique "good" solution of 
the Cauchy problem, satisfying the appropriate entropy conditions. In earlier literature, results in 
this direction were based on three main techniques: 

1 - Comparison principles for parabolic equations. For a scalar conservation law, the 
existence, uniqueness and global stability of vanishing viscosity solutions has been established in a 
famous paper by Kruzhkov [K]. The result is valid within the more general class of L°° solutions, 
also in several space dimensions. For an alternative approach based on nonlinear semigroup theory, 
see also [Cr]. 

2 - Singular perturbations. Let u be a piecewise smooth solution of the n x n system (1.1), 
with finitely many non-interacting, entropy admissible shocks. In this special case, using a singular 
perturbation technique, Goodman and Xin [GX] were able to construct a sequence of solutions u £ 
to (1.8) e , with u £ — > u as e — > 0. See also [Yu] for further results in this direction. 

3 - Compensated compactness. If, instead of a BV bound, only a uniform bound on the L°° 
norm of solutions of (1.8) e is available, one can still construct a weakly convergent subsequence 
u £ u. In general, we cannot expect that this weak limit satisfies the nonlinear equations (1.5). 
However, for a class of 2 x 2 systems, in [DP2] DiPerna showed that this limit u is indeed a 



ut + f(u) x = 



eu. 



'XX ) 



(1-8). 
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weak solution of (1.1). The proof relies on a compensated compactness argument, based on the 
representation of the weak limit in terms of Young measures, which must reduce to a Dirac mass 
due to the presence of a large family of entropies. We remark that the solution is here found in 
the space L°°. Since the known uniqueness results apply only to BV solutions, the uniqueness of 
solutions obtained by the compensated compactness method remains a difficult open problem. 

In our point of view, to develop a satisfactory theory of vanishing viscosity limits, the heart 
of the matter is to establish a priori BV bounds on solutions u(t, •) of (1.8) e , uniformly valid for 
all t G [0, oo [ and e > 0. This is indeed what we will accomplish in the present paper. Our 
results apply, more generally, to strictly hyperbolic nxn systems with viscosity, not necessarily in 
conservation form: 

u t + A{u)u x = e u xx . (l-9) e 

As a preliminary, we observe that the rescaling of coordinates s = t/e, y = x/e transforms the 
Cauchy problem (1.9) e , (1.2) into 

u s + A(u)u y = Uyy, u(0 , y) = u e (y) = u(ey) . 

Clearly, the total variation of the initial data u £ does not change with e. To obtain a priori BV 
bounds and stability estimates for solutions of (1.9) £ , it thus suffices to consider the system 

u t + A(u)u x = u xx , (1-10) 

and derive estimates uniformly valid for all times t > , depending only on the total variation of 
the initial data u. 

The first step in our proof is a decomposition of the gradient u x = Yl scalar com- 

ponents. In the purely hyperbolic case without viscosity, it is natural to decompose u x along a 
basis {ri, . . . , r n } of eigenvectors of the matrix A(u). Remarkably, this choice does not work here. 
Instead, we will decompose u x as a sum of gradients of viscous travelling waves, selected by a 
center manifold technique. 

As a second step, we study the evolution of each component Vi, which is governed by a scalar 
conservation law with a source term, accounting for nonlinear wave interactions. Uniform bounds 
on these source terms are achieved by means of a transversal interaction functional, controlling the 
interaction between waves of different families, and suitable swept area and curve length functionals, 
controlling the interaction of waves of the same family. All these can be regarded as "viscous" 
counterparts of the wave interaction potential, introduced by Glimm [G] in the purely hyperbolic 
case. Finally, on regions where the diffusion is dominant, the strength of the source term is 
bounded by an energy functional. All together, these estimates yield the desired a priori bound on 
||m x (£, OHl 1 ' independent of t G [0, oof . 

Similar techniques can also be applied to a solution z = z(t, x) of the variational equation 

z t + [DA(u) ■ z]u x + A(u)z x = z xx , (1-11) 

which describes the evolution of a first order perturbation to a solution u of (1.10). Assuming that 
the total variation of u remains small, we shall establish an estimate of the form 

HV)|| L1 <L \\z(0, -)|| L1 for all t>0, (1.12) 
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valid for all solutions of (1.11). As soon as this estimate is proved, as in [Bl] a standard homotopy 
argument yields the Lipschitz continuity of the flow of (1.10) w.r.t. the initial data, uniformly in 
time. 

By the simple rescaling of coordinates t i-> et, x i— > ex, all of the above estimates remain 
valid for solutions u £ of the system (1.9) e . By a compactness argument, these BV bounds imply 
the existence of a strong limit u £m — > u in L^, at least for some subsequence e m — > 0. In the 
conservative case where A = Df, it is now easy to show that this limit u provides a weak solution 
to the Cauchy problem (1.1)-(1.2). 

At this intermediate stage of the analysis, since we are using a compactness argument, it is 
not yet clear whether the vanishing viscosity limit is unique. In principle, different subsequences 
e m — > may yield different limits. Toward a uniqueness result, in [B3] the second author introduced 
a definition of viscosity solution for the hyperbolic system of conservation laws (1.1), based on local 
integral estimates. Roughly speaking, a function u is a viscosity solution if 

• In a forward neighborhood of each point of jump, the function u is well approximated by the 
self-similar solution of the corresponding Riemann problem. 

• On a region where its total variation is small, u can be accurately approximated by the solution 
of a linear system with constant coefficients. 

For a strictly hyperbolic system of conservation laws satisfying the standard assumptions (H), 
the analysis in [B3] proved that the viscosity solution of a Cauchy problem is unique, and coincides 
with the limit of Glimm and front tracking approximations. The definition given in [B3] was 
motivated by a natural conjecture. Namely, the viscosity solutions (characterized in terms of local 
integral estimates) should coincide precisely with the limits of vanishing viscosity approximations. 

In the present paper we adopt an entirely similar definition of viscosity solutions and prove 
that the above conjecture is indeed true. Our results apply to the more general case of (possibly 
non-conservative) quasilinear strictly hyperbolic systems. In particular, we obtain the uniqueness 
of the vanishing viscosity limit. 

As in [B3] , [BLFP] , the underlying idea is that a semigroup is entirely determined by its local 
behavior on piecewise constant initial data. Namely, if two semigroups yield the same solution 
to each Riemann problem, then they coincide. In our proof of uniqueness, a basic step is thus 
the analysis of the vanishing viscosity solution to a general Riemann problem. The construction 
given here extends the previous results by Lax and by Liu to general, non-conservative hyperbolic 
systems. As in the cases considered in [Lx], [LI], for a given left state u~ there exists a Lipschitz 
continuous curve of right states u + which can be connected to u~ by z-waves. These right states are 
here obtained by looking at the fixed point of a suitable contractive transformation. Remarkably, 
our center manifold plays again a key role, in defining this transformation. 

Our main results are as follows. 

Theorem 1. Consider the Cauchy problem for the hyperbolic system with viscosity 

u t + A{u)u x = e u xx u(0,x) = u(x). (1-13) £ 

Assume that the matrices A{u) are strictly hyperbolic, smoothly depending on u in a neighborhood 
of a compact set K C M n . Then there exist constants C,L,L' and 5 > such that the following 
holds. If 

Tot.Var.{u} <5, lim u(x) € K , (1.14) 

X — > — oo 
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then for each e > the Cauchy problem (1.13) £ has a unique solution u s , defined for all t > 0. 
Adopting a semigroup notation, this will be written as t i— > u 6 (t, •) = Sfu. In addition, one has: 



BV bounds : Tot. Var.{S £ t u] < C Tot. Var.{u} . (1.15) 



L 1 stability : \\Sfu - Sfv\\ L1 < L \\u - u|| Ll , (1.16) 

p^u- S £ s u\\ hl < L' (\t - s\ + \Vei - ■ (1.17) 

Convergence: ^4s e — » 0+, £/ie solutions u E converge to the trajectories of a semigroup S such 
that 

\\S t u — 5 s v|| x < L ||u — tj|| L i + L' \t — s\ . (1-18) 

These vanishing viscosity limits can be regarded as the unique vanishing viscosity solutions of 

the hyperbolic Cauchy problem 

u t + A(u)u x = 0, u(0, x) = u(x) . (1-19) 

In the conservative case A{u) = Df(u), every vanishing viscosity solution is a weak solution 

of 

u t + f(u) x = 0, u(0, x) = u{x) , (1.20) 

satisfying the Liu admissibility conditions. 

Assuming, in addition, that each field is genuinely nonlinear or linearly degenerate, the vanish- 
ing viscosity solutions coincide with the unique limits of Glimm and front tracking approximations. 

Notice that in the above theorem the only key assumptions are the strict hyperbolicity of the 
system and the small total variation of the initial data. It is interesting to compare this result with 
previous literature. 

1. Concerning the global existence of weak solutions, Glimm's proof requires the additional as- 
sumption (H) of genuine nonlinearity or linear degeneracy of each characteristic field. This as- 
sumption has been greatly relaxed in subsequent works by Liu [L2] and Liu and Yang [LY], but 
never entirely removed. The underlying technical reason is the following. In all papers based on 
the Glimm scheme (or front tracking), the construction of approximate solutions as well as the BV 
estimates rely on a careful analysis of the Riemann problem. In this connection, the hypothesis 
(H) is a simplifying assumption, which guarantees that every Riemann problem can be solved in 
terms of n elementary waves (shocks, centered rarefactions or contact discontinuities), one for for 
each characteristic field i = 1, . . . ,n. At the price of considerable technicalities, this assumption 
can be replaced by some other condition, implying that all solutions of the Riemann problem can 
obtained by piecing together a finite (but possibly large) number of elementary waves. 
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On the other hand, our present approach is based on vanishing viscosity limits and does not 
make any reference to Riemann problems. Global existence is obtained for the whole class of 
strictly hyperbolic systems. 

2. Concerning the uniform stability of entropy weak solutions, the results previously available 
for n x n hyperbolic systems [BC1], [BCP], [BLY] always required the assumption (H). For 2x2 
systems, this condition was somewhat relaxed in [AM]. Again, we remark that the present result 
makes no reference to the assumption (H). 

3. For the viscous system (1.10), previous results in [L3], [SX], [SZ], [Yu] have established the 
stability of special types of solutions, such as travelling viscous shocks or viscous rarefactions, 
w.r.t. suitably small perturbations. Taking e = 1, our present theorem yields at once the uniform 
Lipschitz stability of all viscous solutions with sufficiently small total variation, w.r.t. the L 1 
distance. 

Remark 1.1. It remains an important open problem to establish the convergence of vanishing 
viscosity approximations of the form 

u t + A(u)u x =e(B{u)u x ) x (1.21) e 

for more general viscosity matrices B. In the present paper we are exclusively concerned with the 
case where B is the identity matrix. For systems which are not in conservative form, we expect 
that the limit of solutions of (1.21) e , as e — > 0, will be heavily dependent on the choice of the 
matrix B. 

The plan of the paper is as follows. Section 2 collects those estimates which can be obtained by 
standard parabolic techniques. In particular, we show that the solution of (1.10) with initial data 
u 6 BV is well defined on an initial time interval [0, t] where the L°° norms of all derivatives decay 
rapidly. Moreover, for large times, as soon as an estimate on the total variation is available, one 
immediately obtains a bound on the L 1 norms of all higher order derivatives. Our basic strategy 
for obtaining the BV estimate is outlined in Section 3. The decomposition of u x as a sum of 
gradients of viscous travelling profiles is performed in Section 5. This decomposition will depend 
pointwise on the second order jet (u x ,u xx ), involving 2n scalar parameters. To fit these data, we 
must first select n smooth families of viscous travelling waves, each depending on 2 parameters. 
This preliminary construction is achieved in Section 4, relying on the center manifold theorem. In 
Section 6 we derive the evolution equation satisfied by the gradient components and analize the 
form of the various source terms. As in [G], our point of view is that these source terms are the 
result of interactions between viscous waves, and can thus be controlled by suitable interaction 
Junctionals. In Sections 7 to 9 we introduce various Lyapounov functionals, which eventually allow 
us to estimate the integral of all source terms. The proof of the uniform BV bounds is then 
completed in Section 10. 

In Section 11 we study the linearized evolution equation (1.11) for an infinitesimal perturbation 
z, and derive the key estimate (1.12). In turn, this yields the Lipschitz continuity of the flow, 
stated in (1.16). Some of the estimates here require lengthy calculations, which are postponed to 
the Appendices. Section 12 contains an additional estimate for solutions of (1.11), showing that, 
even in the parabolic case, the bulk of a perturbation propagates at a finite speed. This estimate 
is crucial because, passing to the limit er — >■ 0, it implies that the values of a vanishing viscosity 
solution u(t, •) on an interval [a, b] depend only on the values of the initial data u(0, •) on a bounded 
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interval [a — (3t, b + (3t\. In Section 13 we study the existence and various properties of a semigroup 
obtained as vanishing viscosity limit: S = limS £m . At this stage, we only know that the limit 
exists for a suitable subsequence e m — ► 0. In the case of a system of conservation laws satisfying the 
standard assumptions (H), we can show that every limit solution satisfies the Lax shock conditions 
and the tame oscillation property. Hence, by the uniqueness theorem in [BG], the limit is unique 
and does not depend on the subsequence {e m }. This already achieves a proof of Theorem 1 valid 
for this special case. 

Toward a proof of uniqueness in the general case, in Section 14 we construct a self-similar 
solution co(t, x) = uj(x/t) to the non-conservative Riemann problem, and show that it provides the 
unique vanishing viscosity limit. A definition of viscosity solution in terms of local integral estimates 
is introduced in Section 15. By a minor modification of the arguments in [B3], [B5] we prove that 
these viscosity solutions are unique and coincide with the trajectories of any semigroup S = lim S Sm 
obtained as limit of vanishing viscosity approximations. Since this result is independent of the 
subsequence {e m }, we obtain the convergence to a unique limit of the whole family of viscous 
approximations Sfu — > S t u, over all real values of e. This completes the proof of Theorem 1. 

Finally, in Section 16 we derive two easy estimates. One is concerned with the dependence 
of the the limit semigroup S on the coefficients of the matrix A in (1.19). The other estimate 
describes the asymptotic limit of solutions of the parabolic system (1.10) as t — > oo. 



2 - Parabolic estimates 

In classical textbooks, the local existence and regularity of solutions to the parabolic system 
(1.10) are derived by regarding the hyperbolic term A(u) first order perturbation of the 

heat equation. This leads to the definition of mild solutions, characterized by the representation 

u(t) = G(t) * u(0) + f G(t- s)*A(u(s))u x (s)ds 
Jo 

in terms of convolutions with the standard heat kernel G. 

In this initial section we collect all the relevant estimates which can be achieved by this 
approach. In particular, we prove various decay and regularity results for solutions of (1.10) as 
well as (1.11). Given a BV solution u = u(t,x) of (1.10), consider the state 

u* = lim u(t,x) , (2.1) 

X— > — oo 

which is clearly independent of time. We then define the matrix A* = A{u*) and let A*, r*, I* 
be the corresponding eigenvalues and right and left eigenvectors, normalized as in (1.4). It will be 
convenient to use "•" to denote a directional derivative, so that z • A(u) = DA(u) ■ z indicates the 
derivative of the matrix valued function u i-> A(u) in the direction of the vector z. We can now 
rewrite the systems (1.10) and (1.11) respectively as 

ut + A*u x -u xx = {A* -A(u))u x , (2.2) 
zt + A*z x - z xx = (A* - A(u))z x - (z» A(u))u x . (2.3) 

In both cases, we regard the right hand side as a perturbation of the linear parabolic system with 
constant coefficients 

w t + A*w x -w xx = 0- (2.4) 
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We denote by G* the Green kernel for (2.4), so that 



w(t, x) = j G* (t, x-y) w(0, y) dy . 



The matrix valued function G* is easily computed. Indeed, if w solves (2.4), then its i-th component 
Wi = I* ■ w satisfies the scalar equation 

Wi, t + Km, x - m,xx = . 
Therefore Wi(t) = G*(t) * Wi(0), where 

G*{t,x) = -^=exp{- K 



2^t 1 4t 

Looking at the explicit form of its components, it is clear that the Green kernel G* = G*(t,x) 
satisfies the bounds 

||G*(i)|| Ll < k, \\Gl{t)\\^ < A , ||GL(t)|| L1 < - t , (2.5) 

for some constant k and all t > 0. It is important to observe that, if u is a solution of (2.2), then 
z = u x is a particular solution of the variational equation (2.3). Hence all the estimates proved 
for 

^X 5 ^XX 

certainly valid also for the corresponding derivatives u xx , u xxx . Assuming that the 
initial data u(0, •) has small total variation, we now derive some estimates on higher derivatives. 
In particular, we will show that 

• The solution is well defined on some initial interval [0, t\, where the L°° norm of all derivatives 
decays rapidly. 

• As long as the total variation remains small, the solution can be prolonged in time. In this 
case, all higher order derivatives remain small. Indeed, waiting a long enough time, one has 
||«xxx(*)|| L i << ||«xx(*)|| L i << ||«x(*)|| L i = Tot.Var.{u(t)}. 

Proposition 2.1. Let u, z be solutions of the systems (2.2)-(2.3), satisfying the bounds 

K(*)|| L i <*o, ||z(t)|| Ll <5 , (2-6) 

for some constant Sq < 1 and all t £ [0,t], where 

and k is the constant in (2.5). Then for t € [0,t] the following estimates hold: 

2k5q 



:(t)|| Ll , |M«)|| L i<-^. ( 2 -8) 



, II^WUli < — — , (2.9) 
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>(*)|| L -, IM*)|| L .<^- (2-10) 



Proof. The function z x can be represented as 

z x (t) = G* x (t) * z(0) + J* G* x (t - s) *{{A(u) - A*)z x (s) - (z . A{u))u x (s)) ds . (2.11) 
Using (2.5) and (2.6) we obtain 

f G* x (t - s) * \(A* - A(u))z x (s) - {z . A(u))u x (s)\ds 

J L 1 

<j o M^^c* — - { II^^C^)!!^^ Il^^lli.— ll^c-s)!!^ -i- II^C^)!!^— Il^^-lli.— M-"-rC-»)ll^ } 

< 28 k\\DA\\ ■ -== 11^(8)11, da. 
Jo V * — s 

Consider first the case of smooth initial data. We shall argue by contradiction. Assume that there 
exists a first time r < t such that the equality in (2.8) holds. Then, observing that 

f 1 1 
: ds = / — : da = tt < 4 



we compute 



y/s(t — s) Jo y/cr(l - a) 



Mr)|| L1 < -^=8o + 2k5 \\DA\\ ■ [ 
v r Jo 



1 25 k , 
ds 



y/T — S y^S 



K&0 9 ^oii^.ii 2K(5n 

<-£ + 16k 2 k a S%\\DA\\ < — ^, 
v r V r 

reaching a contradiction. Hence, (2.8) is satisfied as a strict inequality for all t G [0,t]. Observing 
that this estimate depends only on the L 1 norms of u x and z, by an approximation argument we 
obtain the same bound for general initial data, not necessarily smooth. Since z = u x is a particular 
solution of (2.3), the bounds (2.8) certainly apply also to z x = u xx . 

A similar technique is used to establish (2.9). Indeed, we can write 
z xx (t) = G* x (t/2)*z x (t/2)- [ G* x (t-s)\(z.A(u))u x (s) + (A(u)-A*)z x (s)} ds . (2.12) 

Jt/2 k Jx 



We will prove (2.9) first in the case z xx = u xxx , then in the general case. If (2.9) is satisfied as an 
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equality at a first time r <i, using (2.12) and recalling the definitions (2.7) we compute 
||^x(r)|| Ll < -j= ■ + J ^ ~y== ■ j||**»^(«K(s)|| L i + \\z» (uj.* A(ti))« a .(a)|| Ll 

+ \\z» A(u)u xx {s)\\ hl + \\u x » A(u)z x (s)\\ L1 + (A(u) - A*)z xx {s) ^ ds 

<^r + f ^^•{5oII^I|l»||^( s )|| l1 +<5oI| j d 2 ^I|l»||^( s )||J j1 

+ S \\DA\\^\\ u xxx {s)\\ L1 + 6 \\DA\\^\\ 

z xx( s ) || L i 

+ 6o\\DA\\^\\ 

2 xx( s )|| L i 

<^^ + K < 5 (4« 2 <5g|| J D 2 ^|| L »+20K 2 < y || J DA|| L ») T * da, 

T J T /2 S^T - S 

< + 20k 3 k a <5 2 • — == < , 

r y/r/2 r 

reaching a contradiction. 

Finally, using (2.12) and (2.8)-(2.9), the bounds in (2.10) are proved by the estimate 

||^*(r)|| Loo < —}L= • + ! /I — • j \\z x • A(u)u x (s)\\ Loa + ||z • (u x • A(tt))tt x (s)|| 

y/r/2 t/z j T /2 vt - s ^ 

+ || z • A(u)u xx (s) || LOO + • A(u)z a; (s)|| LOO + ||(j4(u) -^o)^x(s)|| loo |^ 

^ |jD2j4||+4 |jDA r i ds 



< — — h 46k kao • — j- < 



° r/2 t^T ' 

□ 

Corollary 2.2. In the same setting as Proposition 5.1, assume that the bounds (2.6) hold on a 
larger interval [0,T]. Then for all t 6 [t, T] there holds 

K*(*)|| L i> K(t)|| L o=, ||^(*)|| l1 = (2.13) 

K**(t)|| L i> Kx(*)|| L o=> = o(i)-<5g, (2.i4) 

^IIl- |M*)|U = 0(1)-S*. (2.15) 



?/, 



Proof. It suffices to apply Proposition 2.1 on the interval [t — t, t\. □ 

Proposition 2.3. Let u = u(t,x), z = z(t,x) be solutions of (2.2), (2.3) respectively, such that 

Tot.Var.{u(0 r )}<^, |K0)|| Ll <^. (2.16) 
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Then u,z are well defined on the whole interval [0,t\ in (2.7), and satisfy 

<5o II /.Ml ^ #0 



U 



:WL<T' fe M- (2 - 17) 



Proof. We have the identity 

z(t) = G*(t)z(0)- J G*(t-s)(z»A(u)u x (s)-(A(u)-A )z x (s)^ds. (2.18) 

As before, we first establish the result for z = u x , then for a general solution z of (2.3). Assume 
that there exists a first time r < t where the bound in (2.17) is satisfied as an equality. Estimating 
the right hand side of (2.18) by means of (2.5) and (2.8), we obtain 

\\z(t) t1 <-^-+ PA|| L » 



4k Jo v 7 * 
< ^ +4/c/c v i5ov / t : < y , 

reaching a contradiction. D 

To simplify the proofs, in all previous results we used the same hypotheses on the functions 
u x and z. However, observing that z solves a linear homogeneous equation, similar estimates can 
be immediately derived without any restriction on the initial size ||z(0)|| L1 . In particular, from 
Proposition 2.3 it follows 

Corollary 2.4. Let u = u(t,x), z = z(t,x) be solutions of (2.2), (2.3) respectively, such that 
||« X (0)|| L1 < <5o/4k. Then u,z are well defined on the whole interval [0,t\ in (2.7), and satisfy 

K(*)|| L i <2/e||ii x (0)|| Ll) |K*)|| L i <2k||*(0)|| l1 , t€[0,t\. (2.19) 

A summary of the main estimates is illustrated in fig. 1. On the initial interval t € [0, t\, with 
1/Sq we have 

K(*)|| L i <^o, (2.20) 
while the norms of the higher derivatives decay: 

IK*|| L i = 0(1) • 5 /Vi, \\u xxx \\ hl = 0(1) • S /t. 

On the other hand, for t > i, as long as (2.20) remains valid we also have 

|K*|| L i = 0{l) ■ Si , ||«*xx|| L i = o(i) • <&o • 

These bounds (the solid lines in fig. 1) were obtained in the present section by standard parabolic- 
type estimates. The most difficult part of the proof is to obtain the estimate (2.20) for large times 
t G [t, oo [ (the broken line in fig. 1). This will require hyperbolic-type estimates, based on the 
local decomposition of the gradient u x as a sum of travelling waves, and on a careful analysis of 
all interaction terms. 
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figure 1 



3 - Outline of the BV estimates 

It is our aim to derive global a priori bounds on the total variation of solutions of 

u t + A(u)u x = u xx (3.1) 

for small initial data. We always assume that the system is strictly hyperbolic, so that each matrix 
A{u) has real distinct eigenvalues Xi(u) as in (1.3), and and dual bases of right and left eigenvectors 
ri(u), k(u) normalized as in (1.4). The directional derivative of a function <j> = <j>(u) in the direction 
of the vector v is written 

( . . , 6(u + ev) — 6(u) 
v • 6(u = D6 ■ v = lim Z±± (3.2 

while 

[rj,r k ] = Tj »rjfc - r/b • r j 

denotes a Lie bracket. In order to obtain uniform bounds on Tot.Var.{u(t, •)} for all t > 0, our 
basic strategy is as follows. We choose 5o > sufficiently small and consider an initial data 
u(0, •) = u satisfying the first inequality in (2.16). By Proposition 2.3, the corresponding solution 
is well defined on the initial time interval [0, i] and its total variation remains bounded, according 
to (2.17). The main task is to establish BV estimates on the remaining interval [t, oof. For this 
purpose, we decompose the gradient u x along a suitable basis of unit vectors fi,...,r n , say 



i=l 



^Vih- (3.3) 
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Differentiating (3.1), we obtain a system of n evolution equations for these scalar components 

Vi,t + {\vi) x -v itXX =4> i i = l,...,n. (3.4) 



Since the left hand side is in conservation form, (3.4) implies 

|MV)|| L1 < |MV)|| L1 +J J \fa(t,x)\dxdt (3.5) 
for all t > i. By (3.3) it follows 

Tot.Var. {«(<,■)} = ||u*(*,-)|| L i < IN*' ^IIl* ■ ^ 

i 

In order to obtain a uniform bound on the total variation, the key step is thus to construct the 
basis of unit vectors {fi, . . . , r n } in (3.3) in a clever way, so that the functions fa on the right hand 
side of (3.4) become integrable on the half plane {t > t, x € IR}. 

As a preliminary, we observe that the choice fj = rj(u), the i-th eigenvector of the matrix 
A(u), seems quite natural. This choice was indeed adopted in [BiBl], where the authors proved 
Theorem 1 restricted to the special class of systems where all Rankine-Hugoniot curves are straight 
lines. Unfortunately, for general n x n hyperbolic systems it does not work. To understand why, 
let us write 

u x = k(u) ■ u x (3.7) 

for the i-th component of u x in this basis of eigenvectors. As shown in [BiBl], these components 
satisfy the system of evolution equations 

(<) t +(A;*4) x - {ul) xx 

= k ■ | A / '•/•'•/. >'' "i: + 2 J2(r k • rj )(ui) x u k x + r k • '■.! "' "''. "', | (3.8) 

j'/fc j,k j,k,£ 



Assume that the i-th characteristic field is genuinely nonlinear, with shock and rarefaction curves 
not coinciding, and consider a travelling wave solution u(t,x) = U(x — At), representing a viscous 
i-shock. It is then easy to see that the right hand side of (3.8) is not identically zero. Since it 
corresponds to a travelling wave, the integral 

J \<f>i(t,x)\dx + 

is constant in time. Hence fa is certainly not integrable over the half plane {t > t, x 6 M}. 

The previous example clearly points out a basic requirement for our decomposition (3.3). 
Namely, in connection with a viscous travelling wave, the source terms c/>j in (3.4) should vanish 
identically. To achieve this goal, we shall seek a decomposition of u x not along eigenvectors of the 
matrix A(u), but as a sum of gradients of viscous travelling waves. More precisely, consider 
a smooth function u : 2R i— > M n . At each point x, depending on the second order jet (u,u x ,u xx ), 
we shall uniquely determine n travelling waves U\, . . . , U n passing through u(x) (fig. 2). We then 
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write u x in the form (3.3), as the sum of the gradients of these waves. As a guideline, we shall try 
to achieve the following relations: 

Ui{x) = u(x) i = l,...,n, (3.9) 

Y J u:(x) = u x (x), J2 u iW = u **( x )- ( 3 - 10 ) 

i i 

Details of this construction will be worked out in the next two sections. 

Remark 3.1. For each i 6 {1, . . . , n}, one can find an (n+l)-parameter family of viscous travelling 
waves Ui passing through a given state u £ M n . Indeed, one can assign the speed Oi « A* and the 
first derivative U[ G M n arbitrarily, and then solve the second order O.D.E. 

U? ={A(U i )-a i )Ui 

In all, this would give us n(n + 1) scalar parameters to determine. Far too many, since (3.10) is a 
system of only 2n equations. In an attempt to fix this problem, we could restrict ourselves only to 
globally bounded travelling wave profiles. Assuming that the i-th field is genuinely nonlinear, the 
analysis in [F] or [Se] shows that there exists a 2-parameter family of viscous i-shock profiles trough 
any given point u. Indeed, one can arbitrarily assign a_, a + < and find unique asymptotic states 
u + ,u~ 6 M n which are connected by a viscous shock profile Ui(-) passing through u (fig. 3), and 
such that the following inner products take the prescribed values: 

(r*, u - u~) = a- , (r*, u + - u) = a + . 

In this case, summing over n families, we would end up with the right number of parameters to fit 
the data, namely 2n. Unfortunately, viscous shock profiles yield only negative values of gradient 
components: (r*, U[) < 0. At a point x where (r*, u x } > 0, to achieve the decomposition (3.10) 
we have to consider also viscous rarefaction profiles (which are not globally bounded). In the 
next section, a suitable family of viscous travelling waves will be selected by the center manifold 
theorem. 
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figure 3 

4 - A center manifold of viscous travelling waves 

To carry out our program, we must first select certain families of travelling waves, depending 
on the correct number of parameters to fit the data. Given a state u £ M n , a second order jet 
(u x , u xx ) determines 2n scalar parameters. In order to uniquely satisfy the equations (3.10), we 
thus need to construct n families of travelling wave profiles through u, each depending on two 
scalar parameters. This will be achieved by an application of the center manifold theorem. 

Travelling waves for the viscous hyperbolic system (3.1) correspond to (possibly unbounded) 
solutions of 

(A(U) -a)U' = U". (4.1) 
We write (4.1) as a first order system on the space M n x M n x M: 

{u = v, 
v = (A(u) - a)v , (4.2) 
a = 0. 

Let a state u* be given and fix an index i G {1, . . . , n}. Linearizing (4.2) at the equilibrium point 
P* = (n*, 0, \i(u*)) we obtain the linear system 




figure 4 
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Let {r*, . . . , r*} and {Z|, . . . , Z*} be dual bases of right and left eigenvectors of A(u*) normalized 
as in (1.4). We call (Vi, . . . , V n ) the coordinates of a vector v G M n w.r.t. this basis, so that 

k*| = i, « = E^ V/ • 

J 

The null space TV for (4.3) consists of all vectors (it, v, a) G iR™ x M n x K such that 

Vj = for all j ^ i, (4.4) 

and therefore has dimension n + 2. By the center manifold theorem [V], there exists a smooth 
manifold A4 C iR n+n+1 , tangent to M at the stationary point P* (fig. 4), which is locally invariant 
under the flow of (4.2). This manifold has dimension n + 2 and can be locally defined by the n — 1 
equations 

Vj = <pj(u,Vi,<T) j^i. (4.5) 

We can assume that the n — \ smooth scalar functions ipj are defined on the domain 

V = {\u-u*\ < e, \Vi\ < e, |cr - A;(it*)| < e}. 

Moreover, the tangency condition implies 

^(u,V l ,a)=0(l)-(\u-u*\ 2 + \V l \ 2 + \a-\(u*)\ 2 y (4.6) 

We now take a closer look at the flow on this center manifold. By construction, every trajectory 

t^P(t)=(u(t),v(t),a(t)) 

of (4.2), which remains within a small neighborhood of the point P* = (it*, 0, Ai(u*)) for all t <E M, 
must lie entirely on the manifold M. In particular, M contains all viscous i-shock profiles joining 
a pair of states u~,u + sufficiently close to it*. Moreover, all equilibrium points (it,0,cr) with 
| it — it* | < e and \a — Xi(u*) | < e must lie on M.. Hence 

(fj(u,0,a) = for all j^i. (4.7) 

By (4.7) and the smoothness of the functions ipj, we can "factor out" the component Vi and write 

ipj(u,Vi,a) = ^j(u,Vi,a) • V, 

for suitable smooth funtions ipj. From (4.6) it follows 

ipj as (it, Vi, a) -» («* , 0, A; (it*)) ■ (4.8) 

On the manifold .M we thus have 

v = ^V k r* k = V- lr*+J2*Pj(u,Vi,a)r* J = V, rf(it, ^, a). (4.9) 
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By (4.8), the function r" denned by the last equality in (4.9) satisfies 

rl(u,Vi,a) -> r* as (u, V i: a) -> («*, 0, Aj(u*)). 



(4.10) 



Remark 4.1. Trajectories on the center manifold correspond to the profiles of viscous travelling 
i-waves. We thus expect that the derivative u = v should be a vector "almost parallel" to the 
eigenvector r* = ri{u*). This is indeed confirmed by (4.10). 

We can now define the new variable 

v i = v i {u,V U G) = V i -\r\(u,V u a)\. (4.11) 

As (u,Vi,a) range in a small neighborhood of («*, 0, Aj («*)), by (4.10) the vector rf remains 
close to the eigenvector r*. In particular, its norm remains uniformly positive. Therefore, the 
transformation Vi < — ► Vi is invertible and smooth. We can thus reparametrize the center manifold 
M in terms of the variables (u, Vi,a) G M n x ]R x M. Moreover, we define the unit vector 

fi(u,Vi,a) = . (4.12) 

vi I 

Observe that fj is also a smooth function of its arguments. With the above definitions, instead of 
(4.5) we can write the manifold M in terms of the equation 

v = Vifi. (4.13) 

The above construction of a center manifold can be repeated for every i = 1, . . . , n. We thus 
obtain n center manifolds Mi C M 2n+1 and vector functions fj = fi(u,Vi,(Ti) such that 

\n\ = 1, (4.14) 

Mi = {(u,v,ai) ; u = ■y i fi(«,v i ,(7 i )} , (4.15) 
as (u,Vi,<Ji) G iR n x 1R x 1R ranges in a neighborhood of (u*, 0, Aj («*)). 



We derive here some useful identities, for later use. The partial derivatives of fj = fj(«, Uj, <Tj) 
w.r.t. its arguments will be written as 

„ . d _ _ . d _ _ . 3 _ 

l.U r\ 'It 1 l.V r\ I ) ' l,<7 1 ' 

OU OVi OOi 

Clearly, f ijU is an n x n matrix, while f ijV , fi t<7 are n- vectors. Higher order derivatives are denoted 
as r ijUcr , f iiCrcr . . . We claim that 

fj (u, 0, (jj) = Vi(u) for all u,<jj. (4-16) 

Indeed, consider again the equations for a viscous travelling i-wave: 

u xx = (A(u) - (Ji)u x . (4.17) 
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For a solution contained in the center manifold, taking the derivative w.r.t. x of 

u x = v = v i r i (u,v i ,<Ji) (4.18) 

and using (4.17) we obtain 

Vi, x fi + Vif itX = (A(u) - (TijVifi. (4.19) 

Since |fj| = 1, the vector is perpendicular to all of its derivatives. Taking the inner product of 
(4.19) with fi we thus obtain 

Vi, x = (K ~ <Ti)Vi , (4.20) 
where we defined the speed A, = Xi(u, Vi, <Jj) as the inner product 

A, = (ft , A(u)h) ■ (4.21) 
Using (4.20) in (4.19) and dividing by Vi we finally obtain 

(Aj - ai)vifi + Vi(fi tU fiVi + fi tV (Xi - <Ti)vi) = (A(u) - a^v^i , (4.22) 

Vi(fi, u fi +r itV (\i - Oi)) = (A(u) - Xi)fi . (4.23) 

By (4.23), as Vi — > 0, the unit vector fi(u,Vi,ai) approaches an eigenvector of the matrix A(u), 
while A, approaches the corresponding eigenvalue. By continuity, this establishes (4.16). 
In turn, by the smoothness of the vector field fi we also have 

fi(u,v i ,a i )-r i (u) = 0(l)-Vi, f i ^ = 0{l)-v il 

fi,ua = 0{l) ■ Vi, f ita(T = 0(1) • Vi. 

Observing that the vectors and f ij(T are both perpendicular to fj, from (4.24) we deduce 

\Xi(u,Vi,ai)-Xi(u)\=0(l)-Vi, \ v = 0(1) ■ v u \ a = 0(1) ■ vf. (4.25) 
A further identity will be of use. Differentiating (4.19) one finds 

Vi,xxh + 2v i)X f iyX + Vif iyXX = (A(u)vifi) x - OiV^fi - aiVifi )X . (4.26) 
From the identities 

taking the inner product of (4.19) with fj x we obtain 

(fi, f i>xx )vi = -(f i>x , A(u)fi)vi . (4.27) 
Taking now the inner product of (4.26) with f, we find 

Vi,xx + (?ii Ti,xx) v i = (^ii {A(u)fiVi) x 'j — <TiVi )X . 

Since v ijt + <JiV itX = 0, using the identity (4.27) we conclude 

Vi, t + (XiVi) x - v iyXX = 0, (4.28) 

where A, is the speed at (4.21). 

Remark 4.2. It is important to appreciate the difference between the identities 

(A(u)-Xi)ri = 0, (A(u)-Xi)f l = Vi(fi tU fi + fi tV (Xi-a t )), (4.29) 

satisfied respectively by an eigenvector and by a unit vector fj parallel to the gradient of a 
travelling wave. Decomposing u x along the eigenvectors one obtains the evolution equations 
(3.8), with non-integrable source terms on the right hand side. When a similar computation is 
performed in connection with the vectors fj, thanks to the presence of the additional terms on 
the right hand side in (4.29) a crucial cancellation is achieved. In this case, we will show that the 
source terms (f>i in (3.4) are integrable over the half plane x 6 M, t > i. 
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5 - Gradient decomposition 

Let u : M i— > M n be a smooth function with small total variation. At each point x, we 
seek a decomposition of the gradient u x in the form (3.3), where = ?i(u,Vi, Uj) are the vectors 
defining the center manifold in (4.15). To uniquely determine the fj, we should first define the 
wave strengths Vi and speeds cr, in terms of u, u x , u xx . 

Consider first the special case where u is precisely the profile of a viscous travelling wave of 
the j-th family (contained in the center manifold M.j). In this case, our decomposition should 
clearly contain one single component: 

u x = Vjrj(u,Vj,<Tj) . (5.1) 

It is easy to guess what Vj,aj in (5.1) should be. Indeed, since by construction \fj\ = 1, the 
quantity 

Vj = ±\u x \ 

is the signed strength of the wave. Notice also that for a travelling wave the vectors u x and u t are 
always parallel, since u t = —(JjU x where <jj is the speed of the wave. We can thus write 

u t = u xx - A(u)u x = ujjfj(u, Vj,aj) (5.2) 

for some scalar ujj. The speed of the wave is now obtained as Uj = —Uj/vj. 

Motivated by the previous analysis, as a first attempt we define 

ut = u xx - A(u)u x (5.3) 
and try to find scalar quantities , u>i such that 

u x = ^2viri(u,Vi,ai), 

* ^ = --- (5-4) 

u t = 2_^uJir i (u,Vi,a i ), v * 

i 

The trouble with (5.4) is that the vectors fj are defined only for speeds <jj close to the i-th 
characteristic speed A* = \i(u*). However, when u x « one has Vi ~ and the ratio Ui/vi may 
become arbitrarily large. 
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To overcome this problem, we introduce a cutoff function (fig. 5). Fix <5i > sufficiently small. 
Define a smooth odd function : M>-^ [— 26i, 28\] such that 



We now rewrite (5.4) in terms of the new variable Wi, related to by u>i = Wi — A*«j. We require 
that <jj coincides with —ooi/vi only when this ratio is sufficiently close to A* = \i(u*). Our basic 
equations thus take the form 

i (5-6) 
ut = ^2(wi - X*Vi)fi(u,Vi,ai), 



where 



u t = u xx - A(u)u x , Ui = X* - 6 ( — ) . (5.7) 



Vi 



Notice that Oi is not well defined when Vi = Wi = 0. However, recalling (4.16), in this case we have 
fj = ri(u), regardless of Uj. Hence the two equations in (5.6) are still meaningful. 

Remark 5.1. The decomposition (5.6) corresponds to viscous travelling waves Ui such that 

Ui(x) = u(x), U' i {x) = v i r i , U" = (A(u) - ai)U- . 

From the first equation in (5.6) it follows 

u x {x) = Y J U[{x) . 

i 

If Oi = A* — Wi/vi for alH = 1, . . . , n, i.e. if none of the cutoff functions is active, then 

u xx {x) =u t + A(u)u x 

= - X*Vi)fi + A(u) ^2 Vih 

i i 

i 

i 

In this case, both of the equalities in (3.10) hold. Notice however that the second equality in (3.10) 
may fail \i\wi/vi\ > 8\ for some i. 

Lemma 5.2. For \u — u*\, \u x \ and \u xx \ sufficiently small, the system o/2n equations (5.6) has a 
unique solution (v, w) = (v\, . . . , v n , Wi, . . . , w n ). The map (u, u x ,u xx ) t— > (y, w) is smooth outside 
the n manifolds Mi = {vi = Wi = 0}; moreover it is C 1 ' 1 , i.e. continuously differentiate with 
Lipschitz continuous derivatives on a whole neighborhood of the point (n*,0,0). 
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Proof. Given (v, w) in a neighborhood of (0, 0) 6 M 2n , the vectors u x ,u t are uniquely determined. 
Hence the solution of (5.6)-(5.7) is certainly unique. To prove its existence, consider the mapping 



A : XT x M n x ]R n ^ M 2n defined by 



A(u,v,w) =y^Ai(u,Vi,Wi), (5.8) 



i=l 



\(wi - \*Vi)ri{u, Vi, A* - 0(wi/vi)) J 

This map is well defined and continuous also when vi = 0, because in this case (5.16) implies 
fi = Ti{u). Computing the Jacobian matrix of partial derivatives w.r.t. (vi,Wi) we find 

dAi ( fi 



d(vi,Wi) f iy (51Q) 

, / v i f i ^ v + (w i /v i )e' i f i ^ -Q'ih,o 

\Wifi,v - KVifi,v - \*(Wi/Vi)e^f iya + {Wi/VifO'^i^ A*6^; >CT - (Wi/Vi)6[fi 

Here and throughout the following, by 6i, 6^ we denote the function 9 and its derivative, evaluated 
at the point s = Wi/vt. By (5.10) we can write 

OA 

B (u,v,w) + B 1 (u,v,w), (5-11) 



d(v, w 



Because of (4.24), the matrix functions B ,Bi are well defined and continuous also when Vi = 0. 
Moreover, for (v, w) small, Bq has a uniformly bounded inverse and B\ — > as (v, w) — ► 0. Since 
A(u, 0, 0) = 06 M 2n , we conclude that the map (v,w) i— > A(u; v,w) is C 1 and invertible in a 
neighborhood of the origin. Therefore, given (u,u x ,u xx ), there exist unique values of (v,w) such 
that 

A(u,v,w) = (u x , u xx - A(u)u x ). (5.12) 

The inverse of the map A w.r.t. the variables v,w will be denoted by A _1 (u; p, q). In other 
words, 

A _1 (u; p,q) = (v,w) iff A(u; v, w) = (p, q) . 
Since fj(tt, 0, <Tj) = we have 

A(u,0,w) = ( , ^t^r^u)). 

i 

Therefore, 

A _1 (-u, 0, (/) = (0, w) where W{ = k(u) • q . 

In particular, ^(-u, 0, 0) = (0,0) G iR 2n . Concerning first derivatives (which we regard here as 
linear operators), we have 

d ^d{v w)^ ' ^' ^ = Bo<yU] °' W>> ' ^' ^ = f ^ ' (™» _ ' ( 5 - 13 ) 
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dA (u, 0,q) _ = where 0; = k{u) - pi, u>i = Z;(u) • & . (5.14) 

We shall not compute the second derivatives explicitly. However, one easily checks that 
d 2 A d 2 A d 2 A 



dvidvj dvidwj dwidwj 



= if i + j. (5.15) 



Moreover, recalling (4.24) and (5.5), we have the estimate 

&A d 2 A cPA = (516) 
dv 2 ' dvidwi ' dwf S\ ' 

Since the cutoff function 9 vanishes for |s| > 3<5i, it is clear that each Aj is smooth outside the 
manifold Mi = {(v,to) ; = Wi = 0}, having codimension 2. Since all second derivatives are 
uniformly bounded outside the n manifolds Mi, we conclude that A is continuously differentiable 
with Lipschitz continuous first derivatives on a whole neighborhood of the point (u*,0, 0). Hence 
the same holds for A -1 . □ 

Remark 5.3. By possibly performing a linear transformation of variables, we can assume that 
the matrix A(u*) is diagonal, hence its eigenvectors r*, . . . r* form an orthonormal basis: 

(r*,r*) = 6 ij . (5.17) 

Observing that 

>*\=0(l)-(\u-u*\ + \vi\), (5.18) 



\ri{u,Vi,ai 



from (4.16) and the above assumption we deduce 

(fi(u,Vi,ai), fji^v^aj)) = Sij + 0(1) • (|u-u*| + H + |^|) 

5.19 

= S i:i + O(l) ■ So , 

{f i ,A{u)f j ) = O{l)-S j + i. (5.20) 

Another useful consequence of (5.17)-(5.18) is the following. Choosing 5 > small enough, the 
decomposition (5.6) will satisfy 

\u x \ < \vj\ < 2y/n\u x \ . (5.21) 



We conclude this section by deriving estimates corresponding to (2.13)-(2.15), valid for the 
components Vi,Wi. In the following, given a solution u = u(t, x) of (3.1) with small total variation, 
we consider the decomposition (5.6) of u x in terms of gradients of travelling waves. It is understood 
that the vectors fj are constructed as in Section 4, taking P* = (u*, 0, Xi(u*)) as basic points in 
the construction of the center manifolds Mi. Here u* = u(t, — oo) is the constant state in (2.1). 

Lemma 5.4. In the same setting as Proposition 2.1, assume that the bounds (2.6) hold on a larger 
interval [0,T]. Then for all t 6 [t, T], the decomposition (5.6) is well defined. The components 
Vi,Wi satisfy the estimates 

|K(t)|| Ll , |k(*)|| L i = O(l)-S , (5.22) 
||«i(*)|| LO o, hi(t)\\^,\\vi, x (t)k^ hUVWv = O(l)-S 2 , (5.23) 

II^WLo., R*(*)|U = ^(i)-^- (5-24) 
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Proof. By Lemma 5.2, in a neighborhood of the origin the map (v, w) *— > A(u, v, w) in (5.8) is 
well defined, locally invertible, and continuously different iable with Lipschitz continuous deriva- 
tives. Hence, for So > suitably small, the L°° bounds in (2.13) and (2.14) guarantee that the 
decomposition (5.6) is well defined. From the identity (5.12) it now follows 

Vi,Wi = 0(1) • (|« x | + \u xx \). 

By (2.6) and (2.13)-(2.14) this yields the L 1 bounds in (5.22) and the L°° bounds in (5.23). 
Differentiating (5.12) w.r.t. x we obtain 

^u Ux + d(v~w) ( Vx,Wx ^ = { Uxxi u *xx - A(u)u xx - (u x • A(u))u x ^j . (5.25) 
Using the estimate 

^ = o(i).(H + H), 

since the derivative dA/d(v,w) has bounded inverse, from (5.25) we deduce 

(v x ,w x ) = 0(1) ■ (\u xx \ + \u xxx \ + \u x \ 2 + \u x \(\v\ + \w\)Y 
This yields the remaining estimates in (5.23) and (5.24). □ 



6 - Bounds on the source terms 

We now consider a smooth solution u = u(t,x) of (3.1) and let Vi, to*, be the corresponding 
components in the decomposition (5.6), which are well defined in view of Lemma 5.2. The equations 
governing the evolution of these 2n components can be written in the form 



j V it + {\Vi) x - v iiXX = (pi, 

< - (6.1) 

[ W itt + {\iWi) x - Wi, xx = tpi . 

As in (4.21), we define here the speed Aj = (fj , A(u)fj) . The source terms <pi,ipi can be computed 
by differentiating (3.1) and using the implicit relations (5.6). However, it is not necessary to carry 
out in detail all these computations. Indeed, we are interested not in the exact form of these terms, 
but only in an upper bound for the norms H^Hl 1 and HV^IIl 1 • 

Before giving these estimates, we provide an intuitive explanation of how the source terms 
arise. Consider first the special case where u is precisely one of the travelling wave profiles on the 
center manifold (fig. 6a), say u(t,x) = Uj(x — <jjt). We then have 



and therefore 



u t = (wj - X*jVj)rj , Vi = Wi = for i ^ j , 

f v i>t + (\Vi) x - v i>xx = , 
1 w i)t + (XiWi) x - w i)XX = . 
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U=U; 



X x' 



figure 6a 




X 

figure 6b 



Indeed, this is obvious when i ^ j. The identity <f>j = follows from (4.28), while the relation 
Wj = (A* — (Jj)vj implies tpj = 0. 

Next, consider the case of a general solution u = u(t,x). The sources on the right hand sides 
of (6.1) arise for three different reasons (fig. 6b). 

1. The ratio \wj/vj\ is large and hence the cutoff function 6 in (5.7) is active. Typically, this will 
happen near a point xo where u x = but u t = u xx ^ 0. In this case the identity (4.28) fails 
because of a "wrong" choice of the speed: <jj ^ A* — (wj/vj). 

2. Waves of two different families j ^ k are present at a given point x. These will produce 
quadratic source terms, due to transversal interactions. 

3. Since the decomposition (3.10) is defined pointwise, it may well happen that the travelling 
j-wave profile Uj at a point x is not the same as the profile Uj at a nearby point x' . Indeed, 
these two travelling waves may have slightly different speeds. It is the rate of change in this speed, 
i.e. (Jj tX , that determines the infinitesimal interaction between nearby waves of the same family. 
A detailed analysis will show that the corresponding source terms can only be linear or quadratic 
w.r.t. aj tX , with the square of the strength of the wave always appearing as a factor. These terms 
can thus be estimated as 0(1) • v^aj^ + 0(1) • Vj<r? x . 



24 



Lemma 6.1. The source terms in (6.1) satisfy the bounds 

<f>i, ipi = 0(1) ■ ^2 (\vj, x \ + \wj, x \) ■ \wj - 6jVj\ (wrong speed) 



+ o(i)lK, 



X ^ j , X | 







(-) 







■^{iWi/^IO*!} 



(change in speed, linear) 
(change in speed, quadratic) 



+ 0(1) • ^2 {\ V j V k\ + \Vj, x V k \ + \VjW k \ + \v jyX W k \ + \UjW kyX \ + \WjW k \) 

0+k 

(interaction of waves of different families) 



From a direct inspection of the equations (6.1), it will be clear that the source terms depend 
only on the third order jet (n, u x , u xx , u xxx ). Since all functions fa, ipi vanish in the case of a 
travelling wave, for a general solution u their size can be estimated in terms of the distance between 
the third order jet of u and the (nearest) jet of some travelling wave. This is indeed the strategy 
adopted in the following proof. An alternative proof, based on more direct calculations, will be 
given in Appendix A. 

Proof of Lemma 6.1. The conclusion will be reached in several steps. 

1. The vector (u x ,u t ) = A(u,v,w) satisfies the evolution equation 
A{u) 



Ux 



+ 









A(u) 



Ux 
U t 







(u x • A(u))u t - (u t • A{u))u a 
Observe that, in the conservative case A(u) = Df(u), the right hand side vanishes because 

(u x • A(u))u t = (ut • A(u))u x = D 2 f(u) (u x ® u t ) . 
In the general case, recalling (5.6) we deduce 



(6.3) 



Ut.7 



A(u))u t - (u t • A{u))u x = 0(1) • ^ (\ V 3 V k\ + \VjW k \). 



(6.4) 



2. For notational convenience, we introduce the variable z = (v,w) and write A for the 2n x 2n 
diagonal matrix with entries Aj defined at (4.21): 



From (6.3) it now follows 



OA 



dA 



TT-Ut + TT- . 

du dz \ w 



+ 



A = 



A(u) 




d 2 A 



(u x <S> u x ) - 



diag(Ai 






A(u) 
d 2 A 





diag(Ai 



A 



dA f v 
dz \w 



u x 





(u x • A(u))u t - (u t • A(u))u 



u x 
w :1 



- 2 



dA 

du 

d 2 A 

du dz 
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Therefore, 



<9A 



w 



+ A 



w 



dz \ \w 



A(u) 
A(u 

d 2 A 



A) +-^A(u)u a 



d 2 A 

, + 2 Q „ u x ® 



«3 



= £ 

(6.5) 

Since the differential dA/dz has uniformly bounded inverse, the right hand sides in (6.1) clearly 
satisfy the bounds 



<f>i = 0(l)-E, 



i>i = 0(1) ■ E, 



i = 1, . . . , n . 



(6.6) 



3. To estimate the quantity E in (6.5), it is convenient to introduce the function 



Ai(u,Vi,Wi,CTi) 

so that A = Yl Aj and E = X) -^i > where 



Ei 



dAi ( L u 



dz; \ XiW 



+ 



' A(u) 
A(u) 





A, 



Vifi(u,Vi,(Ti) 



<9Aj ~ <9Aj ^-^ 



7Tj((u x • - (U t • A(u))u x ) J dvB 



d<Ji 
d 2 A t 



u x + 



d 2 A % 



d 2 A, , 9 2 A, 



<9cT ? ; <9z ? ; V <Ji, x W. 



<9u dz; V m 



u x + 2 



Wi. 



d 2 Aj 
du daj 



. X "'CI. 



+ 



dAi ( d 2 Oi 2 n d 2 Oi d 2 (T. 



vL + 2 



.,2 "1,2! 



dvidwi 



Vi,xWi^ x 



i w 2 

dw 2 l ' x 



(6.7) 



(6.8) 



Notice that in (6.5) we regarded A as a function of the three independent variables (u, v, w), while 
in (6.8) we think A as a function of the four independent variables (u, v, w, a). Regarding the <7j as 
independent variables, one has the advantage that the maps Aj = Aj(«, Vi,Wi,<Ji) are now smooth, 
while Aj = Ai(u, v, w) in (5.8) was only C 1 ' 1 , because of the singularities of the map (u, Vi,Wi) <— > <7j 
in (5.7). The last term in (6.8) is due to the nonlinear dependence of <jj w.r.t. Vi,Wi. By 7Tj(v) 
we denoted the i-th component of a vector v w.r.t. the basis {rf, . . . , r*}. Also notice that in the 
previous computation we used the identity 

dAi doi ~ dAi doi ~ 

-Ai.xVi + ~ — ~ M. x Wi 



d<Ji dv, 
= A 



doi dw, 
Wi/vi ■ Q'fii^ 
'-' 1 Wi/vi(wi/vi - A*)6^f 



i' i,<? 



4. By Lemma 5.2, the inverse map A 1 sets a one to one correspondence 

(u,u x ,u xx ) i ^ 
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between two neighborhoods of the point (u*,0, 0) G M 3n . This map is C 1 with Lipschitz continuous 
derivative. It can be prolonged to a map 

{u,u x ,u xx ,u xxx ) i ^ (u,v,w, a,v x ,w x ,a x ) (6.9) 

which is one-to-one, but of course not onto. Indeed, (5.6) and the identity ut + A{u)u x = u xx 
together imply 

i ij i i 

A vector («,«, w,a,v x , w x ,cr x ) £ M 7n corresponds to some third order jet (u,u x ,u xx ,u xxx ) pro- 
vided that it satisfies the vector equation (6.10), together with 

x* a ( W A WiVi, x - w i:X Vi , ( wA . 
a i = \ i -6\ — \, a i;X = '■ — g — •■ 9 I — i = l,...,n. (6.11) 



5. By the analysis at (6.2), Ei(u, v ^ , , , v£, w% , = whenever the argument corresponds 
to the third order jet of a viscous travelling 2-wave. This is the case if 



*?* = °> 







.,0 



for all j ^ i, 







(6.12) 



<SS U a? = \*-^-, <t9 x = q. (6.13) 



In order to estimate Ei(u,v,w,a,v x ,w x ,a x ) we proceed as follows. We introduce a new vector 
(u, , , a^, v£, w£, a%) corresponding to the jet of a travelling i-wave, by setting 



V? = Vi 



w- 



Wi 



\ * 

Vi, a? = a i = \ 



(6.14) 



The quantities vf ,wf x , af x are then defined according to (6.13), while the components j ^ i 
are as in (6.12). The above construction implies Ef = Ei(u,v^ ,v£ ,w% = 0. Hence 



E i = E i -Ef. 

6. Taking the inner product of (6.10) with fj, recalling that r, has unit norm and is thus orthogonal 
to its derivatives, we obtain 

Wi + (A; - \*)Vi = V ijX + 6, 

where 



J#« j'/i k 

= O(l)-5 J2\v j \- 



(6.15) 
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The above estimate on is obtained using (5.20) together with the L°° bounds in (5.23)-(5.24) 
and the bound on f J)Cr in (4.24). We can now write 



Vi, x = Wi + (Aj - \*)vi + 0(1) • 5 ^2 \vj\ 
= (Aj - <Ti)vi + {wi - O&i) + 0(1) • 5 



We recall that Bi = 9(wi/vi). The first equality in (6.16) yields the implications 



\wA < 3Si\vi 



\wA > SAvi 



0(1)- Vi + O(l)-6 I 



v i = O(l)-v i;X + O(l)-S J2 



Moreover, from the second equality in (6.16) we deduce 



(A; - Ui)wi = (Aj - <Ji) [v iiX - (Aj - X*)vi] + 0(1) • S ^ 



(Aj - <7;)?;; >a; - (Aj - A*)(w ijX - (to* - 0^)) + 0(1) • 5o^N 



(6.16) 



(6.17) 



(6.18) 



— Vi,x ~ ( — - 0i ) ^ + (A; - A*)(u>; - 0^) + 0(1) • <5 ^ N , 



and hence, by (6.18), 



Wi,* - (Aj - crO^i 



+ O(l)-\w i -B i v i \+O(l)-S J2\v j \. 



From the definitions (6.13)-(6.14), using the above estimates we obtain 



Wi -w?\ = w 



Owl +0(1)- S J2\ Vj \, 
0(1)- \ Wi -B iVi \ + 0(1)- S 



J i,x\ 



0(1)- K-^| + 0(1)-<5 O ^|^| 

3^i 



(6.19) 



28 



7. We now compute 



E t = E t - E? 

_ <9Aj / X^i \ _(\A(u) 

" dzi \hwi) x \[ A(u) 
dAi ~ dAi tt-^ _ 



A, 



dA? ( \?v? 



+ 



+ 



+ 





■Ki((u x • A(u))u t - (u t • A(u))u x ) 
d 2 A r ( Vi -\ \ d 2 A? 



d 2 A i 
+ I -tt^tU x &> u x 



A{u) 
A(u) 



A? 



^f 1 

d 2 A 



i,x 



ll'i 



dz, 



[2] 



Vi 



u x -2 



+ 2 



3 2 A,; 



dAi fd 2 a. 



V 

w 

d 2 Af l v 
dudzi \ w 
d 2 a l 





i,x 

l.X 





i,x 



i , x 



Viff + 2 



d 2 A l 



du d(Ji 



dvidwi 



8 2 a 



% 2 



-Wi 



(6.20) 

Observe that the quantities u, V{, &i, fj, \ remain the same in the computations of Ei and Ef . 
Moreover, all the terms involving derivatives w.r.t. vanish when we compute Ef . 

In the remaining steps, we will examine the various terms on the right hand side of (6.20) and 
show that they can all be bounded according to the lemma. As a preliminary, we observe that by 
(6.7) and (4.24) the derivatives of the smooth function Aj = Ai(u, 2*, <7j) satisfy 



dAi d 2 Ai d 2 Ai 



du ' du& ' dzidoi 



(H + H), 



(6.21) 



dAi d 2 Ai d 2 Ai ,. „ 



(6.22) 



8. We start by collecting some transversal terms. Using (6.4) and (6.21)-(6.22) we obtain 



dAi y—^ 

■ ^A(u)v jrj + 



■Ki((u x • A(u))u t - (u t • A(u))u x ) 



d 2 A t _ _ 2 _ \ d 2 Ai ^ 

+ ~ ~ I ^VjVkrj <S>r k -v i r i <g>r i \+ ^j—- ^ 



Wjfj (6.23) 



3\ k J J#* 



Here and in the following, by "transversal terms" we mean terms whose size is bounded by products 
of distinct components j ^ k, as in (6.23). 



9. We now look at terms involving derivatives w.r.t. <7j. One should here keep in mind that, if 
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&i,x 7^ 0, then both sides of the implication (6.16) hold true. Using (6.21) we obtain 

d 2 A, 

"i,x w i,x 

' 2 'n,x t ■ yujiu iyX - 'w itX ui 

2 



u z,x v i,x 



0(1) • Vi(\v iiX \ + \wi, x \)a itX 
0(1) ■ v 2 a i>x + 0(1) • (wiVi yX — Wi yX Vi)ai )X + transversal terms 

0(1) • \WiV itX - Wi !X Vi\ + 







Vi ■ 


(-) 







■Y r , , , + transversal terms . 



An application of (6.22) yields 



(6.24) 



dAi ~ d 2 Ai 2 d 2 Ai , n . , 9 ,. . 9 

\<?i,x + -K-T^i, x + 2 o o &i,xU x = 0(1) • vt{\(Ji, x \ + \<J i)X \) + transversal terms 



<9cr. 



da} ~ 1 > X ' 'dudai 

0(1) • \ViW itX - WiV itX \ + 0(1) 









(-) 







•Y r , , , r , + transversal terms . 



Next, we observe that the quantity 



(6.25) 



d 2 



9*Vi %x dv.dwi 



Vi. x Wi, x 4 



a 



w r 



i ...2 



vanishes in the special case where Wi iX = (wi/vi)vi tX . In general, recalling (6.16) a direct compu- 
tation yields 



dAi (d 2 a. 



'i 2 



d 2 



Gi 



d<Ji \ d 2 Vi l ' x dvidwi 



Vi, x Wi x -\- 



& 

w 

d 2 Wi l ' x 



1 „..2 



dAi 



i 4 I w i x 

vf ' 



Vi 



(6.26) 



0(1) • |vi«;i, x - WiV iiX \ + 0(1) 













\ V i)x 



• y , , r , + transversal terms . 



10. We now complete the analysis of the remaining terms. As a preliminary, we observe that 
the only difference between A^ and A, is due to the fact that one may have wf ^ lOj. The first 
estimate in (6.19) thus implies 



\Af-Ai\, \DAf-DAi\, \D 2 A? -D 2 A l \ = 0(1) ■ \ Wl - 9 lVl \ + 0(1) ■ 5 ^\ Vj \ . 



(6.27) 



By (6.17) and (6.27), if we compute A, or its partial derivatives at the point (u,Vi,Wi,<Ti) instead 
of (u,vf ,wf ,af) = (u,Vi,wf ,crj), the difference in each of the corresponding terms in (6.20) will 
have magnitude 

0(1) • \wi — 9iVi\ ■ + |«>i,x|) + transversal terms . 

For example, 



d 2 Ai d 2 A <> \ 

dvM ~ ~dvW I ® = ®^ ' ][Wi ~ Vi + transversal terms 

= 0(1) • \wi — 9iVi\ v 2 x + transversal terms . 



(6.28) 
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Indeed, if Wi ^ 9iVi, then both sides of the implication (6.17) hold true. 
Observing that d 2 Ai/dw 2 = and recalling (6.17), we have 



d 2 A t 



dz. 



[2] 



Vi, 



If 



t.x 



W; 



l,X 



■J? 
i,x 



= 0(1) ■ (v 2 x - (vtf) + 0(1) ■ (v hxWl , x - v&v,*,) 

= 0(1) • Vi tX \wi — 9iVi\ + 0(1) • Wi tX \wi — 9iVi\ + 0(l)\wiVi tX — ViWi tX \ + transversal terms . 

(6.29) 

In a similar way, using (6.17) and (6.19) one derives the estimate 



d 2 Ai 



<> ) <S> Vih = 0(1) • Vi, x (wi - QiVi) + 0(1) • w iiX \wi 



(6.30) 



+ 0(1) • \wiVi tX — ViWi tX \ + transversal terms . 



Using the identity (4.23), we now compute 



OA, ( X 



dzi \ XiW 



+ 



h o 



A{u) 
A(u) 



A{u) 
A(u 

OA. 



A, 



I-A(u) Vifi,v 

I-A(u)\ [(wi - X*Vi)fi,v 

fi,uh + {\ - <?i)ri,v 
■K(ri,uh + - af)f itV ) f i)U fj + (Xi - <Ji)r itV 



da. 



■en,. 



E 



DA(u)fj 







DA(u)fj 



A(u) 
A(u) 



XiVi 
XiWi 

Wi x 



dAi _ 



With similar arguments as above, we obtain 



(XiVi) x Vi(wi - X*Vi)r ijV 



C*?v?) x v?(w?-X*v?)rl 



0(1) • \viW i:X - v i>x Wi\ ■ X{\ w . /v .\ <3Sl } + C> ( 1 ) ' v i,x\ w i ~ @i v i\ + transversal terms, 



Xi.x.Wi 



X^ v? 

i,x i 
Xfwf 



0(1) • \viWi jX — Vi iX Wi\ + 0(1) • Vi jX \wi — 9iVi\ + transversal terms . 



ViWi.,. 



% l,X 



- 0(1) • \viU>i yX — Vi yX Wi\ + 0(1) • Vi, x \wi — 0iV%\ + transversal terms . 
The above estimates together imply 



dAj ( X^i 
dzi \XiWi 



A{u) 
A(u) 



A, 



5Af / Xf v? 



A{u) 
A(u) 



A? 



= 0(1) • \wi — 6iVi\ (\vi, x \ + l^i.xl) + 0(1) ' \wiVi, x — ViWi tX \ + transversal terms . 
This completes the proof of Lemma 6.1. 



(6.31) 



□ 
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7 - Transversal wave interactions 



The goal of this section is to establish an a priori bound on the total amount of interactions 
between waves of different families. More precisely, let u = u(t, x) be a solution of the parabolic 
system (3.1) and assume that 



K(i)|| L1 < S 



t€[0,T\ 



(7.1) 



In this case, for t > i, by Corollary 2.2 all higher derivatives will be suitably small and we can 
thus define the components Vi,W{ according to (5.6)-(5.7). These will satisfy the linear evolution 
equation (6.1), with source terms (f>i,tpi described in Lemma 6.1. Assuming that 

r-T r 

\(f>i(t,x)\ + \il)i(t,x)\dxdt<5 , i = l,...,n, (7.2) 

and relying on the bounds (5.22)-(5.24), we shall prove the estimate 

/ 1^2 (\ v J v k\ + \ v j,xVk\ + \vjW k \ + \vj, x w k \ + \vjW k , x \ + \wjW k \) dxdt = 0(1) • 5% . (7.3) 



As a preliminary, we establish a more general estimate on solutions of two independent linear 
parabolic equations, with strictly different drifts (fig. 7). 



Lemma 7.1. Let z,z^ be solutions of the two independent scalar equations 

Zt + (\(t,x)z) x - z xx = <p(t,x) , 



+ (\t(t,x)z*) x -zl x = <(S(t,x), 



defined for t £ [0, T]. Assume that 



infA^^x) - supA(t,x) > c > 0. 

*> x t, x 



Then 



\z{t,x)\ |z"(t,x)| dxdt 



(7.4) 



(7.5) 



\z(0,x)\dx + J J \<p(t,x)\ dxdt J (J \z i (0,x)\dx + J J \^{t,x)\dxdt 



(7.6) 




figure 7 
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Proof. We consider first the homogeneous case, where ip = = 0. Define the interaction potential 

Q(z,J) = j J ' K(x-y)\z(x)\\z*{y)\dxdy, (7.7) 

with 

f 1/c if s > 0, oX 

K ( 8 ) = ; f - n ( 7 - 8 ) 



l/c-e cs / 2 if s<0 

Computing the distributional derivatives of the kernel K we find that cK' — IK" is precisely the 
Dirac distribution, i.e. a unit mass at the origin. A direct computations now yields 



d_ 

It 



Q Wf ),,« W )<-/|*x)||,«(M)kx. 



Therefore 



1/ 



z(t,x)\ ^{t^x^dxdt < Q(z(0), z tt (0)) < - ||^(0)|| Ll ||^(°)|| L i • ( 7 - 9 ) 



proving the lemma in the homogeneous case. 

To handle the general case, call T, r" the Green functions for the corresponding linear ho- 
mogenous systems. The general solution of (7.4) can thus be written in the form 



:(/../•) - / T(t,x,0,y)z(0,y)dy + J J T(t,x,s,y)<p(s,y)dyds, 



(7.10) 



z*(t, x) x,s,y)ip$(s,y)dyds. 
From (7.9) it follows 

j jT(t,x,s,y)-T i (t,x,s',y')dxdt<- (7.11) 

for every couple of initial points (s,y) and (s',y'). The estimate (7.6) now follows from (7.11) and 
the representation formula (7.10). □ 

Remark 7.2. Exactly the same estimate (7.6) would be true also for a system without viscosity. 
In particular, if 

z t + (A(t, x)z) x = 0, z\ + (A» (t, x)z*) x = 0, 

and if the speeds satisfy the gap condition (7.5), then 



d_ 
~dt 



x<y 



In the case where viscosity is present, our definition (7.7)-(7.8) thus provides a natural counterpart 
to the Glimm interaction potential between waves of different families, introduced in [G] for strictly 
hyperbolic systems. 

Lemma 7.1 allows us to estimate the integral of the terms |fi^fc|, \vjWk\ and \vjjWk\ in (7.3). 
We now work toward an estimate of the remaining terms \vj >x Vk\, \vj :X Wk\ and \vjWk, x \, containing 
one derivative w.r.t. x. 
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Lemma 7.3. Let z, be solutions of (7.4) and assume that (7.5) holds, together with the estimates 



J J \tp(t,x)\ dxdt < So, J J | yfi (t, x) | dxdt < So , 

N*)L> II^WL <*o, 

||A x (t)|| L i <C*5 , 
for all t £ [0, T]. Then one has the bound 



lim \(t,x) = 



J J \z x (t,x)\\z i (t,x)\dxdt = O(l)-5 2 . 



(7.12) 

(7.13) 
(7.14) 



(7.15) 



Proof. The left hand side of (7.15) is clearly bounded by the quantity 



T-T 



Z{T)= sup / / \z x {t,x)z\t + r,x + i)\dxdt < (C*5%) 2 -T, 
the last inequality being a consequence of (7.13). For t > 1 we can write z x in the form 

z x (t, x) = j G x (l,y)z(t - 1, x - y) dy + J j G x (s, y) [<p - (Xz) x ] (t - s, x-y) dyds , 
where G(t,x) = exp{— x 2 /4t}/2v / vrt is the standard heat kernel. Using (7.6) we obtain 

fT-T 



J J ' \z x (t,x)z i (t + T,x + £)\dxdt 



+ 



G x (l,y)z(t- 1, x-y) z%t + T,x + £) 

T-T ,■ ,-\ 




A 





T-t r r l 



J JO 
T-t r rt 



I L° 



L°° 



dydxdt 

G x (s, y) z(t — s, x — y) z^(t + t,x + £) dydsdxdt 
G x (s, y) z x (t — s, x — y) + r, x + £) dydsdxdt 

dydsdxdt 



+ J J J J\G x (t-s,x-y)<p(8,y)z i (t + T..r-( > ) 

< {^j \G x (l,y)\dy+ ||A X || L °° J J \G x (s,y)\dyd& 

• sup (f ( \z(t- s, x-y)\\z^(t + T,x + £)\dxdt\ 
s,y,T,i \Ji J J 

\G x (s,y)\ dyds) ■ sup / / \z x (t - s, x - y)\ \z*(t + r, x + £)| dxdt 
/ \s,y,T,$Ji J / 

■J j \G x (s,y)\dsdy ■ j J \<p(t,x)\dxdt 




+ PH L° 



< 2 • -25 -2-5 + 2||A|| L =oJ(T) + 0*8% ■ 25 . 



(7.16) 
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On the initial time interval [0, 1], by (7.13) one has 



/ / '\z x (t,x)z*(t + T,x + Z)\dxdt< I |Mi)|| L1 \\z\t + T)\\^dt < (C*5 2 ) 2 . (7.17) 
Jo J Jo 



Because of (7.14) we have 

2||A|| L oo < 2||A|| L i < 2C*5 « 1. 
From (7.16) and (7.17) it thus follows 

J(T) < (C*5l) 2 + — + 2C*5 ■ 1(T) + 2C*5l 



c 



For 5 sufficiently small, this implies T{T) < 95q/c, proving the lemma. □ 

Using the two previous lemmas we now prove the estimate (7.3). Setting z = Vj, z^ = Vk, 
A = Xj, A^ = Afc, an application of Lemma 7.1 yields the desired bound on the integral of |fjVfc|- 
Moreover, Lemma 7.3 allows us to estimate the integral of \vj tX Vk\- Notice that the assumptions 
(7.13)-(7.14) are a consequence of (5.22)-(5.23). The simplifying condition A(i,-oo) = in (7.14) 
can be easily achieved, using a new space coordinate x' = x — X*t. 

The other terms \vjWk\, \vjjWk\, \vj tX Wk\ and \vjWk, x \ are handled similarly. 



8 - Functionals related to shortening curves 

We now study the interaction of viscous waves of the same family. As in the previous section, 
let u = u(t,x) be a solution of the parabolic system (3.1) whose total variation remains bounded 
according to (7.1). Assume that the components Vi,Wi satisfy the evolution equation (6.1), with 
source terms c/)i,ipi bounded as in (7.2). Relying on the bounds (5.22)-(5.24), for each i = 1, . . . , n 
we shall prove the estimates 

J \wi, x Vi - WiV itX \ dxdt = 0(1) • <5q, (8.1) 

2 

dxdt = 0(1) • Si. (8.2) 




The above integrals will be controlled in terms of two functionals, related to shortening curves. 
Consider a parametrized curve in the plane 7 : M 1— > M 2 . Assuming that 7 is sufficiently smooth, 
its length is computed by 

£(7) = / \lx(x)\dx. (8.3) 
Following [BiB2], we also define the area functional as the integral of a wedge product: 

A{l) = \ll \l x {x)/\ lx (y)\dxdy. (8.4) 

1 J Jx<y 
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figure 8 



To understand its geometrical meaning, observe that if 7 is a closed curve, the integral 

\ J l{y) ^lx{y)dy = ^ j j 7a; (x) A 7a; (y) dx dy 

yields the sum of the areas of the regions enclosed by the curve 7, multiplied by the corresponding 
winding number (fig. 8). In general, the quantity -4(7) provides un upper bound for the area of 
the convex hull of 7. 

Let now 7 = ^y(t, x) be a planar curve which evolves in time, according to the vector equation 

It + A7a; = 7a;x • (8.5) 

Here A = A(i, x) is a sufficiently smooth scalar function. It is then clear that the length C(^(t)) 
of the curve is a decreasing function of time. It was shown in [BiB2] that also the area functional 
A{j(t)^j is monotonically decreasing. Moreover, the amount of decrease dominates the area swept 
by the curve during its motion. An intutive way to see this is the following. In the special case 
where 7 is a polygonal line, with vertices at the points Po, ■ ■ ■ , Pm, the integral in (8.4) reduces to 
a sum: 



-4(7) = ^J]h Av j|' 



2 



P - Pi 



i-i 



If we now replace 7 by a new curve 7' obtained by replacing two consecutive edges v/j, by one 
single edge (fig. 9b), the area between 7 and 7' is precisely A Vfe|/2, while an easy computation 
yields 

A(i) < A(j) - -\v h Av fc | . 

The estimate on the area swept by a smooth curve (fig. 9a) is now obtained by approximating a 
shortening curve 7 by a sequence of polygonals, each obtained from the previous one by replacing 
two consecutive edges by a single segment. 

We shall apply the previous geometric considerations toward a proof of the estimates of (8.1)- 
(8.2). Let v,w be two scalar functions, satisfying 



v t + (Xv) x -v xx = <j), 
w t + (Xw) x - w xx = tp . 

Define the planar curve 7 by setting 

l{t,x)=(l v(t,y)dy, [ w(t,y)dy). (8.7) 
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figure 9a 



figure 9b 



Integrating (8.6) w.r.t. x, one finds the corresponding evolution equation for 7 : 

7t + ^x-lxx=([ (f){t,y)dy, f ip(t,y)dy] . (8.8) 
\j — 00 J —00 J 

In particular, if no sources were present, the motion of the curve would reduce to (8.5). At each 
fixed time t, we now define the Length Functional as 



(8.9) 



C(t) = £{j(t)) = J ^v 2 {t,x)+w 2 (t,x) dx 
and the Area Functional as 

A(t) = A{j(t)) = \ [ I \v(t, x)w(t, y) - v(t, y)w(t, x) \ dxdy . (8.10) 

^ J Jx<y 

We now estimate the time derivative of the above functionals, in the general case when sources 
are present. 

Lemma 8.1 Let v,w be solutions of (8.3), defined fort G [0, T]. For each t, assume that the maps 
x 1 > v(t,x) x 1 ^ w(t,x) and x 1— > X(t,x) are C ' , i.e. continuously differentiable with Lipschitz 
derivative. Then the corresponding area functional (8.10) satisfies 

d f 

j t A(t) < -j v x (t,x)w(t,x)-v(t,x)w x (t,x) dar + ||u(t)|| Ll ||^)|| Ll + |m*)|| Ll ||0(t)|| Ll • (8.11) 

Proof. Differentiating (8.10) w.r.t. time we obtain 
dA _ 1 
~dt ~ 2 



J J sign(v(x)w(y) - w(x)v(y)) ■ {v t (x)w(y) + v(x)w t (y) - v t (y)w{x) - v(y)w t (x)}dxdy 
/ / sign(v(x)w(y) - w(x)v(y)) ■ I (v x (x) - \(x)v(x)) x w(y) - v(y)(w x (x) - X(x)w(x)) x \dxdy 

J Jx<y k J 

+ j j sign(v(x)w(y) — w(x)v(y)) ■ ^4>(x)w{y) + v(x)ip(y)^dxdy 
J j sign(v(x)w(y) - w{x)v(y)) ■ ^-^(v x {x) - \(x)v(x))w(y) - v(y)(w x (x) - \{x)w(x))}dxdy 
+ 11 s\gn(v(x)w(y) - w(x)v(y)) ■ I <f>(x)w(y) + v(x)ip(y) >dxdy . 

J Jx<y L J 



(8.12) 
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To simplify the first term on the right hand side of (8.12), for a fixed y we let be the set 

of points < y such that v(x a )w(y) — w(x a )v(y) = 0. Relying on an approximation argument, we 
can assume that these points are in finite number, say x/v(y) < • • • < x 2 {y) < X\(y) < Xo(y) = y, 
and that the function x i— > v(x)w(y) — w{x)v{y) changes sign across each x a (y). For convenience, 
we define the additional point x^+i = — oo and let i(y) = ±1 be the sign of v(x)w(y) — v(y)w(x) 
on the last interval, i.e. when x 6 [x±, y]. We now can write 

J sign(v(x)w(y) - w{x)v(yj) ■ — \{v x (x) - X(x)v(x))w(y) - v(y)(w x (x) - X(x)w(x))jdx 



i(v)-£(-i)' 



dx 



x a +i(y) 



{(v x (x) - X(x)v(x))w(y) - v(y)(w x (x) - \(x)w(x))}dx 
2i (y) " I 5^(- 1 ) a (^( a; a(y))^(y) - w x {x a {y))v{y)^ + {y x (y)w(y) - w x (y)v{y)) J . 



By direct inspection, one checks that the factor i(y) ■ (—1)° has always the opposite sign of 
v x (x a )w(y) - w x (x a )v(y). Therefore 

f sign(v(x)w(y) — w(x)v(y)) • — |(v x (ar) — X(x)v(x))w(y) — v{y)(w x (x) — X(x)w(x))^dxdy 

x<Cy 



N 



a=l 



v x (y)w(y) - v(y)w x (y)\dy - 2^ / \v x (x a (y))w(y) - v(y)w x {x a (y)) 



dy. 



The bound (8.11) is now an immediate consequence of (8.12) and (8.13). 



(8.13) 
□ 



Lemma 8.2. Together with the hypotheses of Lemma 8.1, at a fixed time t assume thatj x (t,x) ^ 
for every x. Then 



d . . 1 

Jt {t) ~ "(1 + 9<5?)3/2 J lw/v{ < 3Si 



V(t)\ 



W(t) 
V(t) 



dx+||0(t)|| L1 + ||^)|| L1 . (8.14) 



Proof. As a preliminary, we derive the identity 

hxx\ 2 \lx\ 2 - <7x, Txx) 2 = (vl + wl x )(v 2 +w 2 x ) - (vv x +w x w xx ) 2 

= {vw xx +v x w x ) 2 = v i \(w/v) x \ 2 . 
Thanks to the assumption that j x never vanishes, we can now integrate by parts and obtain 

\1xt ^ixt) 



\J(lx, Ix) 



dx 



\7x5 ^YxXX > <7x, (X-f x ) x )_ (jx, 1 dx 



hx 



hx 



= J | \lx\xx - (A 17*1),. - 



|7a;a:| (7a: /1 7a; |> "Ixx) 
\lx\ 



+ 



\lx\ 
(lx, (0,Y>)) 



hx 



dx 



< 



\v\ \(w/v) x \ 



■dx + wmwv + wmwv 



(i + { w / v yf /2 

Since the integrand is non-negative, the last inequality clearly implies (8.14). 
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□ 



Remark 8.3. Let u = u(t, x) be a solution to a scalar, viscous conservation law 



ut + f{u) x -u xx = 0, 

and consider the planar curve 7 = (u, f(u) — u x ) whose components are respectively the conserved 
quantity and the flux (fig. 10). If A = /', the components v = u x and w = — ut evolve according to 
(8.6), with = ^ = 0, hence jt + \~/ x — ~j xx = 0. Defining the speed (fig. 11) s(x) = —ut(x)/u x (x), 
the area functional ^(7) in (8.4) can now be written as 

^(7) = ^ / / \u x (x)u t {y) -u t (x)u x (y)\dxdy 

^ J Jx<y 

= \ i i \u x (x) dx\ ■ \u x (y) dy\ ■ \s(x) - s(y)\ 

1 J Jx<y 

= - [wave at x] x [wave at y] x [difference in speeds] . 

2 J J x <Cy 

It now becomes clear that the area functional can be regarded as an interaction potential between 
waves of the same family. In the case where viscosity is present, this provides a counterpart to the 
interaction functional introduced in [L2] in connection with strictly hyperbolic systems. 




x a b c d 

figure 10 




x y 

figure 11 
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Recalling that the components Vi,Wi satisfy the equations (6.1), we can apply the previous 
lemmas with v = Vi, w = Wi, A = Aj, cj> = fa, ip = ipi, calling Li and Ai the corresponding length 
and area functionals. For t G [t, T], the bounds (5.22)-(5.23) yield 



A(t)<||^(t)|| LOO -||^(t)|| L1 =o(i)-5 3 , 
Ci(t) < ||^(t)|| L1 + \\wi(t)Li =o(i)-s . 



(8.15) 
(8.16) 



Using (8.11) we now obtain 



\ w i x v i ~ WiVi x dxdt < 



/ (||ui(t)|| L i||V'i(*)|| L i + lhi(*)ILi||^(*)ILi) dt 



<Ai{t)+ sup (||ui(t)|| Ll + |hi(*)|| L i) • f j (\fa(t,x)\ + \^(t,x)\)dxdt 



te[t,T] 
= 0{l)-5l, 

(8.17) 

proving (8.1). To establish (8.2), we first observe that, by an approximation argument, it is not 
restrictive to assume that the set of points in the t-x plane where Vi )X {t,x) = Wi yX {t,x) = is at 
most countable. In this case, for almost every t 6 [t, T] the inequality (8.14) is valid, and hence 



t J | Wi /vi I <3<5i 



dxdt < 



dt 



Mt)\\ L i + \\Mt)\\ L i) 



dt 



< Ci(t) + 
= 0(1)- Jo- 
Using the bound (5.23) on 1 1 v 1 1 1 loo ? from (8.18) we deduce (8.2). 



\fa{t,x)\ + \^{t,x)\) 



dxdt 



(8.18) 



9 - Energy estimates 

In the same setting as the two previous sections, we shall now prove the estimate 

J (Kx| + K,x|) \w i -6 i v i \dxdt = 0(l)-5%. (9.1) 

We recall that Oi = 9{wi/vi), where is the cut-off function introduced in (5.5). Notice that the 
integrand can be ^ only when \wi/v i \ > S\. 

Consider another cut-off function rj : iR i — > [0,1] such that (fig. 5) 

*>-{! 2 ™ 

We can assume that rj is a smooth even function, such that 

|»/| < , \v"\<m/6f. 
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For convenience, we shall write r/j = r](wi/vi). As a preliminary, we prove some simple estimates 
relating the sizes of Vi, Wi and t^. It is here useful to keep in mind the bounds 

A, - A* =0(1)- \ri - r*\ = 0(1) ■ 5 , \ Vi \, M = 0(1) ■ 5 2 , (9.3) 

valid for t > i and i = 1, . . . , n. Recall also our choice of the constants 

< 5 « S t < i . (9.4) 

Lemma 9.1. If\wi/vi\ > 3<5i/5, i/ien 

KI<2K Jr |+O(l)-<y X)l u il' \vi\<^\vi, x \+O(l)-S J2\ v j\- ( 9 - 5 ) 

.J A 3& 

On the other hand, if\wi/vi\ < 4#i/5, £/ien 

k*l <5iH + e>(iMoJ>;l- ( 9 - 6 ) 

Proof. We recall the first estimate in (6.16): 

Vi, x = Wi + (Xi - X*)vi + G 

Wi + {Xi - X*) Vi + 0(1) ■ 5 N » 



(9.7) 



with 6 defined as in (6.15). By (9.3)-(9.4), from the condition \wi/vi\ > 3e>i/5 two cases can arise. 
On one hand, if 

l©l<^H, (9-8) 

then 

\Vi,x\ > —\ V i\ - • °0\ v i\ - ^qN ^ ~5~l Wi l ' 

and hence 

5 - i i i i i 

K ^ 77T K* > Kl ^ + tK < 2 k*l ■ (9.9) 

zoi 5 

On the other hand, if (9.8) fails, then by (9.3) we have 

\~X i -X*\\v i \=O(l)-S 2 J2\vj\- (9.10) 



In both cases, the estimates in (9.5) hold. 

Next, if \wi/vi\ < 4(5i/5, from (9.7) we deduce 

K x \ < ^\vi\ + ^\vi\ + O(l)-S J2 M- ( 9 - n ) 

If (9.8) holds, then \vi x \ < 5i\vi\. If (9.8) fails, then (9.10) is valid. In both cases we have (9.6). 

□ 
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Toward a proof of the estimate (9.1), we first reduce the integrand to a more tractable expres- 
sion. Since the term \wi — vanishes when \wi/vi\ < Si, and is < \wi\ otherwise, by (9.5) we 
always have the bound 



if, 



- 0%Vi\ < \ViWi\ < Vi 2\vi, x \ + 0(1) • S 



Therefore 



(\Vi, x \ + K,x|) ' K - °i V i\ ^ (\ V i,x\ + \Wi,x\)Vi 2 \Vi,x\ + 0(1) ■ S ^ \Vj 



< 2 ViVi, x + 2Vi\Vi, X W i>x \ + ^2 {\ V 3 V i,x\ + \ V 3 W i,x\) 



(9-12) 



Since we already proved the bounds (7.3) on the integrals of transversal terms, to prove (9.1) we 
only need to consider the integrals of vf x and wf x , in the region where rji ^ 0. In both cases, 
energy type estimates will be used. 



We start with vf x . Multiplying the first equation in (6.1) by rjiVi and integrating by parts, 
we obtain 

/ rjiV^i dx = {riiViVij + w(A;«;) x - ViViVi, xx } dx 

= J { (ViVi/2) t + CX l T]i)x(Vi/2) - (m jt + 2KVi,x ~ Vi,xx)(vf/2) + r^vf ^ dx . 



Therefore 



J Vivl x dx = J 7]ivf/2dx + J (rii,t + KVi,x-r]i,xx)(vi/2)dx 



- / Xi,xVi(Vi /2) dx + / rjiVifadx. 



(9.13) 



A direct computation yields 



/ / Wi,t V iyt Wi \ ~ , ( w i)X V itX Wi 

Vi,t + XiVi,x ~ Vi,xx = Vi 2~ + XiVi 2~~ 

Vi Vf / V Vi vf 



-Vi 



if, 



Vi,xxWj _ Vj^xWi^x V i,x W i 
2 2 ' 3 

vf vf vf 



Vi{vJi,t + {Kwi) x - w iyXX )/vi - rjiWi(v it t + (XiVi) x - v iiXX )/v. 

+ 2v itX T]'Jv i ■ (Wi/Vi) x - Vi{Wi/Vi)l 
, ( 1/Ji Wi<f>i\ , n , Vi, X ( W. 



(9.14) 



Vi Vi Vi 



'I; — — — ' - > '/,— — 



r, \<k/ x 1 \v iy x 
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Since Aj jX = (\ — X*) x , integrating by parts and using the second estimate in (9.5) one obtains 
~Xi, x Vi(vf/2)dx = 



(Xi - X*)(r]i :X Vi/2 + T]iViV i:X ) dx 

5 



\v'i\ \wi, x Vi - v itX Wi\ dx + 



2^ / Vi v i,x dx + ' 5 ° JYI \viVj\dx \ 

j¥=i ) 



< I \Wi,xVi ~ V itX W, 



I dx + ^ J r]iV 2 i x dx + 5 J \ v i v o I dx ■ 



(9.15) 



Indeed, by (9.3)-(9.4), \\ - A*| = 0(1) • 5 « Si. Using (9.14)-(9.15) in (9.13) we now obtain 



1 f 2 , d 

zjmv^dx < -- 

2 



dx 



+ 



Wi 

r}iViVi, x I — 



1 

+ 2 



// 2 / w i 



dx + / \Wi. x Vi - w. 



\ J KKKV'tl + K&l) dx- 

t> i>a; | dx + <5 J ^\viVj\dx + J\vi4>i\dx. 



dx 



(9.16) 

Recalling the definition of rji, on regions where »]• ^ one has < 4<5i/5, hence the bounds 

(9.6) hold. In turn, they imply 



ViViVi, x 



\ V i J x 

< \SWi\ \w i:X Vi - WiV i:X 



ViViVj 



VjW iyX \ + \VjV iy 



Using the bounds (5.22)-(5.23). (7.2)-(7.3) and (8.1)-(8.2), from (9.16) we conclude 



(9.17) 




rjiV ix dxdt< / r]iV i (t, x) dx + 0(1) 




+ \wi4>i\) dxdt 

+ 0(1) • J_ j \w itX Vi - WiV itX \ dxdt + 0(1) ■ 5 J y~] (|v 3 -«; i)a ;| + l^-tvl) 

+ 0(1)- / / \vi(wi/vi) x \ 2 dxdt + 2<5 / / l^t^jl dxdt + 2 / / dxdt 

Ji J\w i /v t \<S 1 Jt J ... Ji J 



0(1) • 51 



(9.18) 



We now perform a similar computation for wf x . Multiplying the second equation in (6.1) by 
r/iWi and integrating by parts, we obtain 

J rjiWi^i dx = J ^{r] i w 2 i /2) t + {XiT]i) x {w 2 /2) - (rj ijt + 2X^1]^ - r) itXX ) (w 2 /2) +^10^} dx . 

Therefore, the identity (9.13) still holds, with Vi,(f>i replaced by Wi,ipi, respectively: 



• j d 



ijiW 2 /2dx 



+ / {m,t + \Vi,x - Vi,xx) {w 2 /2) dx 



- j X itX r]i(w 2 /2)dx + J riiWiipi dx . 



(9.19) 
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The equality (9.14) can again be used. To obtain a suitable replacement for (9.15), we observe 
that, if rji ^ then (9.5) implies 



\WiWi,x\ < 2\Vi tX W. 

Integrating by parts we thus obtain 
Xi, x rn{w1/2)dx 



Sq ^ \VjWi,x\ < «i,x + W i,x + ■ ^0 \ v 3 w i,x\ ■ 



(A; - X*)(r] itX Wi/2 + r)iWiW itX ) dx 

< \\\ - A*|| L ~ • \ / Wi\ \Wi,xVi - v i:X w 



dx+ / Vivl x dx 



+ / ViW^dx + Oil) ■ 5 / ^2\vjW ijX \dx 
J J &i 

\ dx+ \ j 7 li v i, x dx + ^ J 7]iwl x dx + 5 J Y^\ v j w i,x\dx , 



< / \Wi. x Vi - Vi. x W 



(9.20) 

Using (9.14) and (9.19) in (9.18) and observing that \wf/vf \ < 5\ on the region where 77- 7^ 0, we 
now obtain an estimate similar to (9.16): 



1 f 9 , d 

2 h iWi - dx --jt 



[J 


' dx 


K 






2 








11 2 ( 






/ 


ViVi [ 


, V i)x 



+ y/ + \wi<f>i\) dx + ^ J 



dx + I \wi x Vi — WiVi x \ dx 



ViViVi, x 



dx 



VjWi !X \dx + J\Wi4>i\dx. 

(9.21) 

Using the bounds (5.22)-(5.23), (7.2)-(7.3), (8.1)-(8.2) and (9.17)-(9.18), from (9.21) we conclude 
J J r]iW^ x dxdt< J r]iWi(i,x) dx + 0(1) ■ J J (\viipi\ + \wi4>i\) dxdt 

+ 0(1) ■ / / \w iyX Vi - WiV iyX \ dxdt + O(l) -8 Y {\ v j w i,x\ + \vjVi, x \) dxdt 

Ji J Ji J 

+ 0(1)- / / \vi(wi/vi) x \ 2 dxdt + / r)iV? x dxdt 

Ji J\w i /v i \<5 1 Ji J 

+ -So J Y2 \ w i,xVj\ dxdt + 2 j \u)iipi\ dxdt 
= O(l)-5 2 . 

(9.22) 

Using (9.18) and (9.22) in (9.12), we obtain the desired estimate (9.1). 
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10 - Proof of the BV estimates 



In this section we conclude the proof of the uniform BV bounds. Consider any initial data 
u : M i-> M n , with 

r 

Tbt.Var.fu} < — L- , lim u(x) = u* € K . (10.1) 

We recall that k is the constant defined at (2.5), related to the Green kernel G* of the linearized 
equation (2.4). This constant actually depends on the matrix A(u*), but it is clear that it remains 
uniformly bounded when u* varies in a compact set K C M n . 

An application of Corollary 2.4 yields the existence of the solution to the Cauchy problem 
(1.10), (1.2) on an initial interval [0,i], satisfying the bound 

K(t)|| L1 <^. (10.2) 
This solution can be prolonged in time as long as its total variation remains small. Define the time 

r = sup|r; ^£ J\<j>i(t,x)\ + \4i(t,x)\dxdt< j\ . (10.3) 

If T < oo, a contradiction is obtained as follows. By (5.21) and (10.2), for all t 6 [t, T] one has 

lh(i)|| L1 <Elh(')L 

i 

^E +J. J ' \Mt,x)\dx(ttj (10.4) 

< 2v^||«x(t)|| L i +y < ho- 
using Lemma 6.1 and the bounds (7.3), (8.1), (8.2) and (9.1) we now obtain 

E/ I \Mt,x)\ + \Mt,x)\dxdt = O(l)-5 2 < |, (10.5) 



provided that 5q was chosen suitably small. Therefore T cannot be a supremum. This contradiction 
with (10.3) shows that the total variation remains < So for all t G [t, oof . In particular, the solution 
u is globally defined. 

Remark 10.1. The estimates (8.1) and (9.1) were obtained under the assumption (7.2) on the 
source terms. A posteriori, by (10.5) the integral of the source terms is quadratic w.r.t. So- Using 
(10.5) instead of (7.2) in the inequalities (8.17) and (9.18), (9.22), we now see that the quantities 
in (8.1) and (9.1) are both = 0(1) • <5q. Recalling that So is the order of magnitude of the total 
variation, we see here another analogy with the the purely hyperbolic case [G]. Namely, the total 
amount of interactions between waves of different families is of quadratic order w.r.t. the total 
variation, while the interaction between waves of the same family is cubic. 
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Remark 10.2. Within the previous proof, we constructed wave speeds = A* — 6(wi/vi) for 
which the following holds. Decomposing the gradients u x ,u t according to 



u x = y~]vi fi(u,Vi,ai), 

i 

U t = ^{Wi ~ K V i) ?i(u, Vi, Oi), 



(10.6) 



the components Vi,Wi then satisfy 



Vi,t + (Kvi) x ~ v iyXX = <f>i(t, x) , 
w i>t + (XiWi) x - w i:XX = ipi(t, x) , 



(10.7) 



where all source terms <pi,ipi are integrable: 

bi(t, x) | dxdt < 5q , 



\tpi(t, x) \ dxdt < 5q 



(10.8) 



In general, the speeds Ui defined at (5.7) are not even continuous, as functions of t,x. However, 
by a suitable mollification we can find slightly different speed functions ai(t,x) which are smooth 
and such that the corresponding decomposition (10.6) is achieved in terms of (smooth) functions 
Vi,Wi satisfying a system of the form (10.7), with source terms again bounded as in (10.8). 

We conclude this section by studying the continuous dependence w.r.t. time of the solution 
t ' ^ u(t, •)• By (10.4) we have 



Tot.Var.{u(t)} = ||«x(*)|| L i < $o for all t > 0. 

By the estimate (2.8) in Proposition 2.1, the second derivative satisfies 

II xx\ ,|| L i - \ 2k8g i^I if t > t 

Therefore, from (1.10) it easily follows 

||u t (t)|| Ll <L'(l + ^), 
for some constant V . For any t > s > we now have 

IK*) - n ( s )ll L i - / IM t )IIli dr 

J s 

< L'(\t-s\ + \Vi- y/s\j. 



(10.9) 



(10.10) 



;io.n) 



Remark 10.3. A more careful analysis shows that in (10.11) one can actually take V = 0(1) • 
Tot.Var.jn}. However, this sharper estimate will not be needed in the sequel. 
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11 - Stability estimates 



Let u = u(t,x) be any solution of (3.1) with small total variation. The evolution of a first 
order perturbation z = z(t,x) is then governed by the linear equation 



zt + [A{u)z) 



(u x • A(uj)z — (z • A(u)) 



11, 



(11.1) 



As usual, by "•" we denote a directional derivative. The primary goal of our analysis is to establish 
the bound 

||*(*,-)|| L i < L||z(0,-)|| L i foralli>0, (11.2) 

for some constant L. By a standard homotopy argument [Bl], [BiBl], this implies the uniform 
stability of solutions, w.r.t. the L 1 distance. Indeed, consider two initial data u,v with suitably 
small total variation. We can assume that u* = u(— oo) = v{— oo), otherwise \\u — = oo and 
there is nothing to prove. We construct the smooth path 

6^u e = 6u + (l-6)v, 0€[0,1]. 

Calling 1 1 > u d (t, •) the solution of (3.1) with initial data u 6 , for every (>0we have 



\u(t)-v(t)\\ hl < [ 
Jo 



< L ■ 



du 9 (t) 




dO 


L 1 


■ i 


du°(0) 


) 


de 


l - 


■ V\\-L 


L . 



de 



de 



(11.3) 



Indeed, the tangent vector 



du® 



is a solution of the linearized Cauchy problem 

4 + [DA(u 6 ) • z e ]u e x + A(u e )z e x = z e xx , 
z d (0,x) = z e (x) = u(x) — v(x), 



hence it satisfies (11.2) for every e. The bound (11.3) provides the Lipschitz continuous dependence 
of solutions of (3.1) w.r.t. the initial data, with a Lipschitz constant independent of time. In 
particular, it shows that all solutions with small total variation are uniformly stable. 

Remark 11.1. In the hyperbolic case, a priori estimates on first order tangent vectors for solutions 
with shocks were first derived in [B2]. However, even with the aid of these estimates, controlling 
the L 1 distance between any two solutions remains a difficult task. Indeed, a straightforward use of 
the homotopy argument fails, due to lack of regularity. These difficulties were eventually overcome 
in [BC] and [BCP], at the price of heavy technicalities. On the other hand, in the present case with 
viscosity, all solutions are smooth and the homotopy argument goes through without any effort. 

Throughout the following, we consider a reference solution u = u(t, x) of (3.1) with small total 
variation. According to Remark 10.2, we can assume that there exist smooth functions Vi,Wi<Ji for 
which the decomposition (10.6) holds, together with (10.7) and (10.8). 
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The techniques that we shall use to prove (11.2) are similar to those used to control the total 
variation. By (2.19) we already know that the desired estimate holds on the initial time interval 
[0,t\. To obtain a uniform estimate valid for all t > 0, we decompose the vector z along a basis of 
unit vectors and derive an evolution equation for these scalar components. At first sight, it looks 
promising to write 

z = ^2zifi(u,Vi,ai), 

i 

where fi, . . . , f n are the same vectors used in the decomposition of u x at (5.6). Unfortunately, this 
choice would lead to non-integrable source terms. Instead, we shall use a different basis of unit 
vectors fi, . . . , r n , depending not only on the reference solution u but also on the perturbation z. 
Toward this decomposition, we introduce the variable 

T = z x - A(u)z, 

related to the flux of z. By (11.1), this quantity evolves according to the equation 

T t + (A(u)Y) x — T xx = (u x • A(u))z - (z • A(u))u x - A(u) (u x • A(u))z - (z • A{u))u 

+ (u x »A(u))T-(u t »A(u))z. 

We now decompose z, T according to 

z = ^2h i f i (u,v i ,X* - e(gi/hi)) , 

i 



(11.4) 



where 9 is the cutoff function introduced at (5.5). In the following we shall write 

h = fi(u,Vi, X* - Oigi/hi)), 

to distinguish these unit vectors from the vectors fj(-u,t»j, A* — 6{wi/vij) previously used in the 
decomposition (5.6) of u x . Moreover we introduce the speed 



and denote by 



Xi = (h, A(u)fi) 



= B(9i/hi 



(11.5) 



(11.6) 



the correction in the speed for the perturbation. The next result, similar to Lemma 5.2, provides 
the existence and regularity of the decomposition (11.4). 

Lemma 11.2 Let \u — u*\ and \v\ be sufficiently small. Then for all z,T 6 M n the system of In 
equations (11-4) has a unique solution (hi, . . . , h n , gi, . . . , g n ). The map (z, T) ^ (h, g) is Lipschitz 
continuous. Moreover, it is smooth outside the n manifolds Mi = {hi = gi = 0}. 

Proof. The uniqueness of the decomposition is clear. To prove the existence, consider the mapping 
A : M 2n t— > M 2n defined by 



(11.7) 
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K(hi,9i) = 

Computing the Jacobian matrix of partial derivatives we find 



hifi(u, Vi, X* -9{gi/hi)) \ ( 
(9i - \*hi) fi(u, Vi , X* - Oigi/hi)) )' [ ' ' 



d(hi, gi ) \-X*fi - Xtisi/hi)^ + (gi/hife'th,* h + X*9[h,, - (9i/hiWr 



(11.9) 



By (4.24), rj )(7 = O(l) ■ Wj. Hence, for Vi small enough, the differential DA is invertible. By 
the implicit function theorem, A is a one-to-one map whose range covers a whole neighborhood 
of the origin. Observing that A is positively homogeneous of degree 1, we conclude that the 
decomposition is well defined and Lipschitz continuous on the whole space M 2n . Outside the 
manifolds Mi, i = 1, . . . , n, the smoothness of the decomposition is clear. □ 

Writing the identity T = z x — A(u)z in terms of the decomposition (11.4) we obtain 

^2(9i ~ K h i)?i = ^2 hi > x ^ ~ A ( u ) h i r i + ^2 hifi,u Vjfj + ^ hif i:V V i>x - ^ hiVi,a6i,x ■ 
i i i ij i i 

(11.10) 

Taking the inner product with r, and observing that fi is a unit vector and hence is perpendicular 
to its derivatives, we obtain 

9i = K,x + (A, - X*)hi + 

with 

@ = ^2(fi, A{u)rj)hj + ^2^2(h, r jtU r k )h jVk + ^ (f< , f jtV )h jVjtX + ^ (f t , r^^A* ■ 

(11.11) 

Hence, by (5.20) and (4.24), 

<?* = hi, x + (Xi ~ K)hi + 0(1) ■ 5 J2 (N + N) ■ ( 1L12 ) 

A straightforward consequence of (11.12) is the following analogue of Lemma 9.1. 
Corollary 11.3. If\gi/hi\ > 3<5i/5 ; then 

\9i\ <2\h itX \+0{l)-5 Q y^{\v j \ + \h j \), \hi\ < Al^l+O^-^oVd^l + l^-l). (11.13) 

On the other hand, if\gi/hi\ < 4<5i/5 ; then 

|^|<<5i|^|+C»(l)-<5o^(K| + |^|). (11.14) 

Our eventual goal is to show that the components hi,gi satisfy a system of evolution equations 
of the form 

{hi : t + (Xihi) x — hi^ xx = 4>i, 
(11.15) 
9i,t + (Xi9i) x - gi , xx = ipi , 
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where the source terms on the right hand sides are integrable on [t, oof xM. Before embarking 
in calculations, we must first dispose of a technical difficulty due to the lack of regularity of the 
equations (11.4). 

Since our equations (3.1) and (11.1) are uniformly parabolic, it is clear that for t > all 
solutions are smooth. Moreover, by Remark 10.2, we can slightly modify the speeds o\ occurring 
in the decomposition of u x , so that (10.6)-(10.8) hold and the corresponding functions V\ are now 
smooth. On the other hand, the map A in (11.7) is only Lipschitz continuous, hence the same is 
true in general for the functions hi = hi(t, x) and g\ = gi(t, x). Indeed, at points where hi = gi = 
for some index i, the derivatives hi iX or g ijX may well be discontinuous. In this case, the equations 
(11.15) would make no sense. To avoid this unpleasant situation, we observe that each manifold Mi 
has codimension 2. Given the smooth functions z, T and e > 0, by an arbitrarily small perturbation 
we can construct new functions /i",T" satisfying 

\\z* -z\\ C 2 + ||T«-T|| C 2 < e 



and such that the corresponding decomposition (11.4) is C°° outside a countable set of isolated 
points (t m ,x m ) m >i. A further implementation of this technique yields 

Lemma 11.4. Let z, T be solutions of (3.1) and (11.1) respectively. Then for any e > there 
exists smooth functions z$ , T" such that the corresponding coefficients in the decomposition (11-4) 
are smooth except at countably many isolated points (t m ,x m ), m > 1. Moreover, these perturbed 
functions solve the system of equations 

z\ + (A(v)z*) x - zl x = (u x • A(uj)z* - {z» • A(u))u x + ei (t,x) , 
Tf + (A(«)T») t - T« x = \(u x . A(u))z* - (z* • A{u))u x ] - A(u) \(u x • A(u))z* - (z* • A(u)) 

L J x L 

+ (u x • A(u))r* - (u t • A(u))z* + e 2 (t,x) , 



for some perturbations e±, e 2 such that 

J J \ei(t, x) | + \e<z{t, x) | dxdt < e . 



Thanks to this lemma, we can study the time evolution of the components hi,gi by means 
of a second order parabolic system, at the price of an arbitrarily small perturbation on the right 
hand side. In the remainder of the paper, for simplicity we derive all the estimates in the case 
e\ = e 2 = 0. The general case easily follows by an approximation argument. 

In Section 6 we showed that the source terms in the equations (6.1) could be reduced to four 
basic types. The following result is an analogue of Lemma 6.1, providing an estimate for the 
source terms in the equations (11.15). The proof, involving lengthy calculations, will be given in 
Appendix B. 
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Lemma 11.5. The source terms in the equations (11.15) satisfy the estimates 

4>i(t,x), ipi(t,x) = 0(1) • (\ h i,x\ + \hjVj\ + \gjVj\ + \g jtX \) \wj - 6jVj\ 

j 

+ ■ Yl {\ V 3 k 3,x ~ h i V lA + \ V 3,x9j ~ 9j,xVj\ + \hjW jtX - WjhjJ + \9jW jtX - gj, x Wj\} 

3 

+ o(i)-Etol + N)|»*(^)]-x {h , /kjl<Ml } 

+ C> ( 1 ) ' Y {\ h l v k\ + \hj, x Vk\ + \hjV ktX \ + \hjW k \ + \gjV k \ + \g x ,jV k \ + \gjV kjX \ + \h 3 h k \ + \hjg k \) 

+ o(i) ■ (\ h 3<t>3\ + Mil + M;l + MA)- 



The key step in establishing the bound (11.2) is to prove 
Lemma 11.6. Consider a solution z of (11.1), satisfying 

|Kt)|| L i<<*0 for allt € [0,T], 

and assume that the source terms in (11-4) satisfy 



i(t,x) \ + \4>i(t,x) \ dxdt < 6 



i = 1, . . . , n. 



Then for each i = 1, . . . , n one raas i/ie estimates 

r 



(11.16) 



(11.17) 



(11.18) 



J \4>i(t,x)\dxdt = 0(l)- Sl, J J\$i(t,x)\dxdt = 0(l)-6%. (11.19) 



Assuming the validity of this lemma, we can easily recover the estimate (11.2). Indeed, since 
the equations (11.1) are linear, it suffices to prove the estimate in the case where 



(11.20) 



L1 8^k' 



We recall that k is the constant defined at (3.5). By Corollary 2.4, on the initial interval [0,t\ we 
have 

||*(t)|| Ll < 2k||*(0)|| l1 = A= t€[0,t\. (11.21) 



Define the time 



T = sup I t ; ^ 




(5 

i(t,x) \ + \4>i(t,x) \ dxdt < — 



(11.22) 
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If T < oo, a contradiction is obtained as follows. First, we observe that the inequalities in (5.17) 
remain valid for the decomposition of z, namely 



\z\ < \ h i\ ^ %>/n\z\ . (11.23) 

i 

For every r € [t, T], by (11.22) and (11.23) one has 

<EIM r )L 

i 

^5Z(|M*)Hli+J[ J\4>i(t,x)\dxdtj (11.24) 
<2v^||^)|| L1 +y < do. 
We can thus use Lemma 11.5 and conclude 

Y.C I \Ut,x)\ + \i>i(t,x)\dxdt = O(l)-6 2 < |, (11.25) 

provided that So was chosen suitably small. Therefore T cannot be a supremum. This contradiction 
shows that the bound (11.2) holds for all t > and z G L 1 , with L = 8K^/n. The remainder of 
this section is aimed at establishing the estimates (11.19). 

Proof of Lemma 11.6. By Corollary 2.2, for t G [i,T], as long as ||z(t)|| L1 < So we also have the 
bounds 

||2x(f)|| Ll = O(l) ■ 51 , ||^(t)|| Ll = O(l) ■ Si , \\z xx (t)\\ hoc = 0(1) ■ Si . 

By Lemma 11.1, the map (z, T) t— > (h,g) is uniformly Lipschitz continuous. From the previous 
bounds, for every t 6 [t, T] and all j = 1, . . . , n it thus follows 

||^>W|| L i. IM*)|| L i, IM*)|| L ~> IM*)|| L ~ = 0{l)-8l, (11.26) 

||^,xW|| L =o, IM*)|U = 0{l)-Sl. (11.27) 



Recalling that Vi,Wi,hi,gi satisfy the systems of equations (10.7) and (11.15) with source terms 
bounded by (10.8) and (11.18), we now provide an estimate on the integrals of all terms on the 
right hand side of (11.16). 

The same techniques used in Section 7 yield an estimate on all transversal terms, with j ^ k: 
j (\hjV k \ + \h jtX v k \ + \hjV ktX \ + \hjW k \ + \gjV k \ + \g j>x v k \ + \gjV kiX \ + \hjh k \ + \hjg k \) dxdt 

(11.28) 
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0(1) • sl . 



From (10.8) and (11.26) one easily obtains 
I J + Mi | + M ,| + < ^ / (||Ml~ + HftllL-) • (10,1 + KM) ** 

(11.29) 

A further set of terms will now be bounded using functionals related to shortening curves, as 
in Section 8. At each fixed time t € [t, T], for i = 1, . . . , n consider the curves 



Vi(t,y)dy, / hi(t,y)dy 



(v,h), \ 



By obvious meaning of notations, we also consider the curves ^-"' g \ jj, w ' h \ l\ w,9 \ l\ h ' 3 ^ • By (6.1) 
and (11.15), the evolution of these curves is governed by vector equations similar to (8.8). For 
example, 

J v ' h '> + \^ V 'V _ J v > h ) 



f nX f-X \ 

/ 4>i{t,y)dy, I Mt,y)dy 

\J — oo J — CO J 

As in (8.9)-(8.10), we introduce the corresponding Length and Area Functionals, by setting 

C^ h \t) = C{^ h \t)) = J y/v*{t t x) + h*{t,x) dx, 

Af' h \t)=A{ 1 ( t' h \t))= \ I [ \v i (t,x)h i (t,y)-v i (t,y)h i (t,x)\dxdy. 

* J Jx<y 

Similarly we define c[ v ' 9 \t), A[ v ' 9 \t), etc... A computation entirely analogous to (8.17) now 
yields the bounds 



J J (\vihi !X - hiV ijX \ + \v ijX gi - g i>x Vi\ + \hiW ijX - Wih i}X \ 

+ \9iWi, x ~ 9i,xWi\ + \gih i>x - hig itX \\dxdt = 0(1) • 5% . 
Moreover, repeating the argument in (8.18) we obtain 

alxdt = 0(1) • 5 . 

Using the bounds (5.19) on H^Hl 00 and (11.26) on || /i-i , from (11.31) we deduce 



(11.30) 







n 




It J\w i /vi\<3S 1 





(11.31) 



[ T [ (N + N) 

Jt J\w i /v i \<3S 1 



hi 



alxdt = 0(1) • 5% . 



(11.32) 



The integrals of the remaining terms in (11.16) will be bounded by means of energy estimates. 
For convenience, we write rji = r/(gi/hi), where r\ is the cutoff function introduced at (9.2). In 
Appendix C we will prove the estimates 



Jo J 



dx = 0(1) • 5l 



(11.33) 
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j* J gl x fi l dx = O(l)-5 2 . 



(11.34) 



Using (11.33)-(11.34) we now bound the terms containing the "wrong speed" \wi — 6iVi\. All these 
terms can be 7^ only when \ wi/vi\ > Si. Hence by (6.18) we can write 



{\hiVi\ + \giVi\) \wi - 9iVi\ = O(l) • (\hi\ + \gi\) ( v i>x (wi - 6iVi)\ + ^ \vj (wn - 9iV 

By (7.3), (9.1) and (11.26) it thus follows 

J J {\hiVi\ + \giVi\) \w i -e i v i \dxdt = 0(1) -<^. (11.35) 

To estimate the remaining terms, we split the domain according to the size of \gi/hi\. 
CASE 1: \gi/hi\ > 4<5i/5, \v)i/vi\ > Si. Recalling (9.5) we then have 

{\hi, x \ + \9i, x \) K ~ 0%Vi\ < {\h iyX \ + \g itX \) K| 

= (\hi, x \ + \gi, x \) (2\vi, x \ + 0(1) ■ \ v j\) 

< {hi, x + 9\ x + 2vl x ) + 0(1) • £ KxVjl + 0(1) • £ 
Using (11.33)-(11.34), (9.17) and (11.28), we conclude 

{\hi, x \ + \gi, x \) \wi - 6iVi\ dxdt 
{h 2 i, x Vi + 9i, x Vi + 2i>i, x »fc) dxdt + 0(1) ■ / J {\h itX Vj\ + \gi, x Vj\) dxdt 




It J\g i /h i \>45 1 /5 

<-i; 



= <D(l)-5 2 . 

CASE 2: \gi/hi\ < 4<5i/5, \wi/vi\ > 5\. In this case we have 

4(511/1,1 K| 4 
5 01 5 

Using (11.14) we can write 

\hi, x \ \wi - 9iVi\ < \h i:X Wi\ 



g i -Cx i -X*)h i + O(l)-d J2{\h j \ + \v j \) 



(11.36) 



(11.37) 



Wi 



(11.38) 



< 5x1^1 + 0(1) -do Y (N + 
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By (11.28), the integral of the last terms on the right hand side of (11.38) is 0(1) • 5$. Concerning 
the first term, using (11.36) and then (6.16) and (11.14), we can write 

\\hiWi\ < \hiWi - giVi\ 
5 

< \Kvi, x ~ K, x Vi\ + \hi\ \ wi - v ijX \ + \vi\ \gi - h ijX \ 

= \hiV ijX - h ijX Vi\ + 0(1) • 5 \hi\(\wi\ + ^2 Kl) + ■ &o\wi\{\hi\ + ^ (\hj\ + 

Hence, for 5 small, one has 

\hiWi\ < 6\hiVi tX - h ijX Vi\ + 0(1) • ^2 {\hiVj\ + \wihj\ + \wiVj\) . 

By (11.28) and (11.30) we conclude 

/ / \hi tX \ \w{ — OiVi\ dxdt < / / \hi iX Wi\ dxdt 

Jt J\g t /h t \<48 1 /5 ' ' h J\gi/hi\<4:6i/5, \w i /v i \>8 1 

= O(l)-5 2 . 

Recalling that Si < 1, the last remaining term can now be bounded as 

\9i,x\ \wi - 0iVi\ < \g i:X Wi\ 

< \9i,xWi - giW i:X \ + \giW i>x \ 

< \9i, x Wi - giW ijX \ + —\hiW itX \ 

< \9i, x Wi - giW ijX \ + \hiW ijX - h itX Wi\ + \h ijX Wi\ . 

By (11.30) and (11.39) we conclude 

/ f \gi, x \\wi-0iVi\dxdt = O(l)-Sfi. (11.40) 

Jt J\gi/hi\<4Si/5 

This completes the proof of Lemma 11.6. □ 
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12 - Propagation speed 

Consider two solutions u, v of the same viscous system (1.10), whose initial data coincide 
outside a bounded interval [a, b]. Since the system is parabolic, at a given time t > one may well 
have u(t,x) ^ v(t,x) for all x € M. Yet, we want to show that the bulk of the difference \u — v\ 
remains confined within a bounded interval [a — (3t, b + (3t\. This result will be useful in the final 
section of the paper, because it implies the finite propagation speed of vanishing viscosity limits. 

Lemma 12.1 For some constants a, (3 > the following holds. Let u,v be solutions of (1.10) with 
small total variation, whose initial data satisfy 

u(0,x) =v(0,x) x£[a,b]. (12.1) 
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Then for all x 6 M, t > one has 

\u(t,x) -v(t,x)\dx < ||«(0)-«;(0)|| LCO ■mm{ae /3t -( x - b \ ae^ x ~ a ^ . 

On the other hand, assuming that 

u(0, x) = v(0, x) x £ [a, b] , 

one has 

\u(t,x)-v(t,x)\dx < ||u(0)-™(0)|| Loo • [ae^-^-^ +ae l3t+{x - b) ^. 

Proof. 1. As a first step, we consider a solution z of the linearized system 
z t + [^(u)^]^ + [DA(u) ■ z]u x - [DA(u) ■ u x ]z = z xx 
with initial data satisfying 

J \z(0, x)\ < 1 if x < 0, 
\ z(0,x)=0 if x>0. 

We will show that z(t,x) becomes exponentially small on a domain of the form {x > (3t}. More 
precisely, let B(t) be a continuous increasing function such that 



(12.2) 

(12.3) 
(12.4) 

(12.5) 



(t) dt , 



B(0) = 1. 



One can show that such a function exists, satisfying the additional inequality B(t) < 2e ct , for 
some constant C large enough and for all i > 0. We claim that 

\z(t,x)\ < E(t,x) = B(t) exp |4||L>yl|| LO o J \\u x (s)\\ Loo ds + 1 - x j (12.6) 

for all x 6 iR and i > 0. Indeed, any solution of (12.5) admits the integral representation 

z(t) = G(t) * z(0) - [ G x (t-s)* [A(u)z] (s)ds 
Jo 

+ / G(t-s)* (u x » A(u))z(s) - {h* A(u))u x (s) ds , 
Jo L - 1 

in terms of convolutions with the standard heat kernel G(t,x) = e~ x l u /2\firi. Therefore 
\z(t,x)\< J G(t,x-y)\z(0,y)\dy + \\A\\ L ~, J J G x (t - s,x - y) \z(s,y)\dyds 
+ 2\\DA\\ L ™ J^J\\u x {s)\\ Loo G(t -s,x- y)\z(s,y)\dyds. 
For every t > the following estimates hold (see Appendix D for details): 

-(x-y) 2 /4t 



(12.7) 



J G(t,x-y)\z(0,y)\dy< 



2v^ 



e y dy = e , 



(12.8) 
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\A\\ L ~ I l\G x (t-s,x-y)\E( S ,y)dyds<lE(t,x)-^e t - x , (12.9) 



'o 



2\\DA\\ L ~J^ \\u x (s)\\ LOO G(t - s,x - y)E(s,y)dy ) j ds < ±E(t,x) - \e l ~ x . (12.10) 

The bounds (12.7)-(12.10) show that, if (12.6) is satisfied for all t 6 [0, r[, then at time t = r one 
always has a strict inequality: \z(t, x)\ < E(t,x). A simple argument now yields the validity of 
(12.6) for alU > and x € M. 

2. Recalling (10.10) we have 

I, || (2k5 2kS \ 

|M*)|| L .<max|-^. -^j. 

From the definition of E at (12.6), for some constants a, (3 > we now obtain 

\z(t,x) \ <E(t,x) < 2e c *exp |4||DA|| L oo ■ 2k<5 (2v^ + t/Vi ) } e*" x < ae^~ x . (12.11) 

3. More generally, let now z be a solution of (12.5) whose initial data satisfies 

{|z(0, x) | < p if x < b , 
z(0,x)=0 if x>6. 

By the linearity of the equations (11.1) and translation invariance, a straightforward extension of 
the above arguments yields 

\z(t,x)\ <p-ae^ x - b K 

On the other hand, if 

{\z(0, x)\ < p if x > a , 
z(0, x) = if x < a , 

then 

|z(t,x)| < p ■ ae pt+{x - a \ 

4. Having established the corresponding bounds on first order tangent vectors, the estimates (12.2) 
and (12.4) can now be recovered by a simple homotopy argument. For each 6 6 [0, 1], let u be 
the solution of (1.10) with initial data 

u e (0) = 9u(0) + (l-8)v(0). 

Moreover, call z e the solution of the linearized Cauchy problem 

z° t + [DA(u e ) ■ z e ]u e x + A{u e )z e x = z e xx , 

z e (0,x) = u(0,x) - v(0,x). 
If (12.1) holds, then by the previous analysis all functions z e satisfy the two inequalities 

\z°(t,x)\ < \\u(0) - u(0)|| Lco • ae^-^, 
\z e (t,x)\<\\u(0)-v(0)\\^.ae^ x - a \ 
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Therefore 



\u(t, x) — v(t, x)\ < 

Jo 



du e (t,x) 



d6 



d6 



\t,x)\d6 



< ||m(0)-u(0)|| Loc •min{ae /3t - (a; - b) , ae^'^A . 



This proves (12.2). On the other hand, if (12.3) holds, we consider a third solution w of (1.10), 
with initial data 

-u(0, x) if x < b, 

v(0, x) if x > a. 

For every x 6 M and t > 0, the previous arguments now yield 



w(0, x) 



\u(t,x)-w(t,x)\ < \\u(0) -u»(0)|| Loo ■ae? t -( x - a \ 
\v(t,x) - w(t,x)\ < \\w{0) - u(0)|| LOO • ae /3t+(x - 6) . 

Combining these two inequalities we obtain (12.4). 



□ 



13 - The vanishing viscosity limit 



Up to now, all the analysis has been concerned with solutions of the parabolic system (1.10), 
with unit viscosity. Our results, however, can be immediatly applied to the Cauchy problem 



u £ t +A(u £ )u% = eu £ xx , 



u £ (0, x) = u(x) 



(13.1) 



for any e > 0. Indeed, as remarked in the Introduction, a function u £ is a solution of (13.1) if and 
only if 

u e {t,x) = u(t/e, x/e), (13.2) 
where u is the solution of the Cauchy problem 

u t + A(u)u x = u xx , 



u(0, x) = u(ex) . 



(13.3) 



Since the rescaling (13.2) does not change the total variation, from our earlier analysis we easily 
obtain the first part of Theorem 1. Namely, for every initial data u with sufficiently small total 
variation, the corresponding solution u s {t) = S £ u is well defined for all times t > 0. The bounds 
(1.15)-(1.17) follow from 



Tot.Var. {u £ (t)} = Tot.Var.{n(t/e)} < CTot.Var.{u}, 

1 „ 

£ 
t 



\u £ (t)-v £ (t)\\^=s\\u(t)- 
\\u E (t) - n £ (s)|| Ll < e\\u(t/e) - «(s/e)|| L i < eL' 



i < sL \\u(0) — "u(0)|| L1 = eh -\\u — 

/ t .s 



+ 



Moreover, if u(x) = v(x) for x € [a, b], then (12.4) implies 

' (3t — (x — a) 



\u £ {t,x) - v £ (t,x) \ < 



\U — V L°° - 



(I 

a exp I - 



^ + a exp ^ 



0t + (x- b) 



(13.4) 
(13.5) 

(13.6) 
(13.7) 
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We now consider the vanishing viscosity limit. Call U C Lj oc the set of all functions u : M>-> 
]R n with small total variation, satisfying (1.14). For each t > and every initial condition u £U, 
call S £ u = u £ (t, •) the corresponding solution of (13.1). Thanks to the uniform BV bounds (13.4), 
we can apply Helly's compactness theorem and obtain a sequence e v — > such that 



lim u £u (t, 



u(t, 



in 



(13.8) 



holds for some BV function u(t, •). By extracting further subsequences and then using a standard 
diagonalization procedure, we can assume that the limit in (13.8) exists for all rational times t and 
all solutions u £ with initial data in a countable dense set U* C U. Adopting a semigroup notation, 
we thus define 



S t u = lim Sf m u 



(13.9) 



for some particular subsequence e m — > 0. By the uniform continuity of the maps (t, u) t— > u £ (t, •) = 
S e u , stated in (13.5)-(13.6), the set of couples (t, u) for which the limit (13.9) exists must be closed 
in M+ x U. Therefore, this limit is well defined for all u € U and t > 0. 

Remark 13.1. The function u(t, •) = Stu is here defined as a limit in lij oc . Since it has bounded 
variation, we can remove any ambiguity concerning its pointwise values by choosing, say, a right 
continuous representative: 

u(t, x) = lim u(t, y). 

With this choice, the function u is certainly jointly measurable w.r.t. t,x (see [B5], p. 16). 

To complete the proof of Theorem 1, we need to show that the map S defined at (13.9) is 
a semigroup, satisfies the continuity properties (1.18) and does not depend on the choice of the 
subsequence {e m }- These results will be achieved in several steps. 

1. (Continuous dependence) Let S be the map defined by (13.9). Then 



sup 

r>0 



\\S t u - S t v\\ hl 
For every r > 0, the convergence in L^ oc implies 



(S t u)(x)-(S t v)( 



x 



dx . 



(S t u)(x)-{S t v)(x) 



dx 



lim 

m— >oo 



{S e r u)(x)-{SI™v)(x) 



dx < L\\u — v\ 



L 1 ' 



because of (13.5). This yields the Lipschitz continuous dependence w.r.t. the initial data: 

1 1 S f u — S f v 1 1 T , < L 1 1 u — v 1 1 T , . 



(13.10) 



The continuous dependence w.r.t. time is proved in a similar way. By (13.6), for every r > 
we have 



(S t v)(x) - (S s v)(x) 



dx 



lim 

m— »oo 



lim e m L' 

m— >oo 

L'\t-s\. 



{S £ r v)(x)-(S £ -v)(x) 
t 



dx 



t 



S 



+ 





S 
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Hence 

\\S t v - SsV^ < L'\t- s\. (13.11) 
Together, (13.10) and (13.11) yield (1.18). 

2. (Finite propagation speed) Consider any interval [a, b] and two initial data u, £>, with 
u(x) = v(x) for x 6 [a, b}. By (13.7), for every t > and x € ]a + (3t, b — (3t[ one has 



(S t u)(x) - (S t v)(x) 



< limsup 



(S^u)(x)-(S^v)(x) 



< lim — 

m-»oo 

= 0. 



• | a exp ^ 



/ft-fr-ah , ( Pt + {x-b) 
I + a exp I 



)} 



(13.12) 

In other words, the restriction of the function S t u 6 L i 1 oc to a given interval [a', b'] depends only 
on the values of the initial data u on the interval [a 1 — (3t, b' + f3t]. Using (13.12), we now prove a 
sharper version of the continuous dependence estimate (13.10): 



dx < L • 



b+pt 



a- fit 



u(x) — v(x)\ dx . 



[ (S t u)(x)-(S t v)( 

J a 

valid for every u, v and t > 0. Indeed, define the auxiliary function 

. T . M = j u(x) if x £ [a - /3t, b + 0t] , 

{ ' \ v(x) if x i [a - 0t, b + 0t] . 

Using the finite propagation speed, we now have 



(S t u)(x) - (S t v)(x) 



dx 



< L \\w — v\\ti 



(S t w)(x) - (S t v)(x) 

rb+fit 

L ■ / \u(x) — v{x) \ dx . 

Ja-fit 



dx 



(13.13) 



3. (Semigroup property) We now show that the map (t, u) i— ► S t u is a semigroup, i.e. 

S u = u, S s S t u = S s+t u. (13.14) 

Since every S £ is a semigroup, the first equality in (13.14) is a trivial consequence of the definition 
(13.9). To prove the second equality, we observe that 



S s+t u= lim S £ s m S E t m u, S s S t u= lim S £ s m S t u. (13.15) 

m — >oo m — >oo 

We can assume s > 0. Fix any r > and consider the function 

u m (x) = 
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(S t u)(x) if \x\<r + 2(3s, 

(Sl m u)(x) if \x\ >r + 2(3s. 



Since S Sm u — > S t u in L* oc , using (13.7) and (13.5) one obtains 



lim sup f (S £ s m S e t m u)(x)- (S £ s m S t u) (x) 

m— >oo J — r 



dx 



< lim 2r • sup 

m ^°° |x|<r 



(S e ™Sl™u)(x)-(S e s ™u m )(x) 



+ lim \\S E s m u m - S £ s m Sl m u\\ T1 



< lim 2r\\St m u-u m \\ Too • 2ae _/3s/£m + lim L-||n m -S' t £ 

m— too m— >oo 

= 0. 
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By (13.15), this proves the second identity in (13.14). 

4. (Tame Oscillation) We now exhibit a regularity property which is shared by all semigroup 
trajectories. This property, introduced in [BG], plays a key role in the proof of uniqueness. We 
begin by recalling the main definitions. Given a < b and r > 0, we denote by Tot.Var.{n(r) ; ]a, b[ } 
the total variation of n(r, •) over the open interval )a, b[ . Moreover, consider the triangle 

A l.b = {(t,x); t > t, a + [3t < x < b - [3t] . 

The oscillation of u over A^ b will be denoted by 

Osc.jn; A T atb } = sup {\u(t,x) - u(t',x')\; (t,x), (t',x') G A£ >b j. 

We claim that each function u(t, x) = (S t u) (x) satisfies the tame oscillation property: there exists 
a constant C such that, for every a < b and r > 0, one has 



Osc.jn; A T ab } < C ■ Tot.Var.{n(r) ; ]a,b[} 



(13.16) 



Indeed, let o, b, r be given, together with an initial data u. By the semigroup property, it is not 
restrictive to assume r = 0. Consider the auxiliary initial condition 



u(x) 
v(x) = ^ u(a+) 
u(b-) 



if a < x < b, 
if x < a, 
if x > 6, 



(13.17) 



and call v (i, x) = (S t v ) (x) the corresponding trajectory. Observe that 

lim v(t, x) = u(a+) 

x—> — oo 

for every t > 0. Using (1.15) and the finite propagation speed, we can thus write 

Osc.ju; A T a b } = Osc.{v; A T ab ] < 2 sup (Tot.Var.jS^}) 
< 2C • Tot.Var.{t;} = 2C • Tot.Var.{u(r) ; ]a, b[ } , 
proving (13.16) with C = 2C. 

5. (Conservation equations) Assume that the system (13.1) is in conservation form, i.e. A{u) = 
Df{u) for some flux function /. In this special case, we claim that every vanishing viscosity limit is 
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a weak solution of the system of conservation laws (1.1). Indeed, with the usual notations, if is a 
C 2 function with compact support contained in the half plane {x 6 M, t > 0}, one can repeatedly 
integrate by parts and obtain 

J J [u<f>t + f(u)(f) x ] dxdt = Yirn^ J J [u £m fa + f(u £m )4> x ] dxdt 

= - lim // \u £ t m (f) + f{u £m ) x (f\ dxdt = - lim // e m u £ x ^fadxdt 

m-»cxD J J m— >oo J J 



lim / / e m u £m (j) xx dxdt = 0. 

m-»oo 



An easy approximation argument shows that the identity (1.5) holds more generally, assuming only 
<f>eCl 

6. (Approximate jumps) From the uniform bound on the total variation and the Lipschitz 
continuity w.r.t. time, it follows that each function u(t,x) = (Stu)(x) is a BV function, jointly 
w.r.t. the two variables t, x. In particular, an application of Theorem 2.6 in [B5] yields the existence 
of a set of times M C M + of measure zero such that, for every (r, £) € M + x 1R with r ^ Af, the 
following holds. Calling 

u~ = lim u(t,x), u + = lim u(t,x), (13.18) 

there exists a finite speed A such that the function 

for every constant k > satisfies 

lim — / / \u(r + t, C + x) -U(t,x)\dxdt = 0, (13.20) 

1 f Kr 

- \u{t + r, £ + x) - U{r,x)\ dx = 0. (13.21) 

^ J — Kr 



lim 

r^0+ r 



In the case where u~ ^ u + , we say that (r, £) is a point of approximate jump for the function 
u. On the other hand, if u~ = u + (and hence A can be chosen arbitrarily), we say that u is 
approximately continuous at (r, £). The above result can thus be restated as follows: with the 
exception of a null set of "interaction times" , the solution u is either approximately continuous 
or has an approximate jump discontinuity at each point (t, £). 

7. (Shock conditions) Assume again that the system is in conservation form. Consider a 
semigroup trajectory u(t, ■) = S t u and a point (t, £) where n has an approximate jump. Since n is 
a weak solution, the states u~,u + and the speed A in (13.19) must satisfy the Rankine-Hugoniot 
equations 

A (u + — u~) = f(u + ) — f(u~). (13.22) 

For a proof, see Theorem 4.1 in [B5]. 

If u is a limit of vanishing viscosity approximations, the same is true of the solution U in 
(13.19). In particular (see [MP] or [D]), the Liu shock conditions must hold. More precisely, call 
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s i > Si(s) the parametrized shock curve through u and let Aj(s) be the speed of the corresponding 
shock. If u + = Si{s) for some s, then 

Ai(s') < Xi(s') for all a' £ [0, s] . (13.23) 

Under the additional assumption that each characterictic field is either linearly degenerate or 
genuinely nonlinear, it is well known that the Liu conditions imply the Lax shock conditions: 

K(u + ) < A < X l (u~). (13.24) 



8. (Uniqueness in a special case) Assume that the system is in conservation form and that 
each characterictic field is either linearly degenerate or genuinely nonlinear. By the previous steps, 
the semigroup trajectory u(t, •) = Stu provides a weak solution to the Cauchy problem (1.1)-(1.2) 
which satisfies the Tame Oscillation and the Lax shock conditions. By a well known uniqueness 
theorem in [BG], [B5], such a weak solution is unique and coincides with the limit of front tracking 
approximations. In particular, it does not depend on the choice of the subsequence {e m } : 

S t u = lim Sfu , 

i.e. the same limit actually holds over all real values of e. 

The above results already yield a proof of Theorem 1 in the special case where the system 
is in conservation form and satisfies the standard assumptions (H). To handle the general (non- 
conservative) case, we shall need to understand first the solution of the Riemann problem. 



14 - The non-conservative Riemann problem 

Aim of this section is to characterize the vanishing viscosity limit for solutions u e of (13.1), in 
the case of Riemann data 

u{x) = {\ if x<0 > (14.1) 
y ' \u+ if x>0. v ' 

More precisely, we will show that, as e — > 0+, the solutions u e converge to a self-similar limit 
u(t,x) = u>(x/i). We first describe a method for constructing this solution uj. 

As a first step, given a left state u~ and i G {1, . . . , n}, we seek a one-parameter curve of right 
states u + = ^i{s) such that the non-conservative Riemann problem 

u t + A(u)u; x = 0, lo(0,x) = \\ {i x < (14.2) 

[u^ if x > 

admits a vanishing viscosity solution consisting only of i- waves. In the case where the system is in 
conservation form and the i-th field is genuinely nonlinear, it is well known [Lx] that one should 
take 



Ms) = 




if s > 0, 
if s < 0. 



Here Ri and Si are the i-th rarefaction and shock curves through u , respectively. We now describe 
a method for constructing such curve in the general case. 
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Fix e, s > 0. Consider the family T C C°([0, s] ; iR™ x 1R x iR) of all continuous curves 
r i ^ 7(t) = (u(r), V;(r), 0-;(r)), r G [0, s] , 



with 



u(0) = u , |«(r) — « |<e, ki(''")|<e. 



(7 



e . 



In connection with a given curve 7 6 F, define the scalar flux function 

fi{l,r)= f Xi(u(q), Vi(q), <Ji(^)) r G [0, s] , (14.3) 

Jo 

where A, is the speed in (4.21). Moreover, consider the lower convex envelope 
conv/ l ( 7 ,r)=inf{^( 7 ,r') + (l-^/ l (7,r // ); 9 €[0,1], r',r"e[0,s}, r = 9r' + (1 - 9)r"] . 
We now define a continuous mapping % iS : V i-> V by setting 7^7 = 7 = (u, Vi, where 



u(t) = u + fi(u(q), Vi(q), ai(q)) dq , 
Jo 

Vi(r) = /i(7> T ) - conv/i(7,r) , 



(14.4) 



We recall that f, are the unit vectors that define the center manifold in (4.13). Because of the 
bounds 

\u(t) — U~ I < T < S , 

\&i{T) - \i(u-)\ =0(1)- sup A< («(?), fifc). -Ai(«-) = 0(1) -a, 

ce[o,s] 

|ui(r)| =<9(l)-a sup A;(u«), ^(s), a,( ? )) - A,(u") = 0(1) • s 2 , 
se[o,s] 

it is clear that for < s « e « 1 the transformation % tS maps V into itself. We now show that, 
in this same range of parameters, T iiS is a contraction with respect to the weighted norm 

|| (U, Vi,(7i)\\^ = + + e lki||L°= • 

Indeed, consider two curves 7,7' 6 T. For each r G [0, s] one has 

\\u - u'W^oo < / fi(u,Vi,(Ti) - fi(u' ^v'^a'i) dq 
Jo 

= 0(1) • s (\\ u ~ U '|| L - + W v i ~ v iW^ + ll^llL-Iki - Cr<|j L oo^ , 

||^-^||<2||/ l (7)-/ l (7 / )|| L o 

<2 / \i(u,Vi,(Ti) - Xiiu^v'^a-) dq 
Jo 

= O(l) • s (||u - u'|| LOO + - -u-|| L ~ + ||fi||L°=|ki - 0;||l<~) , 
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\°i( T ) - Gill < SU P 
rG[0,s] 



^-conv/i(7,r) - -^-conv fl(j,r) 



< sup 

tG[0,s] 



/i(7,r) 



//(7,r) 



<||A<-AJ|| L .=0(1)-(|| 
For some constant Co, the previous estimates imply 



11 ~ u IIl°° + W Vi 



|7-7 lit < Coeh-Yllt < 2'17-7'llt 



(14.5) 



provided that e is sufficiently small. Therefore, by the contraction mapping principle, the map % t 
admits a unique fixed point, i.e. a continuous curve 7 = (u,Vi,ai) such that 



u{t) = u + fi («(?), crj(?)) efc, 

■/ 

Vj(r) = /i(7,r) - conv/i(7,r) , 
^(r) = J^conv/i(7,r) . 



(14.6) 



Recalling the definition (14.3), from the continuity of u,Vi,<Ji it follows that the maps r i-> u(r), 
r 1 > fj(r) and r i-> fi(j,r) are continuously diferentiable. We now show that, taking u + = u(s) 
corresponding to the endpoint of this curve 7, the Riemann problem (14.2) admits a self-similar 
solution containing only i-waves. 

Lemma 14.1 In the previous setting, let 7 : r 1— > («(t), fi(r), o"i(r)) 6e i/ie /ixed point of the 
transformation T i>s . Define the right state u + = u(s). Then the unique vanishing viscosity solution 
of the Riemann problem (14-2) is the function 



u>(t, x) 




if x/t<ai(0), 
if x/t = ct(t) , 
if x/t > a(s) . 



[14.7) 



lim \\u)(t) — S £ oo 
£^0+ 11 w 1 



Proof. With the semigroup notations introduced in Theorem 1, we will show that, for every t > 0, 

J L i=0. (14.8) 
The proof will be given in several steps. 
1. Assume that we can construct a family v E of solutions to 

v t + A(v)v x = ev xx , (14.9) e 

with 



lim \\v £ (t) -u(t)\ 



L 1 







(14.10) 



for all t 6 [0, 1]. Then (14.8) follows. Indeed, by a simple rescaling we immediately have a family 
of solutions v £ such that (14.9) £ -(4.10) hold on any fixed interval [0, T\. For every t G [0, T], since 
by assumption v £ (t) = S E v E (0), using (1.16) we obtain 

£ hm \\uj(t) - S!uj(0)\\^ < £ hm \\u(t) - « e (t)|| Ll + \im + \\v £ (t) - 5 t e o;(0) || Ll 

<0 + L. £ hm ||^(0) -w(0)|| Ll = 0. 
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2. For notational convenience, call VVL the set of all vanishing viscosity limits, i.e. all functions 
v : [0, 1] x M »-> M n such that 

£ Hm||^(i)-^)|| Ll =0 i€[0,l] (14.11) 

for some family of solutions v £ of (14.9) e . By Step 1, it suffices to show that the function to at 
(14.7) lies in VVL. 

Let us make some preliminary considerations. Consider a piecewise smooth function v = v(t, x) 
which provides a classical solution to the quasilinear system 

tk + A(v)v x = * € [0, 1] , 

outside a finite number of straight lines, say x = Xj(t), j = 1,...,N. Assume that there exists 
5 > and constant states uj , such that 

( uj if Xj (t) — 5 < x < xj (t) , 

U ^ ' X ' ~ \ u+ if Xj(t) < x < Xj(t) + 5 , 

Moreover, assume that each pair of states n~,u^ can be connected by a viscous travelling wave 
having speed Xj. Finally, let v be constant on each of the two regions where x > r or x < — r , 
for some rg sufficiently large. Under all of the above hypotheses, it is then clear that v 6 VVL. 
Indeed, a family of viscous approximations v 6 can be constructed by a simplified version of the 
singular perturbation technique used in [GX] . 

As a second observation, notice that if we have a sequence of functions v m 6 VVL with 

lim ||M*)-^)|| L i =0 t€ [0,1], 

then also v G VVL. 

3. Consider first the (generic) case where the set of points in which /j is disjoint from its convex 
envelope is a finite union of open intervals (fig. 12), say 

N 

{r€[0,a]; / i (7,r)>conv/ i ( 7 ,r)} = lJ]a j ,6 j [. (14.12) 

i=i 

Our strategy is to prove that lo 6 VVL first in this special case. Later we shall deal with the 
general case, by an approximation argument. 




figure 12 figure 13 
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If (14.12) holds, we can make the two following observations. 

(i) For each j = 1, . . . , N, we claim that the left and right states u(aj), u(bj) are connected by a 
viscous travelling profile U such that 

U" = (A(U) - (Tij)U', U(-oo) = u(a 3 ), U(po) = u(bj). (14.13) 

Here <jjj is the constant speed 

(Jij = cn (r) = ^- conv (7, r) r G [aj , &,•] . 

To construct the function J7, consider the variable transformation ]aj, bj[^> M, say r i-> x(t), 
defined by 

j aj + bj \ _ dx(r) _ 1 

V 2 J ' dr ^(r)" 

Let r = t(x) be its inverse. Then the function U(x) = u[t(x)) is the required travelling wave 
profile. Indeed, U obviously takes the correct limits at ±00. Moreover, 

. du dr 

U = — — = ViTi, 

dr dx 
= Vi(v iyT ri + Vif iiT ) 

= Vi(\i - Oij)fi + vf(r iyU fi + ^^(Ai - <Tij)vi). 

Recalling the identity (4.22), we see that U also satisfies the differental equation in (14.13), thus 
proving our claim. 

(ii) On the intervals where /i(7,r) = conv/i(7,r) we have Vi(r) = 0. Hence, by the first equation 
in (14.6) and by (4.16), u T = fj = r^u) is an i-eigenvector of the matrix A{u). 

4. In general, even if the condition (14.12) is satisfied, we do not expect that the function u has 
the regularity specified in Step 2. However, we now show that it can be approximated in L 1 by 
functions u>§ satisfying all the required assumptions. To fix the ideas, let 

= 60 < a i < h < a 2 < b 2 < ■ ■ ■ < a N < b N < o>n+i = s . 

A piecewise smooth viscosity solution can be defined as follows (fig. 13). Fix a small 5 > 0. For 
each k = 0, . . . , N, consider a smooth non-decreasing map 

T fc : [kS, (k + 1)5] i-> [b k , a k+1 ] 

such that 

T k (x) = 

We then define the initial condition 



bk if x < k5 + 5/3, 

a k+1 if x>k5 + 25/3. 



lus(0,x) 



' u(0) if x<0, 

u(r fc (x)) if k5 < x < (k + 1)5 , 

k tt(s) if x>(N + l)5. 
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A corresponding solution of the Cauchy problem can then be constructed by the method of char- 
acteristics: 



co s (t,x) 



' u (0) = u~ if x < <Ti(0)t, 

-u(O) if x = y + (Ti{jk{y))t for some y G ] kS, (k + 1)8[ , 

^u(s) = u + if x > (N + l)S + (Ti(s)t. 



It is clear that the above function ujs satisfies all of the assumptions considered in Step 2. Hence 
u>s G VVL. Letting 8 — > we have ||tc>a(£) — w(i)|| Ll — > for every t. Therefore, by the last 
observation in Step 2 we conclude that also u> G VVL. 

5. To prove the Lemma in the general case, where the set in (14.12) may by the union of infinitely 
many open intervals, we use an approximation argument. For each 8 > 0, by slightly perturbing 
the values of A, we can construct a second matrix valued function A' with 

sup | A' (it) - A(u)\ < 8, (14.14) 

u 

such that the following properties hold. For some right state u + with \u + — u + \ < 8, the non- 
conservative Riemann problem 

u t + A\u)u x = 0, u (0,x) = \ u ~ {i f x< 0> (14.15) 

[ vr if x > 0. 



admits a self-similar solution u' which is limit of vanishing viscosity approximations and satisfies 

-3/3 
-3/3 



I \u'(t,x)-u(t,x)\dx<8 t€[0,l]. (14.16) 

J-30 



Clearly, the fact that lo' is a limit of vanishing viscosity approximations can be achieved by choosing 
A' so that a corresponding transformation T s , will admit as fixed point some curve 7' : r t— > 
(u'(t) ,?4(t) , cr^(r)) for which u'(0) = u~ , u'(s') « u + and with /i(7',r) differing from its convex 
envelope on a finite number of open intervals. 

Call u £ the solution of the viscous Riemann problem (13.1) with initial data (14.1). Using 
(13.7) with v E = u + , a = 0, b = 00, for all t G [0, 1] we obtain 



poo poo 

limy \co £ (t,x) -u(t,x)\dx < limy \u~ - u+\ • ae (/Jt_!B)/e = . (14.17) 

Similarly, 

,-/3 

limy |w £ (t,x)-w(t,x)|dx = 0. (14.18) 
To establish the convergence also on the interval [—(3, ft], call v E the solution of the Cauchy problem 

if x < , 

< + A»§ = et&. , ^(0,x) = <(n+ if 0<x<3ft, 



if x > 3/3 . 



Clearly, 

f/3 



limy \v s (t,x) -u/(t,x)\ dx = 0. (14.19) 



-0 
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because u' is a vanishing viscosity limit and because of the finite propagation speed. Using the 
triangle inequality we can write 



/P eP 
\ui s (t, x)—uj(t, x)\ dx < limsup / \u) £ (t,x) — v £ (t,x)\dx 



+ lim / \v £ (t, x) — co'(t, x)\ dx + / \u'(t,x) — u(t,x)\dx . 



Since t io £ {t) = S £ lo(Q) is a trajectory of the Lipschitz semigroup S £ , recalling (14.14) we have 
the estimate 



\ v ^ s + h)-S £ h v £ (s)\\^ 



\\v £ (t)-u, £ (t)\\ L1 dx<L\\v £ (0)-u; £ (0)\\ L1+ L- ^ Mim + ^ - " ^ j ds 

<3pL\u + -u + \+ L- |y A(v £ (s,x)) -A'(v £ (s,x)) \v £ x (s , x)\ dx^ ds 

< 3PL5 + L8C -Tot. Var.{v £ (0)} 

< C"5 , 

(14.21) 

for some constant C". Estimating the right hand side of (14.20) by means of (14.21), (14.19) and 
(14.16), we obtain 

limsup / \oj £ (t,x) - u(t,x)\dx < C"5 + + 5. 

e^O J-/3 

Since 5 > can be arbitrarily small, together with (14.17)-(14.18) this yields 

lim \\u £ {t) -<j(t)|| Ll =0 for all t€ [0,1], 

completing the proof. D 

Remark 14.2. The transformation % )S defined at (14.4) depends on the vectors fj, and hence 
on center manifold (which is not unique). However, the curve 7 that we obtain as fixed point of 
% iS involves only a concatenation of bounded travelling profiles or stationary solutions. These are 
bounded solutions of (4.2), and will certainly be included in every center manifold. For this reason, 
the curve 7 (and hence the solution of the Riemann problem) is independent of our choice of the 
center manifold. 

For negative values of the parameter s, a right state u + = ^i(s) can be constructed exactly 
in the same way as before, except that one now takes the upper concave envelope of fi : 

cone fi(j, t) = sup {^(7, t') + (l-0)/i(7,T"); 6 6 [0, 1] , t',t"g[0,s], r = Or' + (1 - 9)r"] , 

instead of the lower concave envelope. 

Our next step is to study the regularity of the curve of right states u + = ^i(s). 

Lemma 14.3. Given a left state u~ and i G {l,...,n}, the curve of right states s ^i(s) is 
Lipschitz continuous and satisfies 

lim = r .( u -y (1422) 

s^o ds 
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Proof. We assume s > 0, the other case being entirely similar. For sake of clarity, let us 
introduce some notations. For fixed i and s > 0, let 7*' s = (u l,s , vf , cr|) be the fixed point of 
the transformation % )S in (14.4). Then by definition 

For < s' < s, let 7' = (u', v-, c^) be the restriction of 7*' s to the subinterval [0, s']. Since % iS i is 
a strict contraction, the distance of 7' from the fixed point of % iS i is estimated as 

||Y - 7 M ' || t = oil) ■ || 7 ' - W|| t = 0(1) ■(s-s')s. 

In particular, 

|^' S V) - V) I = 0(1) ■(s-s')s. 

Observing that 

U <.»( S ) _ U i.-( S ') = f f^'fr), *i M) * = (a - *') ■ n(u-) + 0(1) ■ (s - s')s , 

J s 

we conclude 

|*i(a') - ^(a') - (a - s')^(«")| = ■(*-«')*■ ( 14 - 2 3) 

By (14.23), the map s 1— > ^(s) is Lipschitz continuous, hence differentiable almost everywhere, by 
Rademacher's theorem. The limit in (14.22) is again a consequence of (14.23). □ 

Thanks to the previous analysis, the solution of the general Riemann problem (14.2) can now 
be constructed following a standard procedure. Given a left state u~ , call s 1— > *&i(s)(u~) the curve 
of right states that can be connected to u" by i- waves. Consider the composite mapping 

* : (Sl, ...,*„)■-► ^n(Sn) O ■ ■ ■ O l(Sl)(lT) . 

By Lemma 14.3 and a version of the implicit function theorem valid for Lipschitz continuous 
maps (see [CI], p. 253), is a one-to-one mapping from a neighborhood of the origin in M n onto a 
neighborhood of u~ . Hence, for all u + sufficiently close to u~ , one can find unique values s±, . . . , s n 
such that VP(si, • • • ,s n ) = u + . In turn, this yields intermediate states uq = u~,Ui, ... ,u n = u + 
such that each Riemann problem with data Uj_i,Uj admits a vanishing viscosity solution uoi = 
u>i(t, x) consisting only of i-waves. By strict hyperbolicity, we can now choose intermediate speeds 

-00 = x' < a; < x' 2 < ■ ■ ■ < a;_! <x' n = co 

such that all i-waves in the solution u>i have speeds contained inside the interval [X' i _ 1 , A-]. The 
general solution of the general Riemann problem (14.2) is then given by 

u(t, x) = LOi(t, x) for A-_! < I < A- . (14.24) 

Because of Lemma 14.2, it is clear that the function u is the unique limit of viscous approximations: 

lim \\u(t) - 5 t £ u;(0)|| L1 = for every t > 0. (14.25) 
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15 - Viscosity solutions and uniqueness of the semigroup 



In [B3], one of the authors introduced a definition of viscosity solution for a system of con- 
servation laws, based on local integral estimates. Assuming the existence of a Lipschitz semigroup 
of entropy weak solutions, it was proved that such a semigroup is necessarily unique and every 
viscosity solution coincides with a semigroup trajectory. We shall follow here exactly the same 
approach, in order to prove the uniqueness of the Lipschitz semigroup constructed in (13.9) as 
limit of vanishing viscosity approximations. 

Toward the definition of a viscosity solution for the general hyperbolic system 

u t + A(u)u x = 0, (15.1) 

we first introduce some notations. Given a function u = u(t,x) and a point (t,£), we denote by 
U( u . t q the solution of the Riemann problem (14.1) with initial data 

u~ = lim u(t,x), u + = lim u(t,x). (15.2) 

Of course, we refer here to the vanishing viscosity solution constructed in Section 14. In addition, 
we define U^ U . T ^ as the solution of a linear hyperbolic Cauchy problem with constant coefficients: 

w t + Aw x = 0, w(0,x) = u(t,x). (15.3) 

Here A = A(u(t, £)). Observe that (15.3) is obtained from the quasilinear system (15.1) by 
"freezing" the coefficients of the matrix A(u) at the point (r, £) and choosing u{j) as initial data. 

As in [B3], the notion of viscosity solution is now defined by locally comparing a function u 
with the self-similar solution of a Riemann problem and with the solution of a linear hyperbolic 
system with constant coefficients. 

Definition 15.1. A function u = u(t, x) is a viscosity solution of the system (15.1) if 1 1— > u(t, •) 
is continuous as a map with values into L* oc , and moreover the following integral estimates hold. 

(i) At every point (t,£), for every (3' > one has 

rZ+P'h 



lim 

/i->o+ h 



T \ u(r + h, x)-UL Ti Jh, x-0 dx = 0. (15.4) 



(ii) There exist constants C, (5 > such that, or every r > and a < £ < b, one has 
-b-ph 

la+Ph 



1 r b -> 3h z n2 

lim sup - / u(r + h, x)-UL T£) (h,x) dx < C ■ (Tot.Var.{u(r); ]a, b[ } J . (15.5) 
h^0+ h J a+ Bh v ' 



The main result of this section shows that the above viscosity solutions coincide precisely with 
the limits of vanishing viscosity approximations. 

Lemma 15.2. Let S : V x [0, oo[^ V be a semigroup of vanishing viscosity solutions, constructed 
as limit of a sequence S Em as in (13.9) and defined on a domain V C \j} oc of functions with small 
total variation. A map u : [0, T] V satisfies 

u(t) = S t u(0) for all t G [0, T] (15.6) 
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if and only if u is a viscosity solution of (15.1). 



Proof. Necessity: Assume that (15.6) holds. By (13.11), the map t i-> u(t) is continuous. Let 
any f3' > be given and let L, (5 be the constants in (13.13). Then, for any (t, £), an application 
of (13.13) yields 



< L - i 

\ >ti-{f}+f3')h 

<L{(5 + (3')h 



\ u(r + h, x)-m AK x-i) 
u(t,x) - u(t,£ 



dx 



dx + 



I 



u(t,x) - u(t,£+) 



sup 

€-(/?+/?')h<*<€ 



\u(t,x) - u(t,£-)\ + sup 

i<x<z+{0+0')h 



dx j 

(t,x)-u(t,£+)\ 



Hence (15.4) is clear. 

To prove the second estimate, fix r and a < £ < b. Define the function 

u(t, a+) if x < a , 

v(x) = < u(t, x) if a < x < b, 

k u(t, b—) if x > b. 

Call v E ,w E respectively the solutions of the viscous systems 

v £ t +A(v £ )v x = ev : 



J XX 1 



(15.7) 



with the same initial data v E (0,x) = w £ (0,x) = v(x). 

Recalling that S E is a semigroup with Lipschitz constant L, as in [B3], [B5] we can use the 
error formula 



\w £ (h)- ^(^11^ = 11^(^-5^11^ 



< L ■ 



< L 



lim inf 



U £ (t + r)-5> £ (s)|| T1 



dt 



A - A(w £ (t,x)) \w x (t,x)\dxdt 
A(w%0,0) -A(w e (t,x))^ -sup|K(t)|| Ll 
< Ch (Tot.Var.{w}) 2 , 



< L h ( sup 

t ,x 



for some constant C. Letting e — > and using the estimate (13.13) on the finite speed of propaga- 
tion, we obtain 



6-/3/i 



a+f3h 



u(r + h, x)-U\ u . rA) {h, 



dx < — lim 

h e^O 



b-0h 



a+f)h 



\v £ (h, x) — w £ (h, x)\ dx 



< C (Tot.Var.{7}}) 2 = C (Tot.Var.{u ; ]a,6[})' 
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This proves (15.5), with /3 the constant in (13.13). 



Sufficiency: Let u = u(t,x) be a viscosity solution of (15.1). By assumption, the map t i-> u(t) 
is continuous with values in a domain T> C L* oc of functions with small total variation. From 
(15.5) and this uniform bound on the total variation it follows that this map is actually Lipschitz 
continuous: 

||u(t)-u(s)|| Ll <L"\t-s\, (15.8) 

for some constant L" and all s,t G [0, T]. Let L be the Lipschitz constant of the semigroup S, as 
in (13.13). Given any interval [a, b], thanks to (15.8) one has the error estimate 



b-t/3 



a+t/3 



u(t,x)- (S t u(0))(x) 



< L ■ 



lim inf — 

h^0+ h 



dx 

b-(r+h)l3 
a + (r + h)f3 



(15.9) 



u 



(t + h, x) — (5ftit(r)) (x) \dx > dr . 



To prove the identity (15.6) it thus suffices to show that the integrand on the right hand side of 
(15.9) vanishes for all r £ [0,T]. 

Fix any time r 6 [0, T] and let e > be given. Since the total variation of u(r, •) is finite, we 
can choose finitely many points 

a + t(3 = x G < x\ < ■ ■ ■ < x N = b — t/3 

such that, for every j = 1, . . . , TV, 

Tot.Var. {u(t, •) ; ]xj-i, Xj[} < e. 

By the necessity part of the theorem, which has been already proved, the function w(t, •) = 
St- T u{r) is itself a viscosity solution and hence it also satisfies the estimates (15.4)-(15.5). We 
now consider the midpoints yj = (x?-i + Xj)/2. Using the estimate (15.4) at each of the points 
£ = Xj and the estimate (15.5) with £ = yj on each of the intervals Xj[, taking (3 > 

sufficiently large we now compute 

r b-{r+h)P 



lim sup 



1 f~ 
h Ja+( 

< > lim sup — / 



«(r + h, x) - (S h u(r)) (x) 

Xj +hf3 
hp 



dx 



(ju(T + h,x)-U( u . TXj) (T + h, X) 



+ 



U Ur,x^ + h, x)-(S h u(r))(x) 



dx 



N 



1 r*i-hP 
+ > lim sup — / 

" ' h -'xj-i 



, =1 ^-0+ 



u(t + h,x)- Ul u . Tjyj) (h, x-Xj) 



+ 



U, 



N 

< + J2 C (Tot.Var. {u(t); ]xj-i, Xj[}^j 



{u-,T,yj 

2 



;) (/l, X-Xj) - (S h u(T))(x) 



dx 



3 = 1 

< Ce- Tot.Var. {u(t); ]a + T/3, b-T/3[}. 

Since e > was arbitrary, the integrand on the right hand side of (15.9) must vanish at time r. 
This completes the proof of the lemma. □ 
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Remark 15.3. From the proof of the sufficiency part, it is clear that the identity (15.6) still holds 
if we require that the integral estimates (15.4) hold only for r outside a set of times M C [0, T] of 
measure zero. By a well known result in the theory of BV functions [EG], any BV function of two 
variables u = u(t, x) is either approximately continuous or has an approximate jump discontinuity 
at every point (r, £), with r outside a set M having zero measure. To decide whether a function u 
is a viscosity solution, it thus suffices to check (15.4) only at points of approximate jump, where 
the Riemann problem is solved in terms of a single shock. 

Using Lemma 15.2, we now obtain at one stroke the uniqueness of viscosity solutions and of 
vanishing viscosity limits: 

Completion of the proof of Theorem 1. What remains to be proved is that the whole family 
of viscous approximations converges to a unique limit, i.e. 

lim S e t u = S t u, (15.10) 

where the limit holds over all real values of e and not only along a particular sequence {e m }- If 
(15.10) fails, we can find v, r and two different sequences e m ,e' m — > such that 

lim S e T m v^ lim Sr' m v . (15.11) 

m— too m-»cxD 

By extracting further subsequences, we can assume that the limits 

lim S £ t m u = S t u, lim St' m u = S' t u, (15.12) 

m— >oo m— >oo 

exist in L* oc , for all t > and u <EU. By the analysis in Section 13, both S and S' are semigroups of 
vanishing viscosity solutions. In particular, the necessity part of Lemma 14.2 implies that the map 
t 1 > v(t) = S t v is a viscosity solution of (15.1), while the sufficiency part implies v(t) = S' t v(0) for 
all t > 0. But this is in contradiction with (15.11), hence the unique limit (15.10) is well defined. 

□ 

Remark 15.4. The above uniqueness result is obtained within the family of vanishing viscosity 
limits of the form (1.13) e , with unit viscosity matrix. In the more general case (1.21) £ , if the 
system is not in conservation form, we expect that the limit of solutions as e — ► will depend 
on the form of the viscosity matrices B(u). Indeed, by choosing different matrices B(u), one will 
likely alter the vanishing viscosity solutions of the Riemann problems (14.2). In turn, this affects 
the definition of viscosity solution at (15.4). 



16 - Dependence on parameters and large time asymptotics 

We wish to derive here a simple estimate on how the viscosity solution changes, depending on 
hyperbolic matrices A(u). 

Corollary 16.1. Assume that the two hyperbolic systems 

u t + A(u)u x = 0, 
u t + A{u)u x = , 
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both satisfy the hypotheses of Theorem 1. Call S,S the corresponding semigroups of viscosity 
solutions. Then, for every initial data u with small total variation, one has the estimate 

\\S t u - S t «|| Ll = 0(1) • t(sup \A(u) - A(u)\) • Tot.Var.{u} . (16.1) 

Proof. Call S £ , S e the semigroups of solutions to the corresponding viscous problems 

u t + A(u)u x = eu xx , ut + A(u)u x = eu xx . 

Let L be the Lipschitz constant in (1.16) and call w £ (t) = S £ u. For every t > we have the error 
estimate 

||^-5^|| L1 = ||u; £ (t)-^|| L1 

r* [ \\w e {s + h)- s e h w e {s)\lA 

< L ■ / <^ liminf ^ — \J»LL \ ds 

Jo { h ^°+ h J 

< L ■ J J y4(w e (s,x)) — A(w £ (s,x)) \w e x (s,x)\dxds 

< L (sup \A(u) - J |K(s)|| L1 ds 
= O(l) ■ t (sup \A(u) - A(u)\) Tot.Var.{n}. 

Next, we show that some semigroup trajectories are asymptotically self-similar. 



□ 



Corollary 16.2. Under the same assumption of Theorem 1, consider an initial data u with small 
total variation, such that 



/0 /-oo 
\u(x) — u~ | dx + / \u(x) — u + \ dx < oo , 
-oo JO 



(16.2) 



for some states u ,u + . Call uj(t,x) = u)(x/t) the self-similar solution of the corresponding Rie- 
mann problem (14-2). Then the solution of the viscous Cauchy problem 

u t + A(u)u x = u xx , u(0,x)=u(x) (16.3) 

satisfies 

lim / \u(t, ry) — u>(y)\ dy = . (16.4) 

T— tOO J 

Proof. The assumption on u implies that the limit (14.25) holds. For fixed r, call e = 1/r and 
consider the function v E (t,x) = u(tx, rt). Clearly, v £ satisfies the equation 

vf + A{v £ )v x =ev £ xx , v £ {0,x) =u(x/e). 

Therefore, 

J \u(t, ry) - u(y)\dy = J \v £ (l, x) - w(l, x) \ dx 



< \\Sfv £ (0) - SM0)|| Ll + \\S £ M0) ~ w(l)|| Ll . 
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(16.5) 



Observing that 



\Slv%0) - SM0)|| Ll < L ■ \\Sfv%0) - S^(0)|| L1 

= Le (^j \u(x) — u~ | dx + J \u(x) — u + 1 dx^j , 



and using (14.25), from (16.5) we obtain (16.4). 



□ 



Appendix A 

We derive here the explicit form of the evolution equations (6.1), for the variables Vi and Wi 
defined by the decomposition 

u x = ^2vifi(u,Vi,\* - 6(wi/vi)), u t = ^{wi - \*Vi)fi(u,Vi,\* - 9{wi/vi)). (Al) 



By checking one by one all source terms, we then provide an alternative proof of Lemma 6.1. The 
computations are lengthy but straightforward: one has to rewrite the evolution equations for u x 
and u t : 

f (u x ) t + (A(u)u x ) - (u x ) xx = , 

{ K)t + (A(u)u t ) x - (u t ) xx = (u x • A(u))u t - (u t • A(u))u x , 

in terms of Vi,Wi. For convenience, we set 6i = 9(wi/vi). The fundamental relation (4.23) can be 
written as 

Vih,ufi - A(u)?i = -Ajfj + (-A; + A* - 6i)vir i:V . (A3) 
Differentiating (A.l) w.r.t. x and using (A. 3) we obtain 

A(U)U X = ^2 Vi, x fi + ^2 V ifi,x ~ ^2 A ( U ) V i f i 

III 

= ^2 Vi > xfi + ^2 Vi [ v iri,ufi - A(u)f 

i i 

+ ^2vi [v i;X r i;V - 6>- (( 

i 

= ^2(vi, x - Kvi)fi + ^2(-Xi + A* - Vi)v-ri, v 

i i 

+ ^2 v i[ v i,xfi,v ~ Oi((ViWi tX - WiV^ x )/vf)fi,^\ + ^2 v iVjfi, u rj 
i ijij 

= ^2 ( Vi > x ~ ^ iVi ) fi + Vi ^ v + e 'i{ W il V i) f i,o\ + ( Wi < x ~ ~ XiWi ) ~ 

i i 

+ ^2 v i{ X i - e i)?i,v +^2v i v j r i>u r j , 

i ijij 

(A4) 



i,x)/vf)fi,o\ +J2'' 
.2- 
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u tx - A(u)u t = ^2(w ijX - X*v i;X )fi + ^(wi - \*Vi)f ijX - ^(wi - \*Vi)A{u)fi 

i i i 

= ^2{wi, x - X*v i:X )fi ^(wi - X*Vi) [vifi,uh - A(u)fi] 



^2(wi,x - \wi)fi - ^A*(w iiX - XiVi)fi + ^2(wi - X*Vi)(- Xi + X* - 9i)vif i)V 

i i i 

i i^j 

^2 {vi, x - XiVi) Wir i>v + 0' i {w i /v i ) 2 r iia \ + ^ (w i>x - A* to*) [f; - (O-Wi/vi)^ 

i i 

+ V i W i (K ~ &i)fi,v + ^2 
i ijij 

~ ^2K\ ( V i,x ~ XiV?J \ri +V i f i> v + {O'iWi/Vijfi^ + (w itX - XiW^j 



Differentiating (A.l) w.r.t. t one obtains 



(A5) 



U xt = ^ V i>tfi + ^2 ViTi t 
i i 

= 'Y^VijVi + ^ V i [ V i,tn,v ~ 0'i{(Wi,tVi ~ WiV ijt )/vf)f io ^ +^2vi(wj - XjVj) f^ u fj ^ ^ 
i i i,j 

= y~] Vi,t \rj + Vjfi,v + {0' i 'Wi/vi)f i ^ + ^2w itt -9'ifi^ + ^] Vj(wj - X*Vj)r ijU rj , 



U « = ^2( W i,t ~ K V i,t)h + ^2{wi~ X*Vi)f i>t 
i i 

i i 
+ J2{Wz- K V i) ( W 3 ~ X j V j)fi,u?j 



i,3 



^2v itt \wifi,v + 9i{wi/Vi) 2 fi,a\ +^2w iti 

i i 

- ^2K\ V i,t h +Vifi,v + {0'iWi/Vi)f iyC 



fi - O'^Wi/vi)?^ + ^2wi(wj - X*Vj)r ijU fj 

i,3 

- Wi,t9'ih,a + ^2 Vi ( w i ~ X j v 3)ri,urj 



(A.7) 
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Differentiating again u xx — A{u)u x and u tx — A(u)u t w.r.t. x, from (A.4) and (A. 5) one finds 



Utx = [u a 



(A(u)u x ) : 



i i 

+ J^K* - kvi) v 3^,urj + 2v iiX f iiV + ( - 6 iiX + (O-Wi/vi)^^ + J2 

* 3 3 

+ ViV i>x f i>vv + ( - Vi6 i:X + 6iVi, x Wi/vi)r ijV<T + ^ v j Q' i {w i lv i )f i ^ u f j - 0i, x 0' i (w i /v i )f i!C 

3 

i L j 

i i j 

i^j i^j 1 k 



^j,xfi,ufj,v Vi jX Ti tUv rj 6 ^j jX Ti jU T 'j j(T 0i,x?'i,uo-f'j 



{u t ) xx - (A(u)u t ) x 

= ^ (Vi, xx - {\Vi) x ) Wif i>v + O'^Wi/ViffiA + ^ ( w i,xx - (hwi)x) h - 0' i (w i /v i )r i 



(A.8) 



+ ^2(vi, x - Kvi) w ijX r i;V + ^2wiVj f ijVU rj + WiV i;X f i>vv + ( ~ mOi,x + 0'i(w i /v i ) 2 v i ^ f 

i 1 j 

+ (Oi(wi/vi) 2 ) x fi )f7 + ^2v j 9' i (w i /v i ) 2 f iy(7U f j - 9i(wi/vi) 2 e iiX fi jtT 

j 

+ ^2( w i,x - AjtOj) ^Vjfjfi^ + Vi^ri^ - (0 i:X + (0' i w i /v i ) x )r i!a - 1 ^2v j 0' i {wi/vi)f i ^ u f j 



- Vi^OKtVi/Vi)?^ + 0i, x O'i{Wi/Vi)f i>c 



+ J2(wiVi(X* -9i)) x fi,v 



~l~ ^ ] WjVj (A^ ^ ^ VjVijVvT j Vi,xT'i,vv 0i,x^"i,va ~l~ ^ ^ WjVj^ x fi tU fj 



i¥=3 



^ ^ WjVj ^ ^ Vk (j~i,uu ® ^fc) "I - ^i,ufj,ufk) Vj jX Ti u fj v + Vi jX f{ vu rj 9j,x^i,u^j,(7 0i,x^"i,z 



k 



^2K\( V i,x ~ \Vi) fi +Vifi,v +0'i{Wi/Vi)f iy(T +{w iyX -\ i W i ) -9'ifi^ 

i 

+ V i (A* - + XI !''■■« '\i 

3+i J x 



(A9) 
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Substituting the expressions (A.6)-(A.9) inside (A. 2) and observing that 

(u x • A(u))u t - (tit • A(u))u x = y~](wj - \*Vi)vj (fj • A(u))fi - (f, • A{u))fj 

we finally obtain an implicit system of 2n scalar equations, describing the evolution of the compo- 
nents Vi,Wi'. 



i 



fi + Vifi,v + 6'i{wi/vi)fi,cT \ + ^2( w i,t + fti w i) x - 
Vi(v i:X - XiVi) - Vi(wi - X*Vi) 

+ ^ h,ufj (Vi, x - XiVi)vj + (viVj) x - Vi(wj - X*Vj) 

+ ^h,v 2v i:X (v i:X - XiVi) + (vt(X* - 9i)) x 

i 



IV 



i . xx 



i 



'i.vii'i 



Vi(v i:X - XiVi) + vf(X* - 9i) 



iVj{vi, x ~ hvi) + v j v i (X* - Oi 

+ ^2 ^> vv V i V i,x{ V i,x ~ AjVi) + V i>x Vi(X* - 6i) 
i 

+ ^2n,va {Vi,x ~ \Vi) (-Vi0 itX + O'iVi^Wi/v^j - (w itX - XiWi)v ijX 9'i - vf (A* - Oi)di 

i 

+ ^ ri,aufi (v i>x - XiVi)9'iWi - (lV i:X - XiWi)Vi9i 

i 

+ ^2 h,a U rj (v i>x - XiVijVjB'jWi/Vi - (w i>x - XiW^VjO'i 

y~]r it(Ta [-(vj, x - KvijOi^O'iWi/vi + (w i>x - XiWi)9'i9 iy 



+ 



+ y~] Vk {fj 



(A10) 



79 



WiTi 



v + O'^Wi/vifri^ + (w itt + (A 



fi - e'^Wi/v^fi 



~ ^2 X i\ ( Ui >* + (^^) X ~ 
i *■ 

+ ^r i f ii „ - Aito^Vj - Wi(wj - X*Vj) + (wiVj) 



fi + Vif i>v + 0-(tOi/ui)f ijO .j - ^2(wi, t + (AiWi) x - w i>xx ^ 



+ ^2h,v {Vi,x ~ ^iVi)w i:X + (w i:X - XiWi)v i>x + (wiVi(X* - 6i)) : 
i 

+ ^ - \Vi) (e' i (w i /v i ) 2 ) x - (w iyX - XiWi) (9 iyX + (6iWi/Vi) x ) 
i 

+ ^2 f i,vufi {v iyX ~ \iVi)wiVi + WiV 2 (A* - Oi) 
i 

+ (y itX - XiVi)wiV iiX + WiViV iiX (X* - Oi 

i 

+ ^2fi !VU fj (v i:X - XiViJwiVj + WiViVj(X* - Oi) 

+ (*V ~ ^i V i){ - W i e i,* + 0'i( W i/ V i) 2v i,x) ~ {Wi,x ~ XiW^O-Vi^Wi/Vi - WiVi(X* - 0j)0 i)a 

i 

+ ^n >CTU fi (w i)!r - XiViJOlwi/vi - (wi, x - XiW^O-Wi 

i 

+ ^2fi,aufj (V i>x - XiV^VjO-iWi/Vi) 2 - (W i>x - XiW^VjO-Wi/Vi 

+ ^fi,aa [-(v iiX - XiV^O-iwi/vifOi^ + (w i)a . - XiWi)0'O itX Wi/i 

i 

^ ^ WjVj ^ ^ ^fc (j"i,uu (fj ® fk) ~\~ fi,ufj,ufk) Vj,xTi,ufj,v ^i,xfi,uvfj ^j,xfi,ufj,a Oi,xfi,u<jf j 
i^j k 

+ ^{wi - X* Vi )vj (fj • A(u))fi - (n • A(u))fj 
- 5^ A*ai(t,^) 

i i 

(All) 

Recalling the expression (5.10) for the differential dA/d(v,w), we recognize that the equations 
(A.lO)-(A.ll) provide the explicit form of the system (6.5). The uniform invertibility of the differ- 
ential of A implies the estimates 

<t> j: 4, j = 0(l)-J2{\*i\ + \bi\). 
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To prove Lemma 6.1, it thus suffices to show that all the terms in the summations defining aj,6j 
have the correct order of magnitude. 

First of all, one checks that all those terms which involve a product of distinct components 
i 7^ j can be bounded as 

C> ( 1 ) ' ^2 (\ v .i Vk \ + \ v 3,x v k\ + \vjW k \ + \v jjX w k \ + \vjW k>x \ + \wjW k \). (A.12) 

In most cases, this estimate is straightforward. For the terms containing the factor 9i )X or 9j tX this 
is proved as follows. Recalling the bounds (4.24) we have, for example, 



0,, r,„ = 0(1) • v 3 of^'^^W* = 0(1) • (\ Wj J + |t, i>x |) = 0(1) • 5* , 
because of (5.24). Hence 

ViVj9j !X fi tU fj t(T = 0(1) • ^ok-Ujl • 

Next, we look at each one of the remaining terms on the right hand side of (A. 10) and (A. 11) 
and show that its size can be bounded as claimed by Lemma 6.1. To appreciate the following 
computations, one should keep in mind that: 



1. By (6.16) there holds 



Vi, x ~ (A; - \*)vi -Wi = 0(1) • 5 ^2 N • 



Therefore 



(H + K| + K:c| + K,x|) V itX -(Xi-\l)Vi-Wi = 0(1) -So^ {\ v i v j\ + \ w i v j\ + \ v i,xVj\ + \w i>x Vj\) . 

2. By (5.5) the cutoff functions satisfy 9\ = 9'/ = whenever \wi/vi\ > 35\. 

3. By (4.24) we have f^/vi , r ijaa /vi , f i)(7tt /«i = 0(1). 

4. One can have \wi — 9{Vi\ ^ only when \wi\ > S\\vi\. In this case, (6.18) yields 

v i = O(l)-v itX + O(l)-5 J2\ v j\- 

What follows is a list of the various terms, first those appearing in a^, then the ones in 6j. 
Coefficients of f ijU fi : 



Vi{v i>x - XiVi) - Vi(wi - X*Vi) = Vi [v i>x - (Xi - X*)vi - Wi] , 
Vi(w i>x - XiWi) - Wi(wi - X*Vi) = [viW i:X - v i>x Wi] + Wi[v i>x - (Xi - X*)vi - Wi] . 
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Coefficients of f 



i,v • 



2vi, x (v iiX - XiVi) + (v 2 (A* - 6»i)) x = 2u iiX [v i>x - (Aj - A*)^ - 0jVj] + 0- [t^to* - Vjifv] 

- (A, - X*)vi -Wi+ 2v i>x [wi - 6iVi] + 0- [vi, x «;i - ^w^] 



= 2-y. 

Wi,x( V i,x ~ + V iA W i,x ~ + (™i«i(A- - 0j)) a , 

= 2w i>x [v itX - (Aj - X*)vi - Wi] + 2w^ x [wi - 6iVi] + (A* - 0j — Aj + 9-Wi/vi) [v itX Wi - ViW^ a 



Coefficients of fi j(T /vi : 

Vi(v ijX - Kvi) (-0 ijX + (8'iWi/vi)^ - Vi(w itX - XiWi)e' i x ■- 

Vi(v itX -XiVi) (oKwi/vi) 2 ^ -Vi (w itX -XiWi) (9i, x + (9iWi/vi) x y- 



{viW ijX - WiV ijX )6"(wi/vi) x 

r i 2 

■0'i[vi{wi/vi) x \ , 
■{e , ^w i /v i )+29^ [vi {w i /v i ) x 



Coefficients of f iyVU fi : 



v 2 {v iyX - X^) + vf(X* - Oi) = v 2 v iyX - (Aj - X*)vi - Wi + v 2 [wi - 0j«j] , 
ViWi(v itX - XiVi) + v 2 Wi(X* - Oi) = ViWi v i)X - (A, - X*)vi - w t + ViWi [wi - 0j«j] . 



Coefficients of r 



l,VV ■ 



ViV i)X (v i)X - X^i) - v iyX v 2 (X* - Oi) = ViV i)X [v itX - (Aj - A*>j - w^ + ViV iyX [w t - 0^] , 
WiV i>x (v i>x - X^i) - v i:X ViWi(X* - Oi) = WiV i>x v i:X - (Aj - A*)t>j - w { + WiV i:X [w { - 0. 



Vi . 



W, 



20,' I Vi. x — - Wi., 

Vi 



X* + 0j)0- [w i)X Vi - WiV i)X \ , 



Coefficients of f j iW : 

(v i>x - XiVi)(- Vi9 i>x + 9-Vi^Wi/vi) - (w i:X - XiWi)9lv ijX - v 2 (X* - 0j)0 
= v ijX 9i(vi tX Wi - ViW ijX )/vi - 9 itX Vi(v i;X - (Aj - A* + 9i)v^j 

| [vi, x - (Aj - A*)«j - Wi] + [wi - 0jUj]| + (Aj - 

{vi, x - X^i) (-WiO i)X + 9i(wi/vi) 2 v i)X ) - (w i)X - XiWi)9'iV i)X Wi/vi - ViWi(X* - 9i)9 iy 
= 20-^ (vi,x^- - Wi, x ^j { [vi, x - (Aj - X*)vi - Wi] + [wi - 0jt> j]| 
- (Aj - A* + 9^9'^^ - WiV i)X ]wi/vi . 

Coefficients of f itCru ri/vi : 

(vi, x - XiVi)w i Vi9' i - {w i)X - XiWi)v 2 9'i = 9\vi [v iyX Wi - ViW i)X ] , 
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(vi, x ~ XiV^wfe- - (w iyX - XiW^WiViO- = O'jWi [v i)X Wi - ViW iyX ] . 



Coefficients of f i>crcr /vi : 

~{vi, x ~ Kv^WiQ'fii^ + (w iyX - XiWijviO'idi^ = -(0-) 2 Vi {wi/vi) x 

-(vi, x ~ Kv i )wi(wi/v i )9' i 9 itX + (w i>x - XiWijWiO'iOi^ -- 



Vi - 



This completes our analysis, showing that all terms in the summations that define have the 
correct order of magnitude, as claimed by Lemma 6.1. □ 



Appendix B 

We compute here the source terms fa, Y>j in the equations (11.15) for the components of a first 
order perturbation, and prove Lemma 11.4. We recall that 

z = ^2hiri(u,Vi,X* - 9(gi/hi)), T = - X*hi)fi(u,Vi, X* - 6(g i /h i )), (5.1) 



n = fi(u,Vi, X* - §i), Xi = (r i} A(u)fi). 



As in (A.4)-(A.ll), the computations are lengthy but straightforward: one has to rewrite the 
evolution equations for z and T : 

z t + (A(u)z) x - z xx = (u x • A(u))z - (z • A(u))u x , 



T t + (A(u)T) x - T xx = (u x • A(u))z - (z • A(u))u x - A(u) (u x • A(u))z - (z • A(u))u, 



+ (u x • A(u))T — (u t • A(uj)z . 



in terms of hi,gi. The fundamental relation (4.23) implies 

A(u)fi = Xifi + Vi(r iiU fi + (Xi - X* + Oi)r itV ) . 
Differentiating (B.l) w.r.t. x and using (B.3) we obtain 

z x - A(u)z = ^ h i,x?i + ^2 h ^ifi^fi - A(u)r^ + ^2hi \v itX ri, v - (0'i(hig iiX - gih itX )/h 2 )f it 

i i i 

-\- ^ ^ hiV jTi u T j 

= ^2{hi, x ~ KhijVi + ^2h i v i r ijU (f i - r*) +^/i»[t>i, x - (K - X* + §i)vi f i>v 



(5.2) 
(5.3) 



(gi,x - Aiflfj) - j^{hi,x - Xhi) 



fi,a ^ ^ hiVjVi jU Vj, 



= ^2{hi, x - Xihi) h + 0'i(gi/hi)ri^ - J2(gi, x - kg^d'A,* 

i i 
+ ^ h iVifi,u{h - h) + ^2 ki [ Vi > x ~ ^ ~ X i + ^ v + ^2 h i V 3 r i,nfj , 



(5.4) 
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- A(u)T = ^2(gi, x - X*hi tX )fi + Y( 9i - K h i) *W, u fj - A(u)r^\ 

i i 



Y( 9i > x - Kgijh - ^ A* (h i>x - Xih^n + - A*^ i )-y i r ii „(f i - 



J, I) 



A* 



{gi,x — Aift) — —(hi,x ~ \hij 



^» - 8'i(gi/ h i)fi,, 



+ ^ ft [«i,x - (A — A* + 0;)^ ri,„ + ^2 9iVih,u {h - h) + 9iVjr ijU rj 

i i i^j 



Vi, x ~ (Aj - A* + 0;)^ 



3& 



(5.5) 



Next, differentiating (B.l) w.r.t. t we obtain 

i i 

= ^2h i>t ri + ^2hi Vijh,v- {0i(9i,thi- gih iit )/hf)ri,v +^2 h i{ w j ~ XjV^ri^fj ( 56 ) 

i i J i,j 

i i i i,j 

i i 

= ~ x *i h i,t)fi + Y( 9i ~ X *i hi ) V i^i, v ~ {0'i{gi,thi - 9ih iit )/hfjfi,. 

i i 

i i i i,j 

- Y K^ h i,t[^ + ^(ftAi)^] - Q'i9i,th,a + hiV i)t f iyV + Y h i( w j ~ ^j v j)n,ufj^- 



Differentiating again z x — A(u)z and — A(u)T w.r.t. x, from (B.4) and (B.5) one finds 



{B.l) 
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i i 

+ {-o l , x + {e' i g i /h i ) x )f i ,a 

i L j 

+ d'iV i)X {gi/hi)r iyr7V + ^2 / v j 9' i {gi/h i )r i ,uafj - Oi,J'i(gi/hi)r ija 

j 

+ J2(9i, x - K9i) -di{9i/hi) x fi,* -^^VjTi^fj - Vi^r^a + 0-0i, x ri, c 
i L j 

+ ^{hiVifi,u{fi - fi)) x + ^2(hi(vi, x ~ (Aj - A* + 0i)uj)^ h,v 

i i 
^ ^ (Aj Aj + 6{)Vj^ ^ ^ Vjfi vu fj + r Vi,x^i,vv @i, X T'i,v 



i,vTj,vTk + ?i,uu(Tj <S> ^fc)) 

i/j i/j L fc 



"I - "V j )X T i ^ u f j >v -\- Vi^ x fi^ vu rj 6 'j x f "j u f 'j j(T 0i^ x Ti tO - u f j 



(B.8) 



(A(u)T) x = ^2(h i:XX - (\ihi) x ) 0-(gi/hi) 2 r i:a \ +^(g itXX - (Xigi) x ) h - d'i(9i/hi)h,a 



+ ^2v j 9' i (g i /h i ) 2 f ij(TU r j - 9i(gi/hi) 2 §i tX fi 



y^^jh,ufj +v itX f ijV - (9 ijX + (6'igi/hi) x )f ita 

i L j 

-^2 l v A{9i/h i )ri,crufj - v^JKgi/hiYi^ + h, x K^9il^r^ a 
j 

+ ^2[gi{vi, x - (Ai - A* + 0j)vj)j f itV + ^2^giVif ijU (fi - fifj 

i i 

+ ^Zdi{vi, x - (Ai - A* + Oi)Vi) 'Y^v j f i ,vfj + v i:X f i>vv - i:X f i:VO . + ^2(giVj) x ri, u rj 

i 1 j J i^j 



+ hi(v i>x - (Ai - A* + 9i)vi)r i)V + h^f^^i - h) + ^v i v j r i>u f j I . (5.9) 

j/i > * 
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Substituting the expressions (B.6)-(B.9) inside (B.2) we obtain an implicit system of 2n scalar 
equations governing the evolution of the components hi,gf 

^2(hi,t + {\hi) x - h itXX ^j ri + Q'iigi/hiYi,^ + ft i9i ) x ~ 9i > xx ) ~ ^ fi > CT 

i i 

= ^2r itU ri[vi(hi >x - Khi) - h^Wi - A*) 

i 

+ ^2h, u rj (h i:X - Xihi)vj + (hiVj) x - h^wj - \*Vj) 
+ ^2?i,v {hi, x ~ \hi)v itX + (hi(v itX - (Aj - A* + §i)vif) -hi 

i L 

+ 5Z fi > CT ( hi > x ~ ^i k i)(~ + (®'i9i/ h i)x) - {9i,x ~ K9i)9'i, a 



hiVi(v ijX - (Aj - A* + 0*)^) 
M j {vi, x ~ (Aj - A* + 0*)^) 
hiV i!X (vi !X - (Aj - A* + 0;H) 



+ X^>« CT ~ Khi)9iVi !X gi/hi - (g i;X - Xg^Vi^ - hi(v ijX - (Aj - A* + 9i)vi)0 ija 



+ ^2r ijau fi (h ijX - Xih^ViO-gi/hi - (g ijX - Kg^vfy 

i 

+ ^2h,a U rj (h ijX - Kh^VjO^gi/hi - [g iyX - Kg^vfy 

+ ^2 ^> aa [~ ( ki > x ~ ^i h i)d'A,xgi/ h i + {9i,x - \gi)6'iKx\ +^2{hiVifi tU (ri - fi)j 

i i 

+ ^ [(r,- • ^0))rj - (fj • ^(u))f 



^ai(t,x), 



(B.10) 
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+ Q'i(9i/hi)ri,a\ + ^2(gi,t + (\9i)x ~ 9i,xx^j 



i 

^2 ri,uh [(gi, x - \9i)vi - 9i(wi - X*Vi) 

i 

+ ^2 ^> u ^' ( 9i > x ~ X i9i)vj - 9i{wj - X*Vj) + (giVj) 
+ ^2^i,v (gi{vi, x - (A; - A* +0i)vi)^j + (gi, x - \9i)vi, x - gwi, 

i L 

+ X^> CT ~ ^i h i){9i(.9i/hi) 2 ) x - {gi, x - \9i){0i,x + (0'i9i/hi) x ) 

i 

+ ^ fi,vuh (v i>x - (Aj - A* + 9i)Vi)gi' 

i 

+ ^hw (vi, x - (A, - A* + 9i)vi)giV i} 

i 

+ ^ n^f-,- (v i)X - (Aj - A* + 9i)vi)gi 

+ ^2h,va (hi, x - Khi)0'i(9i/hi) 2 Vi tX - (g i>x - Xg^Vi^gilK - gi(v i:X - (A, - A* +^)^)^ 

i 

+ ^2 / r itau fi (h ijX - \ihi)0' i {g i /hi) 2 Vi - (g ijX - Xigi)0' i v i gi/h i 

i 

+ ^2ri,a U rj {h iyX - XihijvjS-igi/hi) 2 - (g ijX - Kg^vfygi/hi 
+ ^2h,aa \~{hi, x - \ihi)9 l i {g i /h i ) 2 9^ x + (g itX - Kgi)9' itX gi/h^ + ^2{giV i r i ^(f i - r*)) 

i i 

~l~ ^ ^ g%^j ^ ^ Vk i j ,-uT k ~t~ ^i,uu{Xj ® ^fc)) ~\~ V j t x^i,uT j,v Vi tX Ti jVu fj 0j,xfi,u^j,cr 9i,x^i,ouT j 
i^j 1 k 

+ ^{wihj - Vigj) {fi • A(u))fj + ^2 Vihj(jfi • A(u))fj - (fj • A(u))r^j 
- J^AJ Oi(t,x) 

i i 

(£•11) 
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Recalling the expression (11.11) for the differential dA/d(h,g), we can write (B.lO)-(B.ll) in 
the more compact form 



dA / [h iit + {\hi) x - h iyXX ] 



E 



a. 



d(h,g) V [gi,t + {\9i) x ~ 9i,xx] J Y^'" X * 1 ^ 

By the uniform invertibility of the differential of A, to prove the estimates stated in Lemma 11.4, 
it suffices to show that, for every i = 1, . . . , n, the four quantities 

can all be bounded according to the right hand side of (11.16). 

We start by looking at all the terms in the expressions (B.lO)-(B.ll) for ctj and 6j. First of 
all, one checks that all those terms which involve a product of distinct components % ^ j can be 
bounded as 

For convenience, quantities whose size is bounded as in (B.12) will be called "transversal terms". 
More generally, quantities whose size is bounded according to the right hand side of (11.16) will 
be called "admissible terms" . We denote by A the family of all admissible terms. We now exhibit 
various additional terms which are admissible. 

1. By (6.16) it follows 

{\hi\ + \gi\+\hi,x\ + \9i,x\)\v itX - (A; - X*)vi - ^ 

= 0(l)-So^2(\hiV j \ + \g i v j \ + \h itX v j \ + \g i:X v j \) e A. 



2. Two other other admissible terms are 



hi[w itX Vi - WiV iiX ] = [hiW itX - Wih itX ]vi + Wi[h itX Vi - hiV iiX ] £ A, 
9i[wi, x Vi - WiV i:X ] = [giW i>x - Wig i:X ]vi + Wi[g i:X Vi - giV i>x ] G A. 



(£.14) 



3. We now consider terms that involve the difference between the speeds: 9i — Qi. We claim that 
the following four quantities are admissible: 

h i v i (9 i -9 i ), giVi(6i-9i), h iyX Vi(9i- 9i), g itX Vi{9i - 6i) e A . {Bib) 
Indeed, from the definitions and the bounds (4.24) it follows 

|Ai - Ai| = 0(1) • \h - h\ = 0(1) • Vi \§i - 0i\ = 0{l) ■ 6 \0i - 0i\. (B.16) 
Since \9'\ < 1, one has 

\9i-9i\ < \(gi/hi)-(wi/vi)\. 



Using (6.16) and (11.12) we now obtain 



\hiVi 



< \9iVi -Wihi\ 

= (hi,x + C\ ~ K)hi + Oil) ■ S J2 {\hj\ + \vj\))vi 



(«i,x + (A< " K) V i + O(l) • 5 J] K-l)^ 



- fi.x^i) + (Aj - \)vihi + 0(1) • <5 X + \hjVi\) 



3+i 



Hence 



< \hi, x Vi - v i)X hi\ + 0(1) • 5 \§i - 9i\\hiVi\ + 0(1) • <5 ^ (|ujt>fc| + l^fcl) • 



IflfjUi - Wihi\ < 2\h iyX Vi - v iyX hi \ + 0(1) • 5 ^2 {\vjV k \ + \hjV k \) € .A, (-B.17) 



showing that the quantity hiV^Oi — Oi) is admissible. 

Observing that Oi — Oi ^ only if either |<7i/foj| < 65i and |iUi/t>j| < 3<5i , or else \gi/hi\ > 6<5i 
and \wi/vi\ < 3<5i, we can write 



9iVi(0i 



< h/hl \hMOi-0i)\ -X{ |lB/fci| < Ml } +^\g iVi \ -X { ^ /h ^ | WiM |< Ml } 

< 6(5i|ft//ii| |/iifi(6*i - 6*^ ) | + id^giVi - Wihi\. 



Hence giVi{9i — Oi) € A. In turn, using (11.12) we obtain 

hi, x Vi(§i - Oi) = g t v i (§ i - 6i) + (Ai - X^hiV^i - Oi) + 0(1) • <5 ^ (\v iVj \ + |^-|), 

showing that the term hi tX Vi{9i — Oi) is also admissible. Finally, using (6.16) one can write 



9i, x Vi{9i - 9i) = (9i - Oi) [gi, x Vi - v i>x gi] + giV i:X (9i - Oi) 
= {9i - Oi) [gi, x Vi - v iyX gi] + giWi(9i - t ) 

+ CXi-X*)gM0i-0i)+O(l)-S J2 



ViV 



3 I • 



To estimate the term giWi(9i — Oi), we observe that Oi — Oi = if \wi/vi\ and \gi/hi\ are both > 3<5i. 
Hence, using again (6.16), we can write 

\9iWiipi -0i)\= M^giVitfi -6i)\- X{ |u>lM|<35l } + 35 1 \h l w l (9 l - 9 % )\ ■ X{ |fli/hi|<35i } 
= 3<5i|^Vi||^ - 9i\ + \hi(v i>x - (K ~ K)vi) & - 0i\ + O(l) • \ hiV i\ ■ 

3& 
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By the previous estimates, this shows that gi :X Vi(9i — 0i) € A, completing the proof of (B.15). By 
(B.16), the following terms are also admissible: 



hi(Xi — Xi), gi{Xi — Xi), hi tX (Xi — Xi), gi, x (Xi — Xi) £ A. 



(B.18) 



4. Next, we claim that 

hi{fi-fi) x , gi(ri-fi)x, hi(\i-\i) x , ^(Aj - A;) x e A. 
Indeed, one can write 



(£.19) 



hi(fi - fi) x = hiVi(9i - 9i) 



EVi.u TijUjTj 



i.v 'i.v 



Vi(6i-6i) ' Vi(0i-0i) 



+ 9' i {g i /h i )[v itX h i - h ijX Vi](fi^/vi) + - giVi iX ](ri ya /vi) 

+ (wi/v^e'i [v i)X hi - Vih i)X ] (h,cr/vi) + 9'i [h i)X Wi - hiW iyX ] (f i)0 ./uj), 



9i{fi-Ti)x = 9iVi0i - 9i)\j2 v j ^Zl 3 + Vi ' x ~ 2 7k 

I j Vi(9i - 9i) Vi{9i 

+ 9' i (g i /hi) 2 [vi tX h i - h iyX Vi](r iyCr /vi) + 9' i (g i /hi)[v i g^ x - giV i)X ] (r iyCr /vi) 
+ (w i /v i )9' i [v ijX gi - Vig ijX ](fi jrT /vi) + 0- [g itX Wi - giW ijX ](f ita /vi). 

By (4.24), the above expressions within braces are uniformly bounded. Hence the first two quan- 
tities in (B.19) are admissible. To prove the admissibility of the last two terms it suffices to repeat 
the above computation, with fj and fi replaced by by Xi and Aj. 

In a similar way as in Appendix A, we are now ready to check one by one all the (non- 
tranversal) terms in the expressions of Sj, 6j in (B.lO)-(B.ll), showing that all of them are admis- 
sible. 

Coefficients of r itU fi : 

Vi(h i>x -Xihi) - hi{wi - X*v{) 

= [vih i)X - hiV i)X ] + [vihi(Xi - X^] + hi(v iyX - (A; - A*)^ - Wi) 



Vi{gi, x -Xi9i) ~ 9i(wi - X*Vi) 

= [gi, x Vi - v i:X gi] + [vigi(Xi - A,)] + gi(v i:X - (Xi - X*)vi - 



Coefficients of f 



v i>x {h iiX - Xihi)+(hi(v itX - (Xi - X* +9i)vi)^ - 



hi, x (v itX - (K - X*)vi - Wi) +2[h i;X (wi - 9iVi)] 

+ [hiV i:X - h i>x Vi] (A* -§i-Xi + 9-gi/hi) + 0- [vig i>x - v i:X gi 

+ 2[h i:X Vi(Xi + 9i - Xi - 9i)] + hi{(Xi - Xi)vi) x - M>» , 
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(gi{vi, x - (A; - A* + + v ijX (g itX - Kgi) - giv it t 



9i,x{vi, x - (Aj - A*)«j - tOj) +2 g i;X (wi - 0^) 
+ (A* -9i-Xi + d'i9i/hi) [v i:X gi - Vig i:X ] + 9'^gi/hi) 2 [vih i>x - v i>x hi] 



+ 2[g i>x Vi(Xi + 6 i -X i - §i)] + gi{(Xi - Xi)vi) 

Coefficients of r i>cr /vi : 

Vi(h i:X - \ihi)(-§i, x + (6'igi/hi) x ) - Vi(g i:X - Xigi)9' ix = - 



-9i 



v. 



(hi, x - Xihi) ^0' i (g i /h i ) 2 ^ + Vi(g itX - A^;) (§' i x - (0-ft//i;) J = -1 



hi 



Coefficients of r ? - „„f ? - 



ihi(v itX - (Aj - A* + = Vihi(v itX - (A, - A*)t>; - Wj) + Vihi(wi - QiVi) 



+ 



^^(Ai + ^-Ai-^) 



- (Aj - A* + = Vigi(v i:X - (Aj - A*)^ - to*) + - 0^. 



+ 



vfgi&i + Qi - K-Oi) 



Coefficients of f 



hiV ijX (v ijX - (A; - A* + 6>j)^) = hiV ijX (v ijX - (A; - A* + 0;)^) + hiViV ijX (Xi + $i - Aj - 0*) 

= Vih i>x (v i:X - (A; - A*)t>; - lOj) + hiViVi^Xi + 0; - A; - 0^ 

+ - 0jWj)] + (u i>a; + (Aj - A*)vj - 0jVi) [hiv i>x - Vihi i3 

giVi, x (v i>x - (A; - X*)vi + 0j-yj) = Vi gi, x (v i>x - (A, - A*>; - w;;) + [swuj^A; + 0; - A; - 0;)] 

+ Uj [gi,x( W i ~ °i V i)] + { V i,x ~ (Ai - K) V i ~ °i v i) [flW,* ~ 

Coefficients of fi jV(T : 

(hi, x - XihijO'^gi/h^Vi^ - (g i;X - A i ft)0^ iiX - hi(v itX - (A, - A* + 9i)vi)6 i;X 
= Vi^O'^hi^gi - hig ijX )/hi - 9 ijX hi(v itX - (A, - A* + 0;)^) 
= 29' i (h itX (g i /hi) - g iiX ) {v i>x - (A; - X*)vi - Wi) + (wi - 0^) 

+ (2(\ - X* + 9i) - (Ai - A* + 0;)) [v iiX gi - Vig itX ] - 0-(&//ii) [v itX hi - Vihi 
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(hi, x - Xih^e'iigi/hifvi^ - (g i>x - Xi9i)0'iVi yX gi/hi - gi(v i:X - (A; - A* + 0i)Vi)0 ita 



na'9j_ I i 9_ 



(vi, x - (Aj - X*)vi - Wi) + (wi - Ovi) 



+ (2& - A* + 00 - (A, - A* + 0-i) [e'iigi/hi) 



x^i Villi x 



Coefficients of f itCru fi/vi : 

(h ijX - Xhi)fyv?gi/hi - {g i>x - A;ftK ? 0- = Ofa [v i>x gi - Vig ijX ] + K v A9ilK) [vih ijX - h iV i tX \ 

(hi, a - \hi)9'ivf{gi/hi) 2 - (g iiX - X l g l )9' i v 2 (g l /h l ) 

= 9'iVi(gi/hi)^[vi tX gi - Vig itX ] + (ftMO [vih i>x - j. 

Coefficients of r itCrcr /vi : 

~(hi, x - hhijv^gi/hijO'iOi^ + (g ijX - Xg^v^'fii^ = (0'i) 2 Vihi[(gi/hi) x ] 2 , 

~{hi, x - Xih^Vitgi/hifO'fii^ + (g i>x - A^O^z^ftMO = (^) 2y iPi [(ft AO*] • 



There are a few remaining terms in (B.lO)-(B.ll) which we now examine. Recalling (B.14) 
we have 

(hiVifi lV (fi - f l )) x =0(1) • h l)X v 2 {0i - 0i) + 0(1) • hiViVi^Oi - §i) 

+ 0(1) ■ hivUOi - 9i) + 0(1) ■ hiViifi - fi) x , 

(giVif i>u (fi - f0) =0(1) • gi,xV 2 (0i - §i) + 0(1) • giViV itX (6i - §i) 

+ 0(1) • 9i vU0i - 9i) + 0(1) • giViifi - h)x , 



hiVi (fi»A(u))fi-(fi»A{u))fi = 0(1) ■ hiViifi - fi) , 

hiViA(u) (fi»A(u))ri-(fi»A{u))fi = 0(1) ■ Kv^fi - fi) . 

(wihi - Vigi){fi • A(u))fi = 0(1) • \wi~hi - g^ , 



Vihi 



{{fi* A{u))fi- (fi»A(u))f^ 



0(1) • {\vi, x hi\ + \vih ijX \)\fi - n\ + 0(1) • Vihiiji - fi) x , 



These terms are all admissible because of (B.15)-(B.19). 

We have thus completed the analysis of all terms in (B.lO)-(B.ll), showing that the quantities 
Sj, hi are admissible. The admissibility of the terms ((Aj — AO^i) and ((Aj — A0<7t) follows 
immediately from (B.18) and (B.19). This completes the proof of Lemma 11.4. □ 
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Appendix C 



Aim of this section is to derive energy estimates for the components hi, gi and prove the 
bounds (11.33)-(11.34). We write the evolution equations (11.15) for the components hi,gi in the 
form 

{hi t t + \^ihi)x ~ hi jXX = (f>i , 
(CI) 
9i,t + {Kgi) x - 9i,xx = A ■ 

For convenience, we define fji = rj(gi/hi). Multiplying the first equation in (C.l) by hifji and 
integrating by parts, we obtain 



J rjihicpi dx = J {fjihih iit + f)ihi(Xihi) x - fnhih iiXX } dx 

= J {vi{ h2 i/ 2 )t - r\i\KK,x ~ fii,x\hl + r)ih\ x + f) itX h itX hi} dx 

= J {{Vihi/ 2 ) t + Mi)x(hi/2) - (ih,t + 2~KVi,x - fn, xx )(h 2 /2) + f)ihl x } dx 

J rjih 2 /2dx + J (fii t t + KVi,x-Vi,xx)(h 2 /2)dx 



Therefore 



dx = — 

dt 



\,xVi(h 2 /2)dx+ I fjihifadx. 
As in (9.14), a direct computation yields 



(C.2) 



Vi,t + Krii,x - m 



xx Vi 



-t ( i>i 9i 4>i \ . hi,x ( 9 



hi hi hi 



+ 2fi[ 



hi V h 



(C.3) 



Since Aj ;X = (Aj — X*) x , integrating by parts and using the second estimate in (11.13) one obtains 



\ x fn{h 2 /2)dx 



fti ~ K){Vi,xh 2 /2 + fjihihi^) dx 



< ||Aj — A*||l°° 



^ J \Vi\\9i,xhi - gih itX \dx + ^- J fiih 2 ix dx 
< J \9i,xK - h itX gi\dx + ^ J f)ihl x dx + 5 Q ^2 J (\hiVj\ + \hihj\) dx , 



because 

\X i -X*\=O(l)-S «S 1 <l. 
Using (C.3) and (C.4) in (C.2), we now obtain 



(CA) 



\l^ h lxdx<-j t 



/¥-B/ 



+ 2 / \m{\hM + \gi<t>i\)dx + 



Vihih iiX ( j- 



dx 



+ \ J (f 1 ) dx + J \9i,xhi- 9ih iiX \dx 
+ So ^2 J (\ h i v j\ + \hihj\) dx + J \hi<t>i\ dx. 



(C.5) 
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Recalling the definition of fji, on regions where f% 7^ one has \gi/hi\ < 4<5i/5, hence the bounds 
(11.14) hold. In turn, they imply 



< 



2«5i 



~, h 2 f 9i 



+ 0(l)-*o£ 



. / Si 



= O(l) • Ift^/tj - gih itX \ + O(l) • 5 ^2 (\ v i9i,x\ + K'^i,z| + \hjgi, x \ + |^i,z|) • 

(C.6) 

Using (C.6) and then the bounds (11.18), (11.26), (11.28), (11.30) and (11.31), from (C.5) we 
conclude 

J j ' Vi h lx dxdt < J f)ih*(t,x) dx + 0(1) • J (|^| + ctecft 

+ It J \ 9i ' xhi ~ 9ihi,x\ dxdt + 0(1) ■ J ' \hi(gi/hi) x \ 2 dxdt 

+ 0(1) ■ 5 ^2 J J (\vjgi jX \ + \vjhi tX \ + \hjgi !X \ + \hjhi tX \)dxdt 

+ $o^2 [ f (\hiVj\ + \hihj\) dxdt + 2 [ [ \h,fa\dxdt 

■ Ji J Ji J 

= 0(l)-dg, 
proving the estimate (11.33) 



(C.7) 



We now perform a similar computation for g? . Multiplying the second equation in (C.l) by 
rjigi and integrating by parts, one obtains 

J fugiipi dx = J {{rh9i/ 2 ) t + (X^i) ^gf / 2) - (r) i>t + 2Xif] i)X - f] i)XX ) (gf/2) + r\ig\ x } dx . 

Therefore, the identity (C.2) still holds, with hi, fa replaced by gi,fa, respectively: 

d 



Vi9i, x dx = - 



dt 



f] i gj2dx 



+ / {Vi,t + KVi,x ~ r)i,xx) (Si / 2 ) dx 



J KxViiSi/ 2 ) dx + j fug^i dx. 



(C.8) 



The equality (C.3) can again be used. To obtain a suitable replacement for (C.4) we observe that, 
if fji 7^ 0, then (11.13) implies 

|SiSi,x| < 2\h i>x g i:X \ + O(l) • 5 ^2 (\ v j9i,x\ + l^jSi,*!) 

J#i 



and hence 



|SiSi,x| < h\ x + g 2 ix + O(l) • 6 {\ v j9i\ + l%i|) 
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Integrating by parts we thus obtain 

J K, x f}i{g 2 i/ 2 ) dx = J (A; - X*)(fj itX gf/2 + fjigig i;X ) dx 



Wi \ \9i,xhi - 9ihi 



9i_ 

h' 2 



dx + f] i h 2 ix dx 



+ / Vi9i, x dx + ■ <*o ^2 / (\ v j9i\ + \hj9i\) dx 

< J \9i,xhi- h itX gi\dx+^ J f) i h 2 ix dx + ]^ j ^gf^dx + j + dx . 

(C.9) 

Using (C.3) and (C.9) in (C.8) and observing that \gtlh\\ < 5f on the region where fj^ 7^ 0, we 
now obtain an estimate similar to (C.5): 



fa w 2 x dx<- 



d_ 

dt 



Vi9t 



dx 



5 2 

+ -i I \vtt{\hdi\ + \9i$i\) dx + 6f 



fj'ihihi 



dx 



4 



dx + J \gi, x hi - gih iyX \ dx + ^ J f\ih\ x dx 
+ <5 ^y (\v j g i \ + \h j g i \)dx + J \gSi\dx. 



(CIO) 

Using (C.6) and then the bounds (C.7), (11.18), (11.26), (11.28), (11.30) and (11.31), from (CIO) 
we conclude 




Vi9i,xdxdt< J fj i g i (t,x)dx + 0(l) 

T 




{\hii>i\ + \gi4>i\) dxdt 



+ 0(1)-^ J \gi,xhi- gih i;X \dxdt + 0(l)- J J [h^gi/hi)^ 2 dxdt 
+ C(1) ■ S ^2 J J {\vj9i, x \ + \vjhi, x \ + \hjgi tX \ + \hjhi, x \)dxdt 



(Gil) 



= 0(l)-5 2 , 
proving the estimate (11.34). 
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Appendix D 

We derive here the two estimates (12.9)-(12.10), used in the proof of Lemma 12.1. 



PIU-^y \G x (t - s,x - y)\E(s,y)dyds 

= \\Ml~ f I |X ~ yl B{s) 

" " Jo J 4(t-s)y/ir(t-s) 



' exp {~ 4(t-^) + A W DA Wl™ J q \\ u x(<r)\\ L ood<T + a - y^dyds 
< ||A|| L =oexp|4||DA|| L o= jf \\u x (a)\\ Loo da + t - xj 



7 

Jo 



B(s) ( f, . f (y + 2(t - s) - x) 2 



= exp [qDAU- j)\u x (a)\\ Loo da + * " *} jT ^==^ (/ K " v^e"^) 
<exp|4||I>A|| L =o jf \\u x (a)\\ Loa da + t-x^ p|| L ~ {^= g + V*) B ( s ) ds 
< exp |4||I>A|| L » J \\u x (a)\\ Loo da + t - 
= ±E(t,x)-±e X p{4\\DA\\ L ~ J 



x 
t 



B(t) 1 



\u x (a)\\ da + t - x 



2\\DA\\ LOO J \\ux(s)\\ Loo (^J G(t-8,x-y)E(s,y)dyj ds 

= 2\\DA\\ Laa j B( 8 )-exv^\\DA\\ Lao j \\u x (a)\\ L00 da + . 



\u x (s) 



2y/n(t-s) 



exp 



{x - yf 
A(t-s) 



y\dy)ds 



<B{t)e t ~ x 2||DA|| L »||tt x (s)|| Loo exp|4||DA|| L o= jT \\u x (a)\\ Loo da^ ds 
= Bifye*-* \ exp |4||DA|| Lo = j\u x (a)\\ L00 da\ - \ 



<lE(t,x)-±e<-*. 
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